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ABSTRACT 


Future fleet needs for real-time intelligence require an unmanned 
platform capable of operations from small surface combatants 
without the need for extensive Support personnel or equipment and 
without causing disruption to the operations of the ship from which 
it operates. A candidate must not only takeoff and land vertically 
but also be capable of high forward flight speeds and efficient on- 
station performance. The design and initial fabrication of a Vertical 
Attitude Takeoff and Landing (VATOL) Unmanned Air Vehicle (UAV) 
airframe was completed at the Naval Postgraduate School. The 
vehicle, called Archytas, was a combination of two existing UAV's-- 
the AROD and Aquila--as well as locally manufactured components, 
including a canard support structure and wing spar. The objective of 
creating Archytas was to provide a_ proof-of-concept platform for 
research to explore performance trade-offs and_ stability 
augmentation. A_ three-degree-of-freedom simulation was used as 
the focus of the design efforts, to validate design decisions made in 
the fields of propulsion, aerodynamics, structures and flight 
mechanics. Engine tests were conducted to determine thrust and 
contro] power. Structural components were designed, fabricated and 
then tested, making modifications where necessary to ensure 
Sufficient airframe strength. A longitudinal control system was 


designed, validated by simulation, and tested structurally. 


11] 


ACKNOWLEDGEMENTS 


This work is essentially the result of life-long advice, instruction, 
nurturing and motivation from people too many to thank. I am 
certainly thankful to my parents, William and Nancy. My father 
showed me by example the skills and thought processes of an 
“Ingegnere’. I could have no better example of an engineer, or a 
father. My mother raised me to believe in myself--I draw from her 
nurturing every day of my life. And to my other mom, Lenore, 
who's insistence on education as a means to an end has not been lost 
on me. 

My work at NPS would not have been possible without the strong 
support of Professor Rick Howard, a gifted mentor who exemplifies 
the type of teacher we all need....we all want to be. 

But above all I am indebted to my family: Anne, William and 
Allison for the daily support, patience and motivation without which 
this thesis would not have been possible. They reminded me of what 


the most important thing in life is: love of family. 


DUDLEY IKNOX LIBRARY 
NAVAL POSTGRADUATE SCHOO! 
MONTEREY CA 93943-5101 


TABLE OF CONTENTS 


PRR © MO Serttttre sa s.creecaccecrsrests ccoecdionsecsscecsinreesnecsecesecesiniasesssenececancensnerosrovscess ] 
Ok hens) MOUS ISD 6) Os aaa asap a Ra soe eR ee ] 
Peeve PIONS OF CURRENT SYSTEMS ................cccccscsessconnessssceseccenseess 2 
BN PMI AMINO) rere ric cno..nolete cease ca ntteccsccosettth tescusctesderieasensssastsansesascasscscatessecsesoese 3 

EMR NM IN errr eters tein c ee cc see tree oie sate cncsanseddeevdcsensdeactsanceaecedasvsscescasceneceonsevacecorne 5 
A. VERTICAL ATTITUDE TAKEOFF AND LANDING (VATOL)............ 5 

OT tM COM scl TTT ere c cae cseccsiscvces ote tccccasccucsdusssetssscecarescesers ees encecee: 5 
Ee Pm Dem (COMM) CEST OTS ctecs.s...10-Seecrreres.ciGthae es sccnseceeseccesersncsaseerever. 7 
B. NAVAL POSTGRADUATE SCHOOL UAV PROGRAMG .............e eee 12 
il, Gi@ ICTR sgl ee sa ees cos SEE ian tee Le 
2 SENG TT (ia li aces ra | ern Men. IS 

Seema PAS DESIGN--GENERAL.............:.secsscocseneccsocssesvseteeenssesecseeseesaccssesenes led 

DEES TICGISY 10) 61) 5 oi BIg 6S. algae se iy 
MPMEMITCMIRZ IOC rr et tere crete Tere sevteetcscetrsascie ua rostossses sv sseiéinecsererrssecstetitvtesseacsncssees 17 
SDRC TCG Aas 7 

GTN LIVE Jerre ce ttt tess cccccccccessutdcessscstecsssosssssssntecssetsstescossnsessesneatonterts slats: 

C DESIGN PROCESS--ARCHYTAS ATIRFRAME........ccccccccscsccsesscsscrsereeeees a0 

eer OMO RGANIZATION... co.cc. cccccceccsccsssssssvsrsessscsvsssevtecssssssssssilecusenssesevaed 2 

meres  APERIMENTS 0.0... .ccccccsscccssssscessscescsesssssscssssscssosssensosessssscssevevscereeres O38 
oo PUR CIS. © TCU Bt SRS DiS IRN a3 
Se SOR TON OP TEST EQUIPMENT viitrcccccccccccccscsstssssssscssscsssssssesseee 24 

PE MenemEeSCTIPtION and inStallatiOM Aeneseccccieicecsescsescsoeee 24 
EP OUO MOLE ANG COMITO] “EQUIPTIEMt 2.2.0... cccecccecscccecesseceececeocoeecnessosees pei 


3. Test Stand ......scccccccccccocossosscssececeeece crcceeeeee nicer antenna ean 25 


C SCOPE AND METHOD OP TESTS (i 2a 
D. RESULTS AND DISCUSSION ...000.00..00sccccccccssesstesseeieseee 2a 
LD. Geme ral .occc.cccssccceecsseseseescnescunscaecesessscsoeeeeteeenene sates ote Fy 

2. Thrust measurements secs cccsi eee 23 
3. Vane control power measurements... 29 

A, Geme ral cccccccccccccceccstseeeccssceecussceecseeneyee stan el Seana rr ys 

b. Dual-vane roll-control power...........2.0.......-) 32 

c. Four-vane roll-control power... 32 

d. Pitch control power modelling..............).e 34 

4. Effect on thrust from blockage in AROD duct 3). =e 33 

5. Impact of propwash on vane position....c2 a... ee bps 

A. NOISC...ccccccssceccenssiocconnsonnssosecnecssesocosgenesenese0sdeeteet = aa 43 

B. Fuel f1ow....0...........00cces000e00s00000.00000 44 

c. Acceleration/deceleration characteristi¢se.........——— 45 

Gd. Rehiability.........0:.....0:csssssaersessconnven ete nen 45 

V. WING SPAR DESIGN, FABRICATION AND TES]... oe 46 
A. GENERAL 1 cccccccsscscssoseeesssvsscesesassaocassvsaitteciscn2stnn eteeeee =r 46 
B. WING SPAR DESIGN AND FABRICATION .............e 46 
L. Design 10a ....ccccccccccssscostssesescsssce eeeesseees eee er 46 

2. Forward spar design and attachment..........)eeeen 48 

3. Aft spar design and attachment..............-.0:mmly===—nn 50 

4. Integration of NPS spars with Aquila wing spars ee By 

C WING SPAR EXPERIMENT Sitea.c..........ccccnyeee 5 
1. Initial component testing..............0.1..02--2...)seeenes 52 


PCACSICMMOImbASICISpal SUTUCTURCMiauadelccanwetenut tet... osewns... 216 
CPCAT MMO MMTOGILICd SP AL... tees cco.s.c...sccccsnuMteereseessceedeeceasesseteessee.. SL, 
POMOC Nem Sihye SUT CLUIAL | LOS U .,;see-cleiis.e.c..cceo0tssdacsesdedstescdecsescsendiues onl 
vineANARD DESIGN, FABRICATION AND TEST cccciiiiicicccecescesseseseesseeees 66 
O, WULEINIS RUAN Le occu sccasee eee soe le Re ee ee 66 
cM ATINUS MOV BS) 6) 0 Gee perce eee eee ee erate ce ee ee 68 
Pe AID SIZING AINIDT LOCATION correc citi cescsnssssnsnessnacsonsoearsoses Ue 
MPU SaaTeAN COMME ATs © Aig te aout ah Oe Oe Son i vedeekiucctavelesebscveseadeeese ps 

Eo Orr © ell epee ae, -: eee na cc Soe esta... RI. PS 

Hem On comma arte Timmer WO CLS LTT ooo 2. SeccscvcveceGaecscs candeateeens teeeee+s-ceectetertersccorseees. 3 

See eri an canialeeantnimiao 11 tya..:...; Ree sess 74 
he TRS SST, USAR acs cc cacao qe), 
fea ARD BOOM STRUCTURAL DESIGN wo cicccccecciccescdecstssscstecesessseasbienee. 8 5 
i OHCTNCTRS,, each; ns oh ee dee 8 5 
SINE GIO “SUPA Se (ey whe: 18 U0) 1 Oe ener ner oer ee Be 8 6 
ScMmCMESITCSS “ETI VatlOM. Semeeinn nemesis... 90 
ah. [PAULUS 6 Ce) IG Meena ree ocegee scene econ 6 ot eee 70 
pean sboonmloading model ymin... eee... Oe 
OLY UACCIr BU SCI Ire Oe) 1 yet cee een see ahr eS: 
Se oOUmuctinal testing of boom support structure..........c..c0n.... we 
LENDS SI ore 10] 
a ede NU a cy oc oa cc > aa dca d)sde ll phased cub autess MMOGs 1s sas cases vee 101 
B. EQUATIONS OF MOTION--LINEAR MOMENTUM .........ccccsseseseeeee Pou 
ied Walk ADT OTN) CCU ACTOS 2055 2200.9, sore. ss saccde eee potcccisssse++:iedsboveescesresseveesces iO: 

2. Modeling of linear momentum equation elements.................. BOS 


a. Thrust modeling.................:c:sscecccs:0mtNeeenens een eens 10s 


b. Lift modeling............ccusceneee 106 

C EQUATION OF MOTION--ANGULAR MOMENTUM .............ccesscceeeees ] 18@ 
1. Derivation of equation..............:sscseesssssssseueeeess 110 

2. Modeling of angular momentum equation elements............... lnk 
a. Pitching Moment due to wing................ see liz 

b. Pitching Moment due to canard deflection iia 113 

c. Pitching Moment due to vane deflecnhon........2-2--- eee 114 

d. Pitching moment due to piteh Tate... 116 

D. COMPUTER CODE ......ccccccscesetiaceccsccdsnssensoenesslentne: conten leew 
1. Program features..............-sssseusssaieees eee ee ee 

2. Program description.......... Geet. ..050-s.eee laley 

E RESULTS AND DISCUSSION noon. ..cc6.00.5.0-s0 usu rez 
1. Time step determination. ..csms........:..0::-..0eeee b22 

2. Takeoff, transition to forward flight)2ese ee 2, 3 
a. Thrust-to-weight envelope@w......c:23205..... a, 2 

b. Takeoff transition tradeoffs:....2....2 a hes 

3. Horizontal flight .....4...:.:...s.00eseeaes ee Pl 

a. Trim conditions ............: 23 3e.......3 lal 

b. Cruise flight short period Characterisics.......... =e a) 2 

c. Predicted short period characteristics 9.0 L32 

d. LONG PeriOd .o......c.cccc. ccc teehee. cesescousernesesst ennai es) (e: 

4. Transition from horizontal flight to vertical landing.............. 139 
a. General o.com ie 

b. Level deceleration profile......\ came... ee 140 


Vill 


ES PGOMMANO SOMOS DS sacs a es re 142 


| NCSC) EOS COURS ecco coor nee 144 
EMPIRE Se NEM l o.oo cccanscesssvrveacscesssosconscdesssasevensdaccontsesnccssrossacconnduneasdtaennsiis 144 
Se Neg erm es, fash) vase sacevarsexskesbvvar ttt Bulan aLgieceeedsoneaebesesdice neces. 145 

co LE EU Aa SITS a a 145 
meevine spar desicn, fabrication Amd teSt ...i.......cccc..cccsesceseeeceeeeeseeee 145 
aCe IMO CSO TIMER ETT TTeeTO TEER OTT ST es 500. .ckicsacssesensteseesevsnnce coettttteeeseccece 145 
“oy SANTI UCC ee eRe co o-25e eet a a 146 
et yd S CONFIGURATION ecriiteiiertits tree cnscsssteteldacessesceccssceoncosssbetaven 147 

Mt PO MIVIENIDATIONS.........c2.cccsscscsessencssssossconseconseossonssoacconesenscensssssseascsoseransenacens 148 
ee Me NRO tse Sede vous sspca dhe satst sg ROOMY, «. ..:s RMR rs axe seenveraen 148 
rag aN eee eee eo 2ee sada ceoe es sc cvonssasese-nesascnvaseecuiscdbccunsactacns 148 

|, TESTS (Oe eee See See ee 148 

2, VAAN Soaue GIGS 50) a Rese a 149 

BME ATV ARCURGEE SV OT 02 getis ssc. a:ssssccsesdttvsassecccclassocecteassocclecsssssecasscoveveactsvevesseeers So) 

SS UNI! ALI OMMMMNMR 6s eas ss seed 2easeacse-cvesesseccosecssoessovee svonesdadlbvavtesiaesuvstscseene 150 

Me CO TL LN ON eres sees e sul esscescetendissce sete escssovecosossderccerecesecsscacenes Jha 

PmecteialedanGeainectional flrohtemechanics........00..cct-..0--eseereseeeeee ip 2 

OLLI (SIS ele ae 152 

|" TESTS ATES] SS Po 154 
APPENDIX A 

0 SILIE.S3 aaa s6 coe 2 he 


APPENDIX B 


ie POGRAPHS OF FABRICATION AND TEST................0.-c0c0sesecesseseneeeees 164 


APPENDIX C 


VORTEX-LATTICE CAMBER PROGRAM AND DATA. ..... oo ccececceeseetecees 169 
APPENDIX D 

CODE FOR PROGRAMS "XPLOT” AND "INCIDENCE ae 171 
APPENDIX E 

SIMULATION ASSUMPTIONS sisisecceccccsccssiee cone 180 
APPENDIX F 

FLOWCHART FOR SEMULATION PROGRAMS «ee 185 
APPENDIX G 

CODE FOR SIMULATION PROGRAMS .............00....0 ee 2 
APPENDIX H 

SIMULATION PROFILES .00..2..0.00..:c01:t01.000000@ledes so 216 
INITIAL DISTRIBUTION LIST. occ kicst... 1st s.acnsaen ar 


I. INTRODUCTION 


A. MISSION NEED 

Among the many news stories reported during Desert Storm was 
that of a group of Iraqi soldiers waving a white flag to a Unmanned 
Air Vehicle (UAV). This particular UAV, the Pioneer, had been 
providing targeting information to a battleship firing 16 inch shells 
from over the horizon. As the Pioneer helped refine the battleship's 
aim, the Iraqis knew that the small craft which flew over their 
position and the ship raining metal upon them were a team. Lacking 
any direct contact with those operating the big guns, they 
Surrendered to the small package of hardware and _ electronics 
orbiting above them. The utility of the UAV, long known to military 
commanders, was graphically demonstrated on world-wide TV. 

UAV's, formerly known as Remotely Piloted Vehicles (RPV's) can 
serve military commanders in a variety of missions. Carrying optical 
and/or infrared (IR) sensors, they provide targeting, intelligence and 
reconnaissance information. They can also be equipped with 
communications intelligence (COMINT) and electronics intelligence 
(ELINT) packages for relay to a ground or sea-based user. Other 
Systems that can be carried aloft by UAV's include: radar, 
communications relay, electronic countermeasures (ECM), laser 
designators, nuclear/biological/chemical (NBC) contamination 
sensors, and even weapons, all without sending military personnel in 


harm's way. 


B. LIMITATIONS OF CURRENT SYSTEMS 

The utility, then, for the UAV 1s great but there are limitations. 
Many systems, such as the Pioneer, Sky Owl, and the Aquila shown in 
Figure I.1, require special equipment or procedures for launch and 
recovery. This equipment includes pneumatic. rail-launchers or 
rocket assist for takeoff, and nets for capture at the completion of the 


UAV's mission. The landborne systems require a fairly large area for 
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Aquila UAV Shown With Special Launch and Recovery Eaiipiiei 
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takeoff and landing. Seaborne UAV's can impact the battlegroup’s 
Operations due to stringent wind and sea state envelopes for launch 
and recovery. Because of these disadvantages, other concepts have 
been explored as alternatives to the classic fixed-wing UAV. 

Rotary wing UAV's, such as unmanned helicopters and ducted- 
fan types, have the advantage that they can operate from small 
forward-deployed areas or ship decks without special equipment. 
However, systems with exposed rotors present a formidable safety 
hazard to nearby personnel. Additionally, vertical-only UAV's are 
limited by their relatively slow forward flight speeds and generally 
shorter endurance times once on-station. They simply cannot get to 
the target area fast enough or stay there long enough. 

C. MOTIVATION 

Another class of UAV's that has recently received renewed 
attention are those with vertical takeoff and landing (VTOL) 
capability and the ability to transition to wing-borne forward flight. 
These vehicles conceptually have the strengths of both fixed and 
rotary wing vehicles: able to takeoff and land vertically but with 
higher forward flight speeds. This class of aircraft is comprised of 
s€veral types, including: lift fan (eg: Ryan XV-5), tilt rotor (eg: Bell 
XV-15 and V-22), tilt wing (eg: LTV XC-142A), vectored thrust (eg: 
Harrier), Lift engines plus cruise engines (separate engines for each. - 
function, eg: Dasault Mirage IIIV), and Vertical Attitude Takeoff and 


Landing (VATOL), also known as "Tailsitters" [Ref. 1]. 


oS) 


The goal, then, for the next generation of UAV's will be to retain 
the capabilities of current systems but to provide this capability in a 


more flexible form. 


Il. BACKGROUND 


A. VERTICAL ATTITUDE TAKEOFF AND LANDING (VATOL) 
1. General 

As discussed in the introduction, one of the VITOL concepts 
capable of wing-borne flight is the VATOL type. It is this type which 
is the subject of this thesis. In 1949 the first “convertible aircraft 
congress’ was held and from this meeting came several designs, as 
depicted in Figures II.1 and JI.2 [Ref. 2]. Some of the advantages of 
the VATOL design [Ref. 3] are: 

1) Simplicity of design. The same engine is used for vertical 
lift and forward flight, normally with only simple modifications to 
allow hovering control. This simple design generally results in less 
weight penalty because there are no pivoting wings, nozzles or 
engine nacelles with the accompanying structural weight. 

2) High thrust-to-weight ratio for mission use. VATOL 
aircraft, with lower structural weight, have more excess thrust which 
can be utilized during takeoff and landing as well as in the mission 
ane a. 

mecduccas or climimated —deoradation of vertical lift 
capability due to "“suckdown" or exhaust ingestion. 

4) In general, better control at high angles of attack. 

In a 1978 study of several generic designs, the David Taylor 
Research Center [Ref. 4] concluded that VATOL was the best VTOL 


configuration, with a 12% higher payload/operating weight ratio than 





Figure 11.2: Early VATOK concept 
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the next best configuration, as shown in Figure II.3. 

The disadvantages of the VATOL concept in general include: 

1) In the case where special launch-and/or-recovery 
platforms or equipment are required, there is a reduction in the 
operational flexibility compared with conventional fixed wing 
aircraft. 

2) Most VATOL designs have no STOL capability, limiting 
their ability to handle increased takeoff weight. 

2. VATOL historical designs 

During the 1950's and early 1960's and then again during 
the Carter administration, there was much research on VATOL 
designs, including several successful manned and unmanned 
programs. The prop-driven Convair XFY-1, Figure II.4, was a 
successful design flown in 1953. Over 400 flights were made, all 
with transitions to and from horizontal flight. Later came the jet- 
powered Ryan XV-13, Figure II.5, which completed 136 flights, of 
which 104 were VATOL [Ref. 1]. Both programs were completed 
without serious incident. It is important to note the configurations of 
these experimental aircraft: all, except for the least successful 
Lockheed XFV-1 about which less has been written, had highly- 
Swept delta wings and most had reaction control augmenting 
aerodynamic control. A review of the many documents concerning 
these and other VATOL studies revealed several potential problem 
areas for VATOL designs. These include: 


1) Longitudinal trim. The jet designs all suffered from a 
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significant pitch-trim problem, encountered during constant-altitude 
decelerating transitions. According to Girard and Everett [Ref. 5], the 
XV-13 “could not be flown steady state at attitude angles between 32 
and 70 deg due to an incompatibility between thrust required for 
constant altitude, and that required for longitudinal control, the 
latter requirement being the greater... The XV-13 utilized a gimbaled 
nozzle, much like that on a rocket, to vector thrust. The solution to 
the problem was to simply pitch through this region by ‘following a 
memorized schedule of engine rpm versus airplane angle’ (to avoid 
excessive altitude gain) [Ref. 5]. 

2) Longitudinal flying qualities. Because the VATOL design 
is often operated at relatively high angles of attack, it is desirable 
that the vehicle have a high coefficient of lift and a post-stall region 
that is as free of abrupt changes in the aircraft's stability derivatives 
as possible. Longitudinally, this means desiring benign stall 
characteristics. The typical approach has been to utilize fairly low- 
aspect-ratio delta wings, trading off a lower coefficient of lift for 
good high-alpha flying qualities. Not all programs experienced 
longitudinal flying qualities problems, however. In 1961 NASA 
performed tests on a 1/5-scale RPV very similar in configuration to 
the XV-13. A three-view of the model is shown in Figure II.6. The 
report on this vehicle stated: "Transitions from hovering to normal 
forward flight and back to hovering flight could be made smoothly 
and easily in the full-scale tunnel, and the model seemed to have 


stability of angle of attack over most of the speed range. These 


flights...represented slow, constant-altitude transitions and covered a 
range of angles of attack from about 20 deg to 90 deg.” [Ref 6]. 

3) Lateral/directional flying qualities. Another problem 
encountered in the high angle of attack regime is marked changes in 
aircraft lateral-directional flying qualities. In many airplanes at AOA's 


just a few degrees above stall, the sign of directional stability can 
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Figure [1.6 
NASA 1/5 Scale Model of a Jet Powered VATOL Research Airplane 


change and aerodynamic roll control power can be drastically 
reduced. The XV-13 experienced this problem. Roll control required 
to “handle extreme sideslip angles exceeds the control available in 
level flight by a considerable margin at and just beyond the stall. 
Adequate control is available at full power [the XV-13 utilized 
reaction roll control to augment aerodynamic control] but then the 
airplane would begin to climb.” [Ref 5]. 

The preceeding review of flight test reports “en bem 
manned and unmanned VATOL research projects suggests an 
additional design objective: to design the vehicle to be flown with an 
angle of attack (a) in regions of predictable airplane response. This 
means either keeping the a _ low--not more than a few degrees 
beyond stall--or exploring higher a behavior in wind tunnel or flight 


Les UL. 


B. NAVAL POSTGRADUATE SCHOOL UAV PROGRAMS 
I. General 

UAV research at the Naval Postgraduate School (NPS) is 
conducted by the UAV Flight Research Laboratory (FRL). Members 
of this lab have experience in several UAV projects and have 
obtained assets for use in future’ proof-of-concept designs. 
Additionally, the FRL has facilities to manufacture and test UAVs. 
Resident skills include metalworking, composite manufacturing and 


model building. 


2. Aquila UAV 

Developed by the U.S. Army in the 1970's as a surveillance 
platform, the Lockheed MQM-105 is a flying wing design with a 
shrouded pusher propeller, as shown in Figure II.7. The Aquila was 
made from Kevlar composite material which resulted in a strong, 
lightweight airframe. The Aquila was launched by pneumatic rail 
and captured in a net, as shown in Figure L[.1. Lateral and 
longitudinal flight control was achieved through wing-mounted 
elevons. The Aquila had no sustained vertical flight capability. 
Several Aquila airframes were obtained by the FRL for use in UAV 
fesearch. 

3. Airborne Remotely Operated Device (AROD) 

Designed and fabricated by Sandia Research Laboratory 
under a contract for the Naval Ocean Systems Center (NOSC), the 
AROD, shown in Figure II.8, was a ducted-fan UAV built of graphite 
and Kevlar composites. It was powered by a two-cylinder, two-cycle 
engine developing approximately 26 BHP. The engine is more fully 
described in Section IV. Capable of hovering and limited forward 
flight, the AROD was controlled by four vanes mounted in the wake 
of the engine driving a _ three-bladed composite propeller. 
Stabilization was performed by an onboard computer fed by rate 
sensors and gyros. AROD characteristics are shown in Table II.J. All. - 
remaining AROD assets were obtained in 1993 by the FRL for UAV 


flight research. 
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4. Archytas UAV Notional Concept 

The work of theses by Blanchette, Ellwood, Brynstead, Merz 
and Kress [Refs. 7-11] have led to the development of a VATOL UAV 
concept combining the assets of Aquila and AROD with locally 
manufactured components. Figure II.9 shows the notional concept 
for FRL's VATOL UAV, named Archytas. It is designed to be flown in 
a hover in the same manner as the AROD, by using the vanes 
immersed in the wake of the propeller, and then transitioned into 
forward flight and flown as a conventional fixed-wing airplane, much 
like the Aquila. Following return to the landing site it will be 
designed to transition back to a vertical attitude and land from a 
hover. 

In his thesis, Lieutenant Kress designed the spar to connect 
the AROD to the Aquila wing, designed and built an engine test stand, 


and performed initial engine tests [Ref. 11]. 
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Figure IJ.7: Lockheed/US Army MQM-105 Aquila UAV 
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Figure 11.8: Sandia Lab/NOSC AROD Ducted-Fan UAV 
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Figure I1.9: NPS VATOL UAV Notional Concept, the Archytas 
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II. ARCHYTAS DESIGN--GENERAL 


A. DESIGN OBJECTIVES 
1. General 

The objectives of the Archytas program are to design, 
produce and test a proof-of-concept VATOL UAV. This vehicle will 
be used to examine the technologies and techniques applicable to a 
military UAV. A conceptual follow-on UAV to the Archytas program 
would be capable of carrying a variety of military payloads. A key 
feature of this future military UAV must be the ability to conduct 
Operations with minimal support equipment and_ personnel. 
Additionally, these operations must be conducted without adversely 
impacting the operations of other components of the war-fighting 
team. 

2. Specific 

As in any design effort, there are contraints which the 
designer strives to meet. In the case of the Archytas, the primary 
specific design objectives were: 

¢ Minimum weight. To increase the maneuver margins 
during hovering flight, high thrust to weight (T/W) is desirable. 
Clearly, a T/W less than one will prevent flight for a VATOL airplane. 
Historical information [Refs 12, 13, 14] suggest that the minimum 
acceptable T/W is in the range of 1.05 to 1.08 with a desirable T/W 


of 1.15. The thrust of the Archytas was considered to be fixed. In 


Ne 


other words, the AROD's propulsion system--the duct, inlet, propeller 
and engine--would be incorporated into the Archytas design with 
few modifications. This left weight as the design variable--low 
weight was a significant design goal. 

¢ Low angle-of-attack flight. As a result of the review of 
the historical VATOL programs another design goal emerged: that of 
operating the Archytas, if possible, at maximum angles of attack of a 
few degrees beyond stall when at flight speeds beyond slow hover. 
This limit on angle of attack was to ensure that the problems and 
uncertainties of high angle of attack flight were avoided, most 
notably in the area of longitudinal and _ lateral/directional flying 
qualities. 

¢ Structural support for a 100 Ib vehicle at 8 g's ultimate 
normal load factor. This requirement resulted from the work of 
Kress [Ref. 11]. 

¢ Minimum modification of AROD duct, to prevent loss of 
Structural integrity of the AROD. 

¢ Favorable static longitudinal stability for forward flight. 

¢ Favorable longitudinal control power for forward flight. 

¢ Minimum special equipment to launch or recover the 
Venice: 

¢ Minimum modification of Aquila wing. 
B. DESIGN TEAM 

The subject of this thesis, design and fabrication of the Archytas 


airframe, is one part of a larger effort at NPS to produce a 


controllable VATOL UAV capable of autonomous flight, as shown in 
Figure IlII.1. Research is ongoing by members of the NPS 
Aeronautical Engineering, Electrical Engineering and Weapons 
Systems curricula to achieve this goal. The topics of the research 
work include: 

¢ Control, including advanced six degree of freedom simulation. 

¢ Navigation; sensors and integration. 

¢ Communication and data acquisition. Up-link and down-link. 

Design of the airframe, therefore, required interaction with these 
other disciplines. The areas of most interest to all parties included: 
weight and center of gravity management, configuration, vehicle 
performance, aircraft control power and aircraft flying qualities. 
C. DESIGN PROCESS--ARCHYTAS AIRFRAME 

Design is by its mature an iterative process. Design of the 
ArchytaS was no exception to this rule. Several disciplines interacted 
closely in the airframe design process. As depicted in Figure III.2, 
the disciplines included: 

¢ Aerodynamics research included wing placement’ and 
orientation. This involved a strong interaction with Flight Mechanics 
to achieve proper longitudinal characteristics. Additionally, strong 
Interaction with the Structures discipline was necessary to supply 
predicted flight loads. 

¢ Flight Mechanics theory was utilized for proper sizing and 


placement of the canard control surface. 
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¢ Simulation was used to predict results of all displines' effects 
on design. 

¢ Structural analysis was performed to achieve a lightweight 
structure capable of supporting predicted aero-loads. This portion of 
the project included test of key structural components. 

¢ Propulsion tests were conducted to characterize the installed 


propulsion system. 
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D. THESIS ORGANIZATION 


This thesis is organized by subject area as shown in the 
preceeding figure. A chapter is devoted to each of the following 
subjects: 

¢ Section IV: Engine Tests 

¢ Section V: Wing spar Structural Design, Fabrication and test 

¢ Section VI: Canard Design and test 


¢ Section VII: Simulation 


Sections VIII and IX follow with conclusions and 


recommendations. 
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IV. ENGINE EXPERIMENTS 


A. PURPOSE OF TESTS 

The purposes of performing engine tests included: 

¢ Characterization of thrust versus throttle position for use in 
simulation program's equations of motion. This led directly to 
estimation of desirable thrust to weight ratios for vertical takeoff 
and landing--which in turn led to a directed weight management 
effort 

¢ Measurement of noise and comparison of thrust with and 
without factory-installed muffler system. 

¢ Determination of roll control power for use by other design 
team members. These data were also used to model pitch control 
power in the simulation program. 

¢ Characterization of engine acceleration and deceleration. 

¢ Determination of fuel flow rates. 

¢ Survey of flow velocity at AROD duct exit for possible use in 
control power modeling. 

¢ Measurement of the effects on thrust of deflecting the vanes- 
-called "Blockage effect”. 

¢ Characterization of “blowback effect"--changes in_ vane- - 
position caused by propeller wash. 

* Gaining of experience in starting and running AROD in 


anticipation of flight test. 
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B. DESCRIPTION OF TEST EQUIPMENT 
1. Engine description and _ installation 

The engine, manufactured by Herbrandson Enterprises, was 
a two-cylinder, two-stroke, horizontally-opposed, air-cooled engine 
with a displacement of 290 cc (17.4 cu in) developing approximately 
26 BHP, driving a _ graphite-composite three-bladed, fixed-pitch 
propeller. The uninstalled engine with its ignition electronics and 
without mufflers weighed 15.5 pounds. [Ref 15}k Installation 
characteristics were previously stated in Table II.1. The engine was 
fueled with a mixture of two-cycle engine oil and unleaded supreme 
gasoline at a ratio of 80:1. 

The propeller blade angle could be varied by the loosening 
of a blade retaining assembly. Propeller blade angle was left at the 
position in which it was received--14 degrees at the tip. Four vanes 
were immersed in the wake of the propeller and were designed to be 
used for control of the AROD in hovering flight. Each of the four 
vanes had a ten-inch span and eight-inch chord. End-caps 
approximately two inches high were placed on each end of the vanes, 
presumably to reduce spanwise loss of lift. 

Self-contained electrical power was designed to be delivered 
from an engine driven generator. For engine testing purposes, an 
external power source was also available [Ref 16]. The engine and 
engine installation in the test AROD, named "Annie", were considered 


representative of the eight vehicles obtained by NPS. 


2. Support and control equipment 

Control of vane and throttle servos was affected by pulse 
width modulated signals in one of two ways: (1) via an umbilical line 
connected to an IBM PC running a "C" software code, or (2) via radio 
frequency control from an off-the-shelf Futaba® radio control (R/C) 
transmitter sending signals to a receiver mounted inside the vehicle. 
Ignition was also provided in two different manners: (1) external 
power through a 28 VDC power source, or (2) on board engine driven 
generator. Several problems with electronics and ignition were 
encountered during engine tests. These problems and the control 
Set-up are described in thesis work by Lieutenant Commander 
Moran [Ref. 16]. 

3. Test stand 

The test stand was designed to measure thrust and rolling 
moment of the AROD. ~ Test stand design and manufacture was 
accomplished in work by Lieutenant Kress [Ref.11]. Modifications to 
the original design concerned improvement in accuracy of force and 
moment by replacement of the original "fish scale" apparatus with 
sensitive mechanical force gages. Three force gages were obtained. 
Two gages, with a combined capacity of 175 lb, were used in parallel 
to measure force and one 50 Ib capacity gage was mounted to 
measure rolling moment at the end of of a nine-inch moment arm, 
for a moment capacity of 37.5 ft-lbf. Figure 1V.1 shows the engine 
test stand, force and moment measuring equipment, and the test 


article, an AROD. Force was measured by adding the pair of force 
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Figure IV.I: Engine Test Stand for AROD Thrust and 
Rolling Moment measurement 


gages mounted on an adjustable plate which was in-line with the 
thrust shaft. 
C. SCOPE AND METHOD OF TESTS 

The tests were conducted on two separate occasions for a total 
engine run time of approximately four hours. Tests were performed 
at the engine run block-house located at the NPS Laboratory. Tests 
were conducted in accordance with the test matrix shown in 
Appendix A, Table A.I. Data included temperature, pressure, thrust, 
engine speed in rpm, vane position, moment reading, and--when 
applicable--computer inputs. Vane deflection was read directly off 
protractors mounted at the base of each vane and calibrated in one- 
degree increments. Further discussion of vane numbering and 
deflection convention follows. Rpm was measured by an oscilloscope 
connected to the ignition source or by a sensor mounted on the side 
of the R/C transmitter which essentially detected light oscillations 
caused by propeller passage. Velocity of airflow in the area of the 
vanes were measured with a hand-held digital wind-speed meter. 
Double hearing and eye protection were worn by all test personnel. 
Data were recorded on pre-printed data sheets and by a video 


Camefra. 


D. RESULTS AND DISCUSSION 
1. General 
The data for engine tests completed are shown in Appendix 


A, Tables A.II and ALIII. 


2. Thrust measurements 

The first task was to determine what effect the mufflers had 
on thrust and noise output. The engines were apparently operated 
by the previous user with mufflers off and a straight-tube exhaust 
installed. In their production configuration, the Herbrandson 
mufflers would not fit inside the AROD duct. Several sets of brand- 
new mufflers were contained in the spare parts obtained by NPS. 
FRL_ personnel built an exhaust manifold adapter that allowed 
installation of the mufflers in the constrained AROD duct. Thrust 
reading at a full throttle rpm of 7595 with the mufflers on was 118 
Ibf. The engine was shut down and the mufflers replaced by the 
Original straight tube exhaust system. Upon restart and re- 
establishment of full throttle, again at 7595 rpm, the thrust had 
increased over four percent to 123 lbf. Temperature and pressure 
for both runs--separated by less than five minutes--were the same. 
Noise was substantially increased by removal of the mufflers, as 
discussed in a later paragraph, but since initial simulation results had 
Showed a need for the highest available thrust-to-weight ratio 
(T/W), the decision was made to conduct subsequent tests with the 
Straight-tube exhausts. This test was repeated twice to verify a 
consistent result. The surprisingly large difference in thrust at the 
Same temperature, rpm and blade angle may have been due to drag - 
on the mufflers, which were in the wake of the propeller, 5am. 
frontal area and blockage effect of the mufflers was approximately 


32 square inches total. 


Thrust measurements were considered accurate to within 
1/2 lbf. Engine speed versus thrust is shown in Figure [V.2. A 
second order polynomial curve fit is overlaid and shows good 
correlation. Figure IV.3 shows the thrust that resulted at different 
throttle settings as well as a polynomial curve fit. The throttle 
position on this plot was arrived at by taking the computer's throttle 
position and mapping it linearly to a zero to 100 percent scale. This 
thrust versus throttle position curve fit was used in the simulation 
program, as discussed in Section VIJ. As shown in Figure IV.3, the 
thrust output at idle was over 20 lbf. This relatively high minimum 
thrust setting may cause the vehicle to be more difficult to slow 
down in horizontal flight. Further testing may be required with the 

idle adjustment set for a lower rpm. 

3. Vane control power measurements 
a. General 

The rolling moment imparted by deflection of vanes 
waS measured on the test stand. Moment readings tended to 
“wander over a range of approximately + 3 Ibf force gage reading 
which, through a nine-inch test stand moment arm, converted to 
approximately + 2 ft-lbf accuracy for rolling moment. This range of 
readings was inconsistent and at times the readings wandered 5 ft- 
Ibf from the nominal reading. The source of this relatively large 
Variation in readings could not be determined. Engine speed 
fluctuations were ruled out in that aural and measured rpm were 


constant, aS were the thrust readings. Wind in the test area was light 
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Figure IV.2: AROD Engine Test Results--Thrust Versus Rpm 
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Figure 1V.3: AROD Engine Modelling--Thrust Versus Throttle Position 
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and variable. The control vane positions were visually checked as 
steady. One possible explanation for the variable moment readings 
would be unsteady flow within the duct--flow over the engine and 
its support structure is downstream of the propeller but just 
upstream of the vanes. 

In keeping with MILSPEC convention for roll control 
deflection [Ref 17] positive vane deflections--clockwise when viewed 
from the vane tip--were those which produced positive rolling 
moments--clockwise viewed from the rear of the AROD. Zero vane 
deflection was defined as the deflection at which a vane was aligned 
with the longitudinal axis of the AROD. The "top" of the AROD when 
mounted on its side in the test stand as shown in Figure IV.1 was 
considered to be the point where two fuel vents were next to each 


other. Vanes were then numbered according to Figure IV.4 below. 
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Figure IV.4:AROD Vane Numbering Convention for Engine Tests 
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b. Dual-vane_ roll-control power 

Dual-vane roll-control power was obtained by deflection 
of vanes two and four, while one and three were set to zero 
deflection. The rolling moment response to this vane deflection, 
shown in Figure IV.5, was fairly linear and showed no propensity to 
drop off at high deflection angles. The combined effects of engine 
and propeller torque were overcome, for the two-vane case, by 
approximately 18 deg of deflection. 

c. Four-vane_ roll-control power 

Four-vane roll-control power was_ obtained by 
deflection of all vanes in the same direction so as to produce a rolling 
moment. The response to this vane deflection with throttle set to full 
open (7520 rpm) is shown in Figure IV.6 arava as with dual vane 
power, was fairly linear and showed no propensity to drop off at 
high deflection angles. The combined effects of engine and propeller 
torque were overcome, for the four-vane case, by approximately 12 
deg of deflection. 

The change in four vane roll control power with changes 
in throttle setting was investigated for two additional rpm settings: 
7240 and 6580. These rpm levels corresponded to approximately 93 
and 75 percent of full throttle thrust . respectively. As shown in 
Figure IV.7, there was little effective change in roll control power as 
throttle was reduced, even though the thrust was decreased by 25 
percent. As discussed in Figure VI], this surprising result was 


researched and it was discovered that the AROD manufacturer, 
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Figure IV.6: AROD Engine Test--Four-Vane Roll-Control Power 


Sandia National Lab, had similar results during testing [Ref 18]. 
d. Pitch control power modelling 

Data from the engine tests were used by other team 
members to model roll control power [Ref 19]. Additionally, the data 
were used to model the pitch control power due to vane deflection 
for use in the 3-DOF simulation, discussed in Figure VII. Details of 
this transformation from roll control power to pitch control power 
follow. 

Pitch control power from the vanes in the Archytas 
airframe will be implemented by deflection of vanes two and four as 
shown in Figure IV.4. They will both be deflected so as to cause 
pitching motion. In other words, to pitch the “nose” of the Archytas 
down vane two will be rotated counter-clockwise (negative for the 
roll control convention) while vane four will be rotated clockwise 
(positive under the roll control convention). Since this will involve 
deflection of two vanes, the data from the dual-vane roll-control 
power tests were used to model pitch control power. 

Dual-vane roll control was assumed to occur from the 
lift of the two vanes at mid-span, as shown in Figure IV.8, with each 
vane producing an equal lifting force, F. 

The rolling moment generated by these pair of vanes, 


then, is: 


9: 
L= 2F{ =| (L = rolling moment in ft - lbf) 
(IV.A) 


Equating this rolling moment model to the curve fit of 


dual vane rolling moment test data from Figure IV.5: 
L=1.5F =-7.86+(0.400)(VP) 
(Recall that VP is vane deflection) (IV.B) 


Solving equation IV.B for F: 


F = —4.24 + (0.267)(VP) (IV.C) 


Pitching moment will then be a result of twice this same 
lifting force applied over a moment arm from the vane quarter chord 
to the longitudinal center of gravity, as depicted in Figure IV.9: 

Which gives a pitching moment due to vane deflection 
of: 


M ne = ~2F (xCR,) 


vane 


or, Substituting for F from equation IV.C: 


M ane = ~2{-4.24 + (0.267)(VP)}(xcg, ) (IV.D) 


The number -4.24 in ft-lbf from equation IV.C 
Tepresents the rolling moment zero-vane-deflection position, as 
shown on Figure IV.5. Because the propwash characteristics will be 
different when vanes are deflected for pitching moment, and to 


simplify user input to the simulation program--i.e. user will not be 
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View looking down AROD longitudinal axis 
with vanes 1 and 3 removed for clarity. 


eps ey 
a 9" — pele 9" : 


| Vane 2 


S NNSSANSSSSASSSNA ANS AAAANSANASA 5 WANSNASNANSNAAN NASANRASSSSANS|| 


Vane 4 









Vanes 2 and 4 shown 
deflected for roll control 


Figure IV.8:AROD Roll Control Functional Model 






LONGITUDINAL c.g. 


@ oO 
bh 
ot oo / 
ge ee ae 


ve / 


Figure IV.9: Pitch Control Functional Model 
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required to know that a certain number of degrees of vane deflection 
are required to counter the swirl of propwash in the duct--this 
quantity was eliminated from the analysis. Note that equation IV.D 
accounts for the sign change required to reflect the difference in sign 
convention between the roll axis--where positive deflections produce 
positive moments--and the pitch axis--where positive deflections 
produce negative moments [Ref 17]. Note also the change in vanriable 


from VP to 6., used to define vane deflections in a pitching sense for 


the simulation program. Making these two changes to equation IV.D 
as well as inserting the final value for xcg, of 1.5 ft, which was 
arrived after considering issues discussed in Figure WI, the final 
model for pitching moment due to vane deflection, used in Figure 


VII, was: 


M,,.. = —(0.80)(6 


ee) 


where 0,,,, = vane deflection in degrees (IV.E) 


4. Effect on thrust from blockage in AROD duct 
The changes in thrust due changes im wane detlecniog 
termed ‘blockage effect", were documented by noting changes in 
thrust due to vane deflection at three throttle settings. This analysis 
was performed primarily to predict the change in thrust with vane 


deflection for improved computer modeling of thrust, particularly in 


oe) 
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a hover, where engine thrust levels are critical. A secondary goal of 
these tests was to attempt to characterize the swirl of propwash in 
the duct. The effect on thrust due to vane deflection is shown in 
Figure IV.10. Regardless of throttle setting, the maximum thrust 
was achieved at about -5 deg of vane deflection. The differences in 
indicated thrust were likely a result of changes in the local angle of 
attack, and hence vane drag, due to swirl in the duct. Apparently the 
local velocity at the vanes was approximately in-line with the vane's 
symmetrical airfoil section when the vanes were at -5 deg. At vane 
deflections off this peak-thrust-value the thrust fell off steadily, 
decreasing approximately ten percent at full throw in each direction. 
Incorporation of these results will increase the fidelity of future 
simulations. | 
5. Impact of propwash on vane _ position 

Early in engine testing it was observed that, when using 
computer control of vane position, the vanes would be at a different 
deflection when the engine was running than they had been during 
position calibration with the engine shut down. Furthermore, as 
throttle--and presumably local velocity at the vanes--was increased 
the vane position would change with no change in the commanded 
position. This effect, which came to be known as "blowback effect", 
was documented during testing completed on April 16, 1993. Vane 
position was calibrated with the engine off through a program 
written in "C" and run on an IBM PC 386. Reference 16 describes the 


computer control set-up. Figure IV.11 shows the change in vane 


oe 


deflection with the engine running. A general trend in the errors is 
apparent in Figure IV.11, with errors approximately linear from 
minus 15 to plus ten degrees and then dropping off at higher 
commanded vane positions. 

Engine tests performed on April 27, 1993 were conducted 
using a Futaba® R/C set to command vane and throttle positions. 
When using radio control, this blowback error did not occur. 

6. Flow speed at vane _ locations 

The speed of airflow at two locations near the vanes was 
measured with a hand-held digital wind meter for future use in duct 
aerodynamic modeling and/or control power prediction. Readings 
fluctuated approximately + 2 ft/sec. These fluctuations, which may 
have been due to actual fluctuations in flow speed or to wind meter 
inaccuracies, were considered insignificant to the goal of the analysis. 
Flow speed was gathered over a range of speeds at two different 
locations as shown in Figure IV.12. Rotor theory predicts an induced 
velocity at the disk-plane of a rotor according to the following 


formula [Ref 20]: 





"  \2pA (IV.F) 
where: 
v, = Induced velocity in ft/sec 
T = Thrust in Ibf 


p= density in slugs/ft* 3 


A= Rotor area in ft*2 
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Figure IV.10:AROD Engine Tests--Vane Deflection’s Effect on Thrust 
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Survey results are shown in Figure IV.13, with a second 
order polynomial curve fit displayed for survey position #1. Also 
shown in Figure IV.13 is the rotor theory prediction for a 12-inch 
radius rotor at the test day density of 0.002382 slugs/ft‘3. The data 
followed the slope predicted by theory but with an approximate 
average increase in velocity of 12 ft/sec. This velocity above 
prediction was probably due to a combination of three factors: 

¢ The duct was not accounted for by rotor theory Dut 
would have the effect of containing and speeding the flow. 


¢ The value for thrust, T, used in equation (IV.F ) was not 


4 
f 

f 
| 





actual thrust but rather engine thrust minus drag from the engine 
and its supports as well as duct wall and vanes. 
¢ The hub of the propeller, at approximately 15 percent of 


the propeller diameter, was larger than that on most rotor systems. 






Measurement position #2: 
At vane mid-chord and mid-span, 
approximately one inch from 

vane upper surface. 












Measurement position #1: 
At vane mid-chord, approximately 
equidistant between the two vanes 
and between center of duct and 
duct shell. 
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Figure [V.12:AROD Engine Tests--Locations of Flow Velocity Survey 


7. Miscellaneous engine tests 
a. Noise 

The sound level was measured for the purpose of 
determining hearing protection required for test personnel. 
Additionally, the noise level figured into the decision on whether to 
operate the AROD with or without its mufflers, as discussed in a 
previous paragraph. Double hearing protection must be worn in 
areas with sound levels higher than 90 dB. With mufflers not 
installed, and when within about 20 ft of the AROD at full power, the 


noise level exceeded the noise meter limit of 120 dB. Installation of 
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Figure I[V.13:AROD Engine Tests--Flow Velocity Survey 


mufflers reduced the sound levels by about 6 dB. For the mufflers- 

off case, a separation distance of greater than 100 ft was required to 

bring sound levels below the double-hearing-protection limit. These 

measurements were taken with the AROD positioned at the FRL test 

cell, which is surrounded by earth berms. Further testing will be 

required to determine sound levels for flat areas, such as that 
expected for flight test. 
b. Fuel flow 

Fuel usage was recorded for three separate engine runs 

at sustained high power, at or above 80 percent of maximum rpm. 

The average fuel consumption for these high power conditions, 


approximately representative of hovering flight, was 3.6 gallons per 
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hour. Qualitatively, fuel usage at reduced power was significantly 
less. No attempt was made to lean the engine to reduce fuel flow. 
c. Acceleration/deceleration characteristics 
The acceleration and deceleration characteristics were 
evaluated by attempting to record rpm versus time during rapid 
throttle movements. The response of the engine was too fast to allow 
observation/recording of the rpm on an oscilloscope. Qualitatively, 
by listening to engine response, the spool up and down 
characteristics were first order, with a time constant less than one 
third of a second. 
d. Reliability 
During the three data periods and several other practice 
and engine familiarization sessions, no engine related failures 
occurred. However, several failures occurred with the electronics 
which caused the engine to shut down. These failures generally 
occurred as a result of the AROD's high vibration levels. The failures 
included broken solder joints and wiring, and damaged circuit board 
components. Electronics failures and related issues are discussed in 


the thesis by LCDR Moran [Ref.16]. 
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V. WING SPAR DESIGN, FABRICATION AND TEST 


A. GENERAL 

The purpose of designing a spar was to connect the existing 
Aquila wing, and the wing root support structure, which came to be 
known as the "shoulder", to the AROD duct. Previous work by Kress 
had specified a general design of two composite spars in the shape of 
rings that would surround the AROD duct [Ref. 11]. Initial work with 
cardboard mock-ups showed this design to be impractical for 
installation on the AROD duct, due to the AROD's longitudinal struts 
and landing gear, and so the spar was redesigned. The new design 
consisted of a two spar configuration with a shape similar to the 
previous design but spanning only 90 deg of the AROD duct, instead 
of the original design's 360 deg. 

Several structural tests were conducted to determine component 
strength. A final, full configuration test was then conducted to 


determine the strength of a production-like installation. 


B. WING SPAR DESIGN AND FABRICATION 
1. Design load 
The spar was designed to withstand 8 g's ultimate load 
factor and 5.33 g's limit load factor for a 100 Ibf vehicle. This 


converted to 400 Ibf of lift load per wing for ultimate strength and 
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267 lbf for design limit. For purposes of spar design and test, this lift 
load was assumed to occur as a point force at the quarter chord of 
the mean aerodynamic chord, as shown in Figure V.1. The moment 
at the wing root, then, was this lift load times the distance from the 
mean aerodynamic chord (MAC) to the wing root, 30.5 inches. This 
resulted in a predicted moment of 12,200 and 8,128 in-lbf for 
ultimate and design limits, respectively. The load for each of the two 
Spars was then assumed to be carried as a function of its longitudinal 
distance from the point lift force. Ultimate and design limits for each 
Spar are detailed in Table V.I. A sample calculation for the aft spar 
is presented below: 

i ee i 


AFT SPAR... = 
design moment (# WINGS) ya 


sen 21) -5) 
- 16 
= 5588 in - Ib, 


Wt = design weight 
INz = load factor 
y = y coordinate distance to MAC 


= longitudinal distance between spars 


~ 
Tl 


longitudinal distance from design spar 
quarter chord of MAC 
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Table V.I 
ARCHYTAS WING SPAR STRUCTURAL DESIGN LIMITS 


ULTIMATE MOMENT DESIGN MOMENT 
(8 g's) — in-lbf (5.33 g's) in-lbf 


3812 2540 
8388 5588 
TOTAL 12200 $128 





2. Forward spar design and attachment 

Figure V.1 depicts the two-spar system conceived for use on 
the Archytas flight article [Ref. 21]. The basic structure consisted of 
a combination of graphite and fiberglass in an epoxy matrix 
Surrounding a foam core. The two spars were connected to the AROD 
in different ways. The forward spar was affixed to the duct's 
exterior By three methods {refer 10 Fieure mae. 

¢ The spar was cut to fit the duct surface as closely as 
possible and was attached with a special structural adhesive to the 
outside of the duct. Adhesive thickness varied depending on the 
quality of the fit between the spar and the duct--approximately 
1/32 to 1/8 inch. The load path for this connection method--termed 
“peel strength"--was from the surface of one composite component 
(the AROD and the forward spar), through the structural adhesive 


and into the other component. 
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¢ Grip plates were manufactured from 0.053-inch thick 
6061-00 aluminum and configured with 0.25 in nut plates. These 
grips were attached to the forward spar surface with the structural 
adhesive, with a targeted adhesive thickness of 5.0 mm. As specified 
in industrial references, to ensure a good bond between the 
aluminum grip plates and the composite/epoxy spar, the aluminum 
was bead-blasted and then cleansed by a commercially available 
trichloroethane product ("Brake Clean") [Ref 22 and 23]. Research 
was conducted to determine the feasibility of first etching the 
aluminum plates with chromic acid and then phosphoric anodizing 
them, in an attempt to optimize the bonding surface [Ref 22]. 
Discussions with personnel in the NPS Physics department as well as 
commercial anodizing contractors revealed that this process was 
environmentally risky and required special equipment and chemicals 
not readily available [Ref 24 and 25]. A decision was made to utilize 
a mechanical abrasion method (bead blasting). Additionally, the grip 
plates had number 25 holes drilled in them so that the structural 
adhesive could flow through for improved bonding [Ref 26]. The 
grips were attached to the AROD duct by 0.25 in bolts passing 
through their own load-bearing plates, which were adhered to an 
AROD strut with the structural adhesive, then through the strut and 
into the grip's nut plate, as shown in Appendix B, Figure B.1 

¢ The last connection made between the forward spar and 
the duct was through a aluminum grip plate of the same material as 


above but oriented radially. In other words, the base of the grip 
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plate (with an integral 0.25 in nut plate) rested on the duct’s exterior 
surface and this plate was attached to the forward spar at a 90-deg 
angle to the previously described grips. A _ single 0.25 in bolt was 
passed through a plate mounted on the inside of the AROD duct, as 
shown in Appendix B, Figure B.2, and into the nut plate on the grip. 
The bolt was positioned below the propeller to avoid Foreign Object 
Damage (FOD) to engine and propeller components. The purpose of 
this last connection was to improve the spar's "peel strength”. 
3. Aft spar design and attachment 

The aft spar was similar in construction to the forward spar 
except that, owing to reduced area available for mounting in the 
wing root, it was 2.0 in thick, 1.0 in thicker than the forward spar. 
The aft spar had aluminum grips mounted in the same manner as the 
forward spar. These grips were bolted to the AROD's aluminum 
landing gear legs and constituted the entire mounting strength for 
the aft spar. Appendix B, Figure B.3 shows grips installed on both 
spars. 

4. Integration of NPS spars with Aquila wing spars 

The forward and aft spars manufactured at the FRL were 

mated to the Aquila wing's root by a combination of bolts, plates and 


Structural adhesive as shown in Figure V.2. The Aquila wing was 


then attached to its wing root by its original system of bolts and nut- -~ 


plates. 
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Archytas FRL/Aquila Wing Spar Integration 


C. WING SPAR EXPERIMENTS 
1. Initial component _ testing 
Forward and aft spars were first fabricated as detailed in 
previous work [Ref 11], at which point a _ structural test was 
performed on these test specimens to compare actual strength with 
those predicted in the finite element analysis of the previous ile 
work. The structural test at this point of development was on the 


individual spar component. Figure V.3 depicts the set-up for the 


tan 
tr 


structural test on an individual spar. The upper part of the figure, 
part a, shows the spar in a notional installation. The lower part of 
Figure V.3 depicts the experiment to test spar strength. The test 
spar was mounted in a supporting fixture designed to test the 
narrow area of the spar and the bonding of the structural adhesive to 
the aluminum plate, as depicted in Figure V.2. The test spar was 
mounted upside down so that the weight placed on the beam would 
simulate a lift load. Appendix B , Figure B.4 is a photograph of the 
test Set-up. 

The 215 Ilbf weight, a stainless steel bar, was set in place by 
a hand forklift at known distances marked on the steel beam. The 
feeelbn steel] tube was then placed on the steel beam and 
progressively moved farther out the beam. Once the steel tube was 
at the beam’s end it was removed, and the 215 lbf weight was then 
moved out enough to produce the same total moment that had 
existed when the bar was in its former location with the tube at 
beam’s end. The tube was then replaced on the beam, near the test 
spar, and then again progressively moved out the beam. The weight 
of the steel beam was also accounted for and was assumed to act as a 
point force at half its length. This process continued until either a 


failure occurred or a desired moment value was reached. 
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Archytas 


With the combination of tube, beam and bar weights 
simulating an aerodynamic lift load, the moment supported by the 


test spar was calculated as follows (refer to Figure V.3): 


>! Moments =() 
ZS Ra otis) Wy (Meee = 0 
or 


Mspar = 215(1,) + 42(1,.) +W(1) 


The loading was placed slowly, i.e. the weights were not 
dropped on the beam, and were left in position for at least two 
minutes without any failure of the spar. Failure of the spar was 
quantified by listening for cracking noises in the vicinity of the spar. 
Typically, when the load was applied there would be a few small 
cracking noises as, presumably, a few of the fibers and/or portions of 
the matrix broke and "settled in". These "“micro-failures" were minor 
and short lived and, if a load was relaxed and re-applied, generally 
did not reoccur. Failure was generally marked by a slow, continuous 
Series of cracking noises which did not cease, although the time 
interval between cracks was often on the order of a minute, finally 
culminating in an actual failure of the whole structure. The time 
period from application of the load to the actual failure took as long 
as ten minutes. All of this is to say that the structural testing was 
considered to be fairly conservative, in that the moment was not 
rapidly increased in an effort to increase the apparent strength of 


the component. Table V.II shows the loading sequence for the initial 


aes, 


test of the forward spar. Recall from Table V.I that the goal was to 


reach an ultimate moment of 3812 in-lbf for this component. 


Table VII 
ARCHYTAS FORWARD SPAR TEST (BEAM=16.5 Ibf, L=74 in) 


Ly (in)(!) L, (in)() Notes (Mgpa,_ in-lbf--Goal=3812) 


= Not applied or 3 small cracking noises 
(M=2508) 
6 | a2 No racing (23012) 


10-15 small cracking noises, 
















appeared to come from plate 
area (M=3264) 


Small cracking noises continuous, 





every 10-15 seconds. Failure 
after approximately 10 minutes 
(M=3516) 
Notes: (1) Refer to Figure V.3 for dimensions and experiment 


method 


The moment at failure was 92 percent of the ultimate. 
corresponding to a load level of 7.4 g's. Post-experiment analysis of 
the test spar revealed that the failure may have started at a lower 
loading level as the laminates separated, but that final failure was 


partially masked by the test fixture constraining a portion of the 


(Nn 
ON 


spar. The failure occurred at the juncture of the shear web and spar 
cap laminates [Ref 11, pp 31] adjacent to the plate that was designed 


to transmit the loads from the Aquila's wing spar. Figure V.4 depicts 


the location of this failure. 
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2. Redesign of basic spar. structure 

As a result of the failure witnessed in the first structural 
test, a modification was made to the basic spar design in an effort to 
increase its strength. Specifically, the load path between the shear 
and spar cap laminates, which heretofore had consisted solely of the 
relatively weak matrix, was strengthened, as depicted in Figure V.5. 
This modification was made to the tip area of an existing aft spar. 
These modifications changed the weight of this particular test aft 
spar from 0.77 to 1.51] Ibf, an increase of 49 percent, but it was felt 
that the weight of a spar built with these modifications incorporated 
during original construction would not gain as much weight. In any 
case, the original design was clearly unsatisfactory and in need of 


modification. 

Modifications made 
after failure. 1-2 layers 
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graphite fibers wrapped 
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Figure V.5: Modifications Made to Archytas Spar 


3. Testing of modified spar 
The modified spar was tested in the same way as the 
original test spar, except that the beam on which the weights were 
placed was replaced with a shorter version. Table V-.III describes 
the loading pattern. The goal for this test aft spar was to reach an 


ultimate moment of 8388 in-lbf. 


Table VIII 
ARCHYTAS MODIFIED AFT SPAR TEST (BEAM=12 Ibf, L=37 in) 


Notes (Mgpar_in-Ibf--Goal=8388) 
Nothing noted (M=4314) 
Sees e386] CNothing noted 


Small cracking sounds near base 
of test fixture (M=6078) 


pik 60 Nothing noted _(M=6834) 
Nothing noted (M=5604) 


24 Small cracking sounds near base 
of test fixture (M=6078) 
As above (M=7620) 


Failure of spar at attach point to 
fixture (M=7872--94% of goal) 


Notes: (1) Refer to Figure V.3 for dimensions and experiment 



















method 
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During the testing of this modified spar, cracking sounds 
were heard and appeared to be coming from the base of the test 
Specimen, where it was held by the test fixture. As shown in Figure 
V.6, post-test analysis of the spar showed that the modified areas 
were without permanent deformation and uninjured. Appendix B, 
Figure B.5 is a photograph of this test spar showing the failure. The 
failure occurred adjacent to the area that was modified. This area 
had not been modified because it had been assumed that the test 
fixture would provide the necessary strength to support the spar, 
obviating the need to modify the entire test spar. As mentioned, the 
modified area was unaffected by the 93 percent of ultimate moment, 
correlating to 7.5 g's. While the results of this second test appeared 
to shadow the results of the first test it was obvious to test personnel 
that the modified spar was significantly improved because it did not 
bend or deform as had the previous design. Because of this 
confidence in the modified design a decision was made to 
manufacture a series of new spars using the modified technique and 


to proceed to the next stage of testing. 


Modified area (shown 
larger than actual size 
for clarity) No permanent 
deformation noted. 







Failure point a 
(adjacent to Coser 







modified area) = eww 


eee eee = = 


si eS IS 


Figure V.6: Archytas Second Spar Test Failure Description 
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4. Comprehensive” structural test 

After basic design problems were identified and corrected 
during the component testing phase, a test was conducted on an 
entire duct, wing spar, and wing assembly. An AROD duct which had 
been previously partially damaged was utilized as a Structural Test 
Vehicle (STV). The damage on the STV was limited to areas not 
affecting the outcome of structural testing. As depicted in Figure V.7, 
the STV was a standard AROD duct with a prototype right-hand wing 
assembly mounted to it. Components included the NPS manufactured 
wing spars and their attaching hardware as described earlier in this 
chapter, the canard boom with support hole drilled in the forward 
spar for it and boom base mounted to the aft spar (see Section VI), 
and the Aquila wing root bolted to the NPS spars. The STV was 
assembled using the techniques and materials anticipated for use in 
the flight article. The beam used for load application was mounted 
SO as to simulate the assumed point lift force, as depicted in Figure 
V.1; in other words, it was placed five inches in front of the aft spar 
and eleven inches aft of the forward spar. Similar to the structural 
tests conducted on individual components, the STV was loaded by a 
combination of large (215 Ibf) and small (42 lbf) weights. Because 
the wing root ("shoulder") assembly used on this full test was 
intended for use on the flight article, the decision was made to limit 
the load applied to the design limit of 5.33 g's or 8128 in-Ibf. The 


loading sequence and comments thereof are contained in Table V.IV. 
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Figure V.7: Archytas full configuration test set-up 


Table V.IV 
FULL CONFIGURATION STRUCTRUAL TEST (BEAM=12 Ibf, L=46 in) 


Notes (Mgpay_ in-Ibf--Goal=8128) 

Zl Small cracking noises heard for 

first 5-10 seconds. (M=4140) 

Appeared to be coming from 
wing shoulder 


As above (M=4644) 


Small cracking noises lasting 15- 
20 seconds (M=4896) 
Noises continued and wing 
Shoulder was visibly deforming. 
Test suspended (M=5148--63% 
of goal) 
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The wing shoulder was removed and inspected. As shown 
in Figure V.8, the wing shoulder had begun ripping at a location 
adjacent to the wing spar. No failures were noted in the FRL 
manufactured parts. A decision was made to suspend further testing 
of the wing shoulder until it could be modified. To determine the 
strength of the FRL manufactured spars and their attachment with 
the AROD duct, a test fixture was constructed to take the place of the 
wing shoulder and allow loading of just these components. This 
experimental set-up is shown in Figure V.9. Table V.V describes the 
results of this test. The decision was made to limit the load to design 
limit (5.33 g) in order to preserve the capability to do testing with a 


modified shoulder. 


Table V.V 
ARCHYTAS STRUCTURAL TEST OF FRL SPAR INSTALLATION 
(BEAM=12 lbf, L=43 in) 


me} | Nothing noved (M=3096) 

















(M=8196--101% of goal) 
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Figure V.8: Failure of Aquila wing shoulder 





Figure V.9: Structural test of wing spars without wing shoulder 
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During application of the load, there was approximately five 
degrees of bending in what would be the wing. When the load was 
removed, the structure returned to its original position. The spars 
and attachments were inspected after the test and no cracking, 


permanent deformation or loss of strength were noted. 
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VI. CANARD DESIGN, FABRICATION AND TEST 


A. GENERAL 

After formulating the initial concept of combining the wings 
from an Aquila with the AROD ducted-fan UAV, the next issue was 
how longitudinal control and stability would be provided. Use of a 
conventional, tail-aft design was discarded for several reasons, 
including too high a center of gravity for takeoff and landing and the 
Structural complexity of building a combination tail and landing gear. 
The other alternative was the use of a canard. The effort of the 
present thesis, then, was to design a canard and design, fabricate and 
test a Support structure for the canard, that would have the following 
EMiractenistics: 

¢ Good static longitudinal stability with a 15 percent static 
margin, and the aircraft trimmable throughout expected flight range 

¢ Good longitudinal control power for forward flight 

¢ Location of the canard at least one chord length in front of the 
nearest AROD structure to minimize aerodynamic interference 

¢ Lightweight construction capable of withstanding expected 


flight loads 
The process that resulted in a final canard size and location was 


iterative in nature. For example, static margin was easily increased 


by moving the canard forward, but because stronger (heavier) booms 
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would be required, this added weight to a weight-critical design. The 
iterative nature of the canard and boom design process was as 
follows: 

(1) The wing was modeled using vortex lattice techniques to 
determine what moments the canard would be required to counter. 

(2) Knowing the wing characteristics, the canard could then be 
positioned and sized to provide a certain static margin. This was 
accomplished largely through the creation of a MATLAB program 
"xplot", which displayed the aircraft's neutral point and center of 
gravity for a range of canard sizes and locations [Ref. 27]. Note that 
at this point the canard and boom weight were assigned a nominal 
value, per linear foot for the boom and per square foot for the 
canard. 

(3) Once the canard location and size were selected, further 
analysis of the static longitudinal trimmability was carried out 
through creation of another MATLAB program, named "incidence", 
whose output of pitching moment versus angle of attack graphically 
depicted the incidence angles required to fly over a range of speeds, 
or in other words, the stability and trimmability of the chosen 
design. 

(4) With the canard configuration chosen, the boom was 
Optimized by selecting the lightest tubing size capable of 
withstanding the expected loads. 

(5) Components were then manufactured and tested to assure 


that they could indeed withstand the expected loads. 
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B. WING MODELING 

The airfoil of the Aquila wing was a modified NACA 23015 
section. Modification to the airfoil consisted of a slight reflex at the 
trailing edge, affecting the last ten percent of the upper surface and 
30 percent of the lower surface [Ref. 28]. Table VI.1 shows the 
Aquila airfoil coordinates. Figure VI.1 plots these coordinates as well 
as the mean camber line, and shows polynomial curve fits for two 
segments of the mean camber line. 

The Aquila wing also incorporated a linear geometric twist of 
three deg, with the tip at a lower angle of incidence than the root 
[Ref 28]. A planform view of the basic wing geometry used for 
analysis 1s shown in Figure VI.2. The Aquila wing sweep angle of 28 
deg was retained for the Archytas design for several reasons. The 
primary reason was the relative structural simplicity and associated 
weight savings of mating the Aquila wing spar in-line with the spar 
manufactured by FRL personnel. Additionally, lower sweep angles 
had the adverse effect of reducing the static margin for a given 


center of gravity, causing longitudinal stability and control problems. 
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fables VI] 
AQUILA WING AIRFOIL COORDINATES (PERCENT OF CHORD) 


Station Upper surface Lower surface 
(x/c--percent) (z/c--percent) (z/c--percent) 
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Figure VI.1 


NACA 23015 


With Mean Camber Line Curve Fits 


Aquila airfoil, modified 
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An analysis was performed on the basic wing planform using a 
NASA-developed vortex lattice FORTRAN program, named “SUB’, 
resident on the NPS VAX computer system [Ref 29]. The wing was 
divided into 100 panels in a ten-by-ten equally spaced grid. A 
FORTRAN program, named "“camber.for" and shown in Appendix C, 
converted the derivative of the equation for the mean aerodynamic 
chord, given in Figure VI.1, to a local incidence angle for each of the 
100 panels. The output of the program "“camber.for" did not include 
the three deg geometric twist. This twist was calculated manually 
for each panel. The final input file to program "SUB" is also given in 
Appendix C. The results of the vortex-lattice analysis are listed in 


Maple: Valet: 


Table VIL 
RESULTS OF VORTEX LATTICE ANALYSIS OF AQUILA WING 


CHARACTERISTIC VALUE 


Lift curve slope, Sen ay aa} | 


per radian 
0.06581 per deg 
Slope of pitching moment/Cy 
Glin e. 


dCm/dCy = 8 pe 


Intercept of Cm/C, curve for 0.00761 
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C. CANARD SIZING AND LOCATION 
1. Static Margin 
a. General 

With the wing characteristics known, a_ Suitable 
balancing force, in the form of a canard, could now be designed. A 
Static margin of 15 percent was targeted during the longitudinal 
Stability analysis, to provide a good balance between stability and 
pitch response [Ref 30]. Figure VI.3 depicts the coordinate system 
and variables used during the static margin analysis. Note that the 
variables shown on Figure VI.3 are written in the form used in the 
static longitudinal stability program “xplot”. 

b. Computer’ modelling 

The static margin, essentially the distance between the 
whole aircraft's aerodynamic neutral point and its center of gravity, 
was determined by creating a MATLAB program named "xplot". The 
code for this program is shown in Appendix D. The algorithm for this 
program follows: 

(1) The aircraft neutral point was assumed to be due to 
a combination of two items: the Aquila wing and the to-be-designed 
canard. As with the simulation effort described in Section VII, no 
aerodynamic allowance was made for the AROD duct. The following 


equation was used to define the neutral point [Ref. 31]: 


Xow, An 3- | 
BX ac aircraft C ecg amend eI o—< rae a 


“vb da. (VILA) 
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where: 
X y- wrerate = Longitudinal positon of aircraft aero center (neutral point), normalized by c 
Xp = Aerodynamic center of wing/body combination 
C = Mean aerodynamic chord of wing (32.2 in) 
Vy = Tail volume coefficient 
Kean 
eS 
where : 
X,.- = Longitudinal position of canard's aero center (user specified in "xplot” ) 
S. = Area of canard (a vector of values in "xplot’ ) 
5S = Wing reference area (29.2 ft) 
a.) = 3-dimensional lift curve slope of canard 
a,,, = Lift curve slope of wing/body combination 


OF 


— = downwash (assumed zero) 


OM 


(2) With reference to Figure VI.3, allowances were 
made for the weight and location of the following components: AROD 
duct including landing gear, Aquila wings (including spars and 
hardware discussed in Section V), fuel, electronics pod, canard with 
servo and supporting booms, and an optional ballast box at the tip of 
the canard's supporting boom. A full description of these 
components and their weights in contained at the end of program 
“xplot", Appendix D. 

2. Longitudinal  trimmability 
The purpose of program “xplot’, discussed in the preceding 


section, was to specify the canard size, given a desired static margin 
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portion of the analysis a program named “incidence” was written in 
MATLAB to depict the incidence angles required to trim the Archytas 
over a range of angles of attack. A printout of the code for program 
“incidence’ is contained in Appendix D. The germane equations were: 


Cmoaircraft =Cmacwp Vin =D) Me (VI.B) 


where: 


Co ae Aircraft pitching moment at zero deg absolute angle of attack 
Che i Wing pitching moment coefficient at zero 


deg angle of attack (from vortex - lattice analysis) 


1. = canard incidence (measure positive with leading edge down) 


and 
AG.) C | 
—_— es = C, | (h— hac.) — Vz —bacanae 
da, Le ( roby) H a 
La (VI.C) 

where: 
d(C 
eae = Slope of aircraft pitching moment versus absolute angle of attack curve 


a 


C., = Lift curve slope of wing (from vortex lattice) 
h = distance from leading edge of mean aerodynamic chord, C, 
to center of gravity, normalized by ¢ 
hac,, = distance from leading edge of ¢ to wing/body aerodynamic center, 
normalized by ¢ 


C icceag = Lift curve slope Greanane 


The program “incidence” simply calculated equations VI.B and 


VI.C over a range of canard incidence angles, i, . and displayed the 
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results on a plot of pitching moment coefficient versus angle of 
attack. Typical flight angles of attack were calculated for a cruise 
condition, called "CR" (sea level, standard day, speed of 135 ft/sec, 95 
Ibf weight), and a slow flight condition, called “PA” for “powered 
approach" (speed of 59 ft/sec). These angles of attack, 2.3 deg and 
10.6 deg respectively, were also plotted by "incidence" to give the 
user a graphical display of the incidence angles to trim the Archytas 
over a range of angles of attack, and therefore speeds. 
3. Results 

In earlier thesis work the electronics pod, containing 
components for servo control and vehicle stability augmentation, had 
been located below the engine. Figure VI.4 shows this aft pod 
location as "configuration 1". The original AROD had the electronics 
pod located as shown on Figure VI.4 and labelled as “configuration 
oe. Initial longitudinal stability estimations for the Archytas 
configuration were made with the electronics pod in _ the 
“configuration 1" position. The results showed that an Archytas with 
ten Ibf of electronics that far aft required canard support booms on 
the order of 8 ft to achieve the 15 percent static margin goal. Even 
with this long a moment arm for the canard, its chord was less than 3 
in (remember that, due to installation reasons, the canard span was 
considered constrained between the two booms, a distance of 31.5 
in). Use of the program "incidence" showed that the aircraft 
configured with this long boom and small canard required incidence 


angles beyond + 20 deg to trim the vehicle over the operational 


Wi 


airspeed range. 

It appeared at this point that an aft electronics configuration 
was not feasible. However, before abandoning the aft pod placement, 
a tradeoff study was conducted to estimate the amount of ballast 
required to bring the center of gravity for configuration 1 far enough 
forward to satisfy stability considerations. The study compared an 
aft pod (configuration 1), with ballast placed on the tip of the boom, 
versus a pod placed forward. Multiple executions of programs 
“xplot’ and “incidence” were conducted to achieve’ similar 
longitudinal stability and trimmability characteristics. Figure VI.5 
shows pitching moment coefficient plots for each of the two different 
configurations. The two plots on Figure VI.5 are similar in that the 
chord of the canard for each configuration is about the same, 12 in 
and the aircraft can be trimmed over the same range of airspeeds 
with similar incidence angles. But, the similarity ends there in that 
the aft pod configuration required 7 Ibf of ballast! This increased the 
weight of the vehicle a corresponding amount, thereby reducing the 
thrust-to-weight ratio by over seven percent, with no added 
capability, just ballast. As discussed in the simulation section of this 
thesis, that kind of reduction in thrust to weight would likely prove 
to be unacceptable. Therefore, the aft pod configuration was 
abandoned. 

The forward pod configuration was then considered in 
detail. To minimize aerodynamic interference between the pod and 


the canard, a minimum separation distance of one canard chord 
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between the two was chosen as a design constraint. 


Shows the location of the forward pod 


and 


separation constraint. Program ‘xplot” 


for a range of canard locations. 
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Static Longitudinal Stability Comparison Between Pod Configurations 
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Archytas Forward Electronics Pod Configuration--Detail 
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To review, the constraints and goals at this point were: 

¢ 15 percent static margin (analysis tool: "xplot”) 

¢ Minimum weight (analysis tool: “xplot”) 

¢ One chord length separation between canard and pod 

¢ Trimmable over flight range (analysis tool: “incidence” ) 

¢ Canard span limited to 31.5 in by installation restrictions 

The result of runs for various canard locations (variable 
“Canardlocation" in program "“xplot") on the two programs is 
Summarized in Figure WVI.7. As noted on the figure, with a canard 
located 52 in in front of the AROD duct entrance, all constraints were 
satisfied and the weight was minimized. This canard location and 
size was therefore selected. Figures VI.8 and VI.9 show the results 
of running programs "xplot" and "incidence" for the selected canard 


configuration. 


104 


Static margin: 15% 
Source: "xplot" 


103 Canard/pod 


separation 
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Lightest configuration 
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Figure VI.7: Summary of Longitudinal Stability Analysis 
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Figure V1.8: Final results from program "xplot 
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Figure VI.9: Final results from program "incidence" 
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D. CANARD BOOM STRUCTURAL DESIGN 
1. General 

Once the canard had been sized, and placed, by the 
longitudinal stability and control analysis described in previous 
paragraphs, a final design for the boom to support the canard could 
be performed. An initial effort was made to manufacture a graphite- 
epoxy tubular boom in an attempt to minimize weight. Finite 
element analysis for composite structures was unavailable and so a 
trial-and-error method ensued. A test boom was manufactured by 
wrapping graphite in an epoxy matrix around a waxed plastic 
mandrel. The resulting structure was too flexible and easily 
crushable and was considered too weak to support expected loads. 
To prevent delay of other aspects of the project, manufacturing of a 
composite boom was abandoned in favor of aluminum. A _ tubular 
shape was chosen as a good aerodynamic shape and to allow routing 
of wiring and tubing to the canard servo and planned sensors such as 
angle of attack and airspeed. 

The installation for the canard boom is shown in Figure V.1 
and consisted of a support of the boom by a hole in the forward spar 
and a base for the boom attached to the aft spar. Because the 
majority of the boom would be forward of the forward wing spar, the 


boom was modeled as a cantilever tubular beam as show. in Figure: 


Pe 10. 
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Figure VI.11: Archytas Canard Boom Elemental Portion 


With Forces and Coordinate System shown 


The following substitutions are then made: 
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which results in: 


d tidc+dk) Laa=0 
re 


OT: 


dei | “dA 0 


Converting to polar coordinates (i.e. y=rsin@ and dA=trd@) results 


pote 
dti+ {rsin@trd@ =0 


0 


Eliminating t from the above equation and removing all non-@ 


terms from the integral and then evaluating the integral results in: 


a 








"fsinddo = O 
(VIB) 


(ae 6+d0 — 
sr cua a O 
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Consider for a moment the bracketed term to be evaluated: 


[-cos6]¥ + 4° = —cos(8+d6)+cosé 
=—cos@cosd@+sinOsind8+cos@ 


(making a small angle assumption :) (VI1.C) 
=—cos@+sin@d0—cosé 
=sin@d@ 


inserting equation VI.C into VI.B and then solving for dT results in: 


oe 


integrating both sides results in: 


= Fr? sin @d@ 
Jar= _——— 


_ Fr*cos0 


2 ] 


the moment of inertia for a tube (l=nr>1) is then inserted, giving 


the final expression for shear stress: 


_ Fcos@ 
rr (VI.D) 
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os 


3. Tensile stress derivation 
Consider next the term or 


Mzy 


OR Se 7 


Converting to polar coordinates and substituting l=nr-t and 


M,=-Fx results in the final expression for tensile stress: 


Fxsin@ 


O,=- 
‘ mrt 


(VI.E) 


4. Failure’ criteria 

The Von Mises stress criterion was chosen because it is less 
conservative than the Tresca criterion [Ref 32]. A less conservative 
approach was considered appropriate because of the unmanned and 
weight critical nature of the project. A factor of safety was also built 
in because material analysis was based on the yield stress, not the 
ultimate stress. Following is the development of the Von Mises 
stress for the loading shown in Figure VI.10 and developed in 


previous paragraphs. 


akaise 


Ove. Mises =O 


et Aru 4.) +6, frre safir +o 


=20,°+6T° 


oh 8, 


Substituting equations VI.D and VI.E gives the final general 
expression for Von Mises stress: 
Pe Gea Sing Ger scose 


0 SU 
vin rrr? nr’ t? (VI.F) 


To find the position of the highest Von Mises stress, the 
derivative of equation VI.F is taken with respect to 8 and set equal 


to zero and constant terms eliminated. The result: 


to 





qo E -3]sin6coso) =(0 
(VI.G) 


Which is satisfied at 8 = 0, 90, 180 and 270 degrees--and is 
termed the first derivative criterion. To determine if each point 
satisfying the first derivative criterion was a maximum or minimum 
point, the second derivative of equation VI.F with respect to @ was 


taken. The result: 


Always > 0 


2 
vm 7-3 |c0s26 >or<or=0?? 
(VI.H) 


2 | 


The maxima will occur when equation VI.H is less than zero. 
The first factor in equation VI.H is always greater than zero and so is 
excluded from the analysis. 


The second term is less than zero only when 
(x*/r? - 3) <0 (V1.1) 


The largest allowable radius, r, due to physical limitations to 
the installation of the wing spar on the AROD, was about three inches. 
A 3.0 in radius substituted satisfies the inequality of equation VI.] 
only when x is less than 5.2 in, which is clearly an impractically short 
boom length and therefore the second term of equation VI.H is also 
excluded from the analysis. 

The range of 8 that will make equation VI.H less than zero 
are given as follows: 

cos20 <0 when: 

45 desrees < U < 135 deprees 
eit 


225 desrees < G@ < SilSedecnec. 


This range of angles is termed the second derivative 
criterion. There are only two angles which satisfy both the second - 
and first derivative criterion: 90 and 270 degrees. Substituting 90 or 
270 degrees into equation VI.F gives the final expression for the 


maximum Von Mises stress: 


_ 2F*x? 


O — 
VIN maximum Werf. 





(VI.J) 


5. Canard boom loading model 
Figure VI.12 depicts the loading on the canard’s support 
boom and a static diagram of the loads and moment arms involved. 
The canard lift was assumed to be the maximum lift that an airfoil of 


the canard’s size could generate. This value for maximum lift was 


calculated as follows: 


pen —liOlbp 
where: 


p=0.00237688 slugs/ft> (std sea level) 
V=100 knots=169 ft/sec 


Sc =canard area=(31.5)(14) in? 
C L, =canard coeff. of lift 


= iy 
Cin Osta]] =(27)(15 degrees) (VLK) 


The value for velocity came from the simulation work. The 


Canard area was the result of the longitudinal stability analysis. The 


Canard'’s coefficient of lift was considered as a worst-case scenario. 
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a. Actual installation 


L canard 
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16" 59 " 


R base 
b. Static diagram 


Figure VI.12: Archytas Canard-Boom Loading 
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6. Material selection 

Equation VI.J gives the Von Mises stress as a function of the 
applied force (one half of L, from equation VI.K--because there are 
two booms to support the canard’s lift), the distance from force 
application (a maximum value of 59 inches as shown in Figure 
VI.12), and the radius and thickness of the tube. This equation was 
implemented on a spreadsheet for 31 commercially available 
aluminum tubing configurations [Ref. 33]. The spreadsheet is shown 
in Appendix A, Table A.IV and the Von Mises strength of each of the 
configurations is shown on Figure VI.13, along with the yield stress 
mammes for the two different types of aluminum ([Ref. 32]. 
Configurations which fall below the appropriate yield line are those 
that "passed" the Von Mises criterion. The configurations that passed 
this criteria were then compared by weight. The weight per linear 
foot was calculated simply by knowing a configuration's volume 
times its density [Ref 32], as shown in Appendix A, Table A.IV. The 
configuration eventually selected was the lightest one: 6061-T6 
aluminum with a 2.5 inch outer diameter and a wall thickness of 
0.035 inches, as shown in Figure VI.14. 

Prior to purchasing the selected tubing a buckling analysis 
meeeeconaucted. With a thickness ratio, r/t, of 35, the chosen 
configuration was not considered a "thin walled" structure (the 
literature focused on tubes with thickness ratios beginning at 500) 
and therefore comparison with empirical information showed that 


buckling was extremely unlikely [Refs 34 and 35]. 
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Figure VI.14 
Canard Boom Configurations Compared By Weight 


7. Structural testing of boom support structure 
The canard booms were designed to be supported as shown 
in Figures VI.10 and VI.1. They were manufactured as specified in 
reference 21. The strength of these supports, the carry-through hole. 
in the forward spar and the boom base on the aft spar. were each 
tested. The tests were conducted in a similar manner as_ that 


described in Section V for the wing spar. The load was applied 
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gradually and incrementally using a two weight system. Figure 
VI.15 shows the test spar in its fixture. The two separate 
components tested with this one test spar are shown: the canard- 
boom carry-through hole and the canard-boom base. 

With reference to Figure VI.12, a simple statics problem 


showed that at a design load (canard lift) of 85 lb, the force at the 
boom base ("Rpa eg of Figure VI.12) would be 228 1b and the 


support’s reaction (‘"R "of Figure VI.12) would be 313 Ib. The 


support 
boom base was tested, as shown on the upper portion of Figure 
V1I.16, to 243 percent of the design load without any failure or 
permanent deformation. Appendix B, Figure B.6 is a photograph of 
the test of the boom base. The boom's forward-spar carry-through 
hole was tested, as shown on the lower portion of Figure VI.16, to 
207 percent of the design load without any failure or permanent 
deformation. Appendix B, Figure B.7 is a photograph of the test of 
the carry-through hole. The boom base was then retested by 
removing it from the vise shown in Figure VI.15 and adhering it to a 
test spar with standard epoxy. This test spar was constrained and 


the test of the boom base repeated with the same results. 


Os 


BOOM CARRY-THROUGH 
HOLE TEST POINT 





Figure W1.15: Archytas canard-boom carry-through and base test 
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Figure VI.16: Archytas canard-boom test set-up 
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VII. SIMULATION 


A. GENERAL 

The equations of motion were developed for a_ three-degree-of- 
freedom (3-DOF) simulation of Archytas motion--two displacements 
(downrange and vertical) and one rotation (pitch). The simulation 
Gide nOt address either cross-range displacement _ or 
lateral/directional rotations. The simulation program, named AWAI, 
will be discussed following derivation of the equations of motion and 
explanation of the modeling of germane forces and moments. 
Several assumptions made during formulation of the simulation are 


contained in Appendix E. 


B. EQUATIONS OF MOTION--LINEAR MOMENTUM 
1. Derivation of equations 
Figure VII.1 depicts the general coordinate system and the 
nomenclature used during formulation of the conservation of linear 
momentum equations of motion for the Archytas. 
With reference to Figure VII.1, the first two equations of 
motion are simple vector resolutions of the velocity vector into the 


inertial reference frame: 


x= Vcosy (VII.A) 


h = Vsiny (VII.B) 
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Figure VII: Coordinate System and Nomenclature for 


Archytas Conservation of Linear Momentum Equations of Motion 


The next equation of motion is obtained by applying 


Newton's second law for conservation of linear momentum along the 


Velocity vector: 


> Fatong velocity vector — (m)(accelerationalong velocity vecion) 


Summing forces and making the following substitution: 


a= V_ results in: 


Tcos@ - D - mgsiny = mV which, when solved for Y. 
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m m (VIL, 
The next equation of motion is obtained also with Newton's 


second law, this time applied along the vertical inertial axis: 
dF, = mh (VII.D) 


Consider for a moment the expression h: 


oe d(h) 
a dt 


and with h = Vsiny (equation VII.B): 


fe d(Vsiny) 
dt 
_ yy d(siny) .  d(V) 
er a + siny 
= Vcosyy + sinyV (VII.E) 


Substituting equation VII.E into equation VII.D and 


Summing forces: 


Leosy - Dsin + Tsin@ - mg = m(Vcosyy + sinyV) 
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Solving for ye 


,. Leosy _Dsiny , Tsin@ _ ; 
r mVcosy mvcosy mVcosy 


_ L__Dtany i Tsin@ g 
mV mv 





mVcosy  Vcosy V (VII.F) 


Substitution of the following: 
sin@ = sin(y + @) = sinycosa@ + sin@ cosy 


in equation VII.F gives: 





,. L_ _Dtany , Titanycosat+sing) 2 Vtany 
Y= TV ~ mV mV 


Vcosy V 





(VG 


Equations VII.A,B,C and G, then form four equations of 


motion for displacement and velocity in two axes. The variables 


solved for by finite differencing these four equations are as follows: 


X = downrange distance 

i =herehht 

V = Speed along velocity vector 
y flight path angle 
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2. Modeling of linear momentum equation’ elements 
a. Thrust modeling 
The thrust is modeled based on engine tests (see Section 
IV). Data from these tests were used to transform a user input of 
throttle position (on a scale from O% for idle to 100% for full power) 
to engine thrust. Details of the analysis to produce _ this 
transformation are contained in Section IV. When plotted, the model 


appears as shown in Figure VII.2. 
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Archytas Throttle Position vs. Thrust model 
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b. Lift modeling 

The lift of the wing was modeled as a function of angle 
of attack. The pre-stall relationship between the trim coefficient of 
lift and the trim angle of attack was obtained from Aquila wind 
tunnel tests [Ref. 28]. However, no information was available on the 
lift curve past the stall angle of attack. In the absence of high-angle- 
of-attack lift-curve data for the Aquila wing, the post-stall lift versus 
angle-of-attack relationship was assumed. The relationship between 
coefficient of lift and angle of attack is graphically depicted on Figure 
VII.3. This figure has two parts: part (a) shows the lift curve as a 
function of angle of attack in radians; because no one curve fit would 
work over the entire range of angles of attack, the lift curve was 
broken into several sections and polynomial curve fits were 
developed for each segment. Part (b) is the same plot but with curve 
fits removed and angle of attack in degrees. 


The lift was then calculated from the formula: 


f= py SE; 


where: 


p =density 

C= Sne0 

S= reference wing area 
C, =Coeincienter lint 
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SOURCE: Aquila wind tunnel data for unstalled region; assumed profile for stalled regions. 
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b. Lift model with curve fits removed 


Figure WII.3 (two graphs): Archytas Coefficient of Lift Modeling 
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c. Drag _ modeling 
The coefficient of drag, Cp , is traditionally expressed as 


a function of the coefficient of lift: 


Cp =Cp, +KC,? (VII.H) 


But for the case of a vehicle potentially operating at high angles of 
attack, such as the Archytas, this equation predicted a decrease in 
drag at angles of attack beyond stall. Unfortunately this is not the 
reality and so an alternate model for the coefficient of drag was 


assumed. The new model assumed that Ch would continue to 


increase aS a second-order function of angle of attack (a). A further 


assumption was that the relationship between Cr and'a would be the 


same for pre-stall and post-stall regions, allowing use of equation 
VII.H as follows. 

From Figure VII.3a, for angles of attack less than stall, 
the following relationship exists between the coefficient of lift and 


angle of attack: 


C, =0.0952+3.38a (Vita 


The following constants were used: 


l 
Kae = (00936 
7A Re (VIT.J) 
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where: 


AR = aspect ratio = 4.25 
e =span efficiency factor =0.8 


Cpo = 0.06 (VILK) 


Substituting equations VII.I , VII.J and VII.K into VII.H 


results in: 


Cp = 0.0608 + 0.0602a@ +1.07a7 (VILL) 


Equation VII.L is graphically depicted in Figure VII.4: 


Coefficient of drag, Cd 
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Figure VII.4: Archytas Coefficient of Drag Model 
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The drag was then modeled traditionally as: 


| 2 
D=+5pV*SCp (VII.M) 


C. EQUATION OF MOTION--ANGULAR MOMENTUM 


1. Derivation of equation 


Figure VII.5 depicts the general coordinate system and the 
nomenclature used during formulation of the conservation of angular 


momentum equation of motion for the Archytas. 
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Figure VII.5: Coordinate System and Nomenclature for Archytas 


Conservation of Angular Momentum Equations of Motion 


leh O 


With reference to Figure VII.5, the equation of motion for 
the pitch axis is developed. Applying Newton's second law for 


conservation of angular momentum: 


> Pitch Moments aboutc.g. = (Moment of Inertia) * (pitch angular acceleration) 


OT, 


M wine ate M CANARD at: M . VANE a M .PrtcH RATE — ie (VII.N) 


®: 
lI 
cy . 


solving equation VII.N for @ and substituting results 
in: 
g = Meaning + Meacanaro + Megvane * Megrmrci rare 
(VIT.O) 


With reference to Figure VII.5, the individual moments are 


expressed as follows: 


M wing = M,, — (Leosa)(xcgy) (VII.Pa) 


M cecanard = (Lift canagp COS &)(XC8¢ ) (VII.Pb) 


M vane = — (Liftyane (XC8y) (vill Pe} 


M PITCH RATE = —Moqq (VII.Pd) 


Note that, unlike Section VI's formulation for static 
longitudinal stability where dimensions forward of the leading edge 
of the mean aerodynamic chord were negative, distances in 
equations VII.Pa through VII.Pd are absolute distances and the sign 
of each term is accounted for by the moment convention shown on 
Figure VII.5. 

Substituting equations VII.Pa through VII.Pd into equation 


VII.O gives an equation of motion for pitch acceleration: 


M.— (Loos a@)(xcgy,) + (Liftc pp COS) ACE) = Ui) er ie Mg4 
i! 


YY 


qg= 


(VII.Q) 


Combining equation VII.Q with the following equations, 
defines pitch angle, rate and angular acceleration. 


q=6 (VIL.R) 


Oey (VII.S) 
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Modeling of Angular Momentum Equation Elements 
a. Pitching Moment due to wing 

The lift of the wing was modeled as _ previously 
descrmbear “(section <v ly). With reference to equation VII.Pa, the 
pitching moment of the wing about the center of gravity was 
modeled as a function of the moment about the aerodynamic center, 


taken to be the quarter chord of the mean aerodynamic chord, and 


the lift times the distance from the quarter chord to the center of 
gravity. The moment arm for the lift force, xcgw, was a function of 
where the wing was installed with respect to the center of gravity. 
The wing's longitudinal position was limited due to installation 
constraints. 
b. Pitching Moment due to canard deflection 

For the simulation, the canard airfoil was assumed to be 
the same as the main wing, to provide lift as well as longitudinal 
control. The lifting force of the canard was modeled based on the 
coefficient of lift versus local canard angle-of-attack relationship 
shown in Figure VII.6. The canard angle of attack, as shown on 
Figure VII.6, was a combination of the vehicle's angle of attack and 
the canard deflection. Changes in this angle due to pitch rate were 
considered but discounted because, at the highest pitch rate 
observed during simulation runs, the change in canard angle of 
attack due to the pitch rate was less than ten percent. At nominal 
pitch rates this change in local angle of attack was approximately one 
percent. Additionally, the assumption was made that upwash from 
the wing had no effect on the canard. This was assumed because the 
canard was over two chords in front of the wing and the ducted-fan 
was assumed to negate some of the interaction between the wing and 


the canard. 


Canard coefficient of Lift, CLe 





Local canard angle of attack,(ALPHA-Delc), deg 
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Canard coefficient of lift versus local angle-of-attack model 


c. Pitching Moment due to vane _ deflection 

The pitching moment due to vane deflection was 
modeled based on engine test data. These tests produced rolling 
moment data, which were analyzed to predict pitching moment 
effectiveness. Details of this analysis are contained in Section IV. As 
discussed in Section IV, the results of the engine tests and analysis 
showed essentially no change in control power for thrust levels from 
approximately 90 to 120 Ib (approximate thrust-to-weight ratios of 
0.9 to 1.2), which was in agreement with tests performed by the 
AROD's manufacturer [Ref 18]. Because the thrust levels Tequimes 


during hovering and initial transition were anticipated to be in this 


constant-control-power range of 90 to 120 Ib and with the 
assumption that airflow over the vanes would remain approximately 
constant following transition to horizontal flight, with increases in the 
local velocity at the vane due to forward flight speed balanced by 
throttle decrease, control power due to vane deflection was modeled 
as constant. The relationship between vane deflection and pitching 
moment, for all speeds and throttle settings, was therefore modeled 
as shown in Figure VII.7. The vanes are assumed to produce zero 
moment at zero deflection. As discussed more fully in Section IV, 
engine tests showed a swirl of the airflow in the vicinity of the 
vanes, which resulted in some moment being produced by the vanes 
at zero deflection. This swirl effect, however, was difficult to predict 
and no tests were performed to analyze the effect of two vanes, 
deflected in a pitching sense, on the airflow patterns near the vanes. 
The vane control power, therefore, was referenced to an arbitrary 


neutral position. 


Pitching moment about cg, Mcgy, ft-Ib 
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Control Vane pitching moment model 


d. Pitching moment due to pitch rate 
The moment due to pitch rate, also known as pitch 


damping, was modeled conventionally as per equation VII.Pd, using 
the dimensional derivative M_, defined as follows [Ref. 36]: 


q 
_ Ge aoe 
’ 2V (VIED) 
where: 
Cmg = pitching moment coefficient due to pitch rate 


q = dynamic pressure (do not confuse with q, pitch rate) 
5S = Wing area 
C = mean aerodynamic chord 


= spect 


With no better data available, the value for C used 


mq 
was that given for the Aquila airframe of -1.47 per radian. [Ref. 28]. 
D. COMPUTER CODE 
1. Program features 

A main program, called "AWAI1", and _ several function 
programs (similar to subroutines) were written on MATLAB® to 
perform simulation of the three-degree-of-freedom equations of 
motion. A MACINTOSH®  IIX computer was used to ease 
programming, data analysis and processing. The program included 
the following features: 

¢ Plot of vehicle's downrange and height position as well as 
pitch angle were graphically displayed in an easily interpretable 
manner. Data to enable replay of flight were stored and replay was 
available. 

¢ Initial conditions could be specified by user, or default 
values were easily input. 

¢ Program was suspended at user's discretion so_ that 
changes in control inputs could be made if desired. 

¢ Plots of all stored variables were available for later 
analysis. 

¢ Vehicle's flight path and all aerodynamic data were 
Stored for later analysis or re-use. 

2. Program description 
By way of introduction, the programs and their basic 


functions are listed below: 


© "AWAI"--performed variable’ initialization, file 
management, numerical differencing, input/output control and 
overall simulation management. 

¢ "archytas3"--stored aircraft constants, e.g., wing area, 
canard location relative to center of gravity, pitch damping 
coefficient. 

¢ "“yprimes4"--contained equations of motion VII.A, VII.B, 
VII.C, and VII.G; calculated the derivatives of downrange position, 
height, speed, and flight path angle. 

¢ "control3"--contained the model for the conservation of 
angular momentum, described in equations VII.Q and VII.R. Output 
was updated pitch angle, rate and acceleration. 

- "“THROTTLE"--converted user input from throttle 
position, in percent, to thrust in Ibf, based on engine test data. 

¢ "labeller"--provided labels for plots requested by the 
user. 

¢ ‘plotter"--provided a series of eight plots of control input 
and vehicle response. Used following program execution. 

« “Movie"--replayed the vehicle's flight path and pitch 
angle from most recent program execution (or a saved simulation 
hui). 

A flowchart of the program and several of its subroutines is 
given in Appendix F. Appendix G provides a listing of the code for all 
simulation programs. The basic program scheme, as depicted on the 


flowchart, was as follows: 


1) Aircraft constants were imported from “archytas3". 

2) The user was prompted for how many iterations were to 
be performed, what time increment to use, vehicle weight, control 
mode (i.e., use vanes or specify a fixed pitch rate) and whether or not 
the user wanted to suspend the program to make control inputs (if 
answer to last question was "yes" then user was asked “how many 
seconds between suspensions?"). The user also had the option of 
typing “1° which input default values. 

3) The user was prompted for vehicles initial height, speed, 
pitch attitude and flight path angle. 

4) Initial derivatives of the information passed in step 3) 
were calculated. Initial drag was calculated based on drag model. 

5) The output file, called "OUTFILE" which stored virtually 
all important values, was initialized. This file was updated each 
time step to allow later plotting of all variables. 


6) The following vector of primary variables was defined: 


x downrange position 
1/4 be height 
vA S speed 

Ve flight path angle 


(r} 


The equations of motion for downrange distance, height, 
speed and flight path angle were solved using predictor-corrector 
finite differencing. Euler first forward, half step and Richardson 


Extrapolation were combined as follows [Ref 37]: 


YO =Y +Ary, 


n+1 


At 
Y =¥,+> ¥, 
GiG 


Y= Viv 1+ Yat 


n+] 


— Z ] 
Ywi= IA) a 


7) The user was presented with a picture of the vehicle's 


motion, as shown in the sample output, Figure VII.8. This display 


provided the user with a display of the aircraft's position, displayed 
as a circle, and also pitch attitude, displayed as a line drawn through 


the circle. The user also received speed cues by observing the space 


between successive positions in this constant-time-increment 





process. 
100 
os 
= Prk 
a 60 Q : 
—- © 
= 6 
Ob © 
ie 40 p : 
% 
o 6 | 
a —! 
=20) Q 20 40 60 80 100 


Downrange distance, x, ft 


Figure VII.8: Sample Plot From Program AWA1--Continuously 
Updated Position and Pitch Attitude 


8 At user-specified intervals, the user was then allowed to 
view a plot of any variables desired and to make control input 
changes. Plots were presented as requested. The user had control 
over canard and vane deflections and throttle position. 

9) The program returned to step 5) and repeated this 
process until the user-specified program run time had elapsed or the 
vehicle crashed (height less than zero). 

10) Following program execution the user could type 
‘plotter’ which provided eight plots on two pages. The plots 
provided (all versus time except last plot): 

¢ angle of attack 

* speed 

* pitch attitude 

* pitch rate 

¢ canard deflection 

¢ vane deflection 

¢ throttle position 

¢ height versus downrange position 

11) The user also could type "Movie" and be presented with 
a replay of the vehicle's flight path. 

12) Ancillary programs included: 

¢ "“Mplot"--to produce a plot of all pitching moment . 
equation values versus time. Used for troubleshooting. 

¢  ‘filesaver"--to save program runs for future analysis 


*  “resurrect"--to re-load files saved by "filesaver” 


E. RESULTS AND DISCUSSION 
1. Time step determination 

The optimum time step for use during runs of program 
"AWAI" was determined by varying the time increment, called 
“DELTAT” in the program, and observing simulation behavior. To 
speed the simulation a large time step was desirable but too large a 
time increment destabilized the numerical differencing, producing 
erroneous results. Time increments of 0.2, 0.1, 0.05 and 0.025 were 
compared while performing identical maneuvers. In each case, the 
maneuver was initialized with the vehicle in a vertical attitude at a 
speed of 1 ft/sec. After allowing the vehicle to accelerate for a few 
seconds, a gentle pushover transition was initiated. Figure WVII.9 
depicts a comparison of the vehicle's angle of attack response with 
varying time increments. As shown in the figure, at time increments 
of 0.2 and 0.1 the angle-of-attack response showed a tendency to 
“spike”. This was true after any maneuver was initiated. In fact, at a 
Q.2-sec time increment the simulated vehicle would always depart 
controlled flight, as defined by angles of attack exceeding 90 deg in 
one to two seconds. As shown in Figure VII.9, the response was 
smooth at time increments less than or equal to 0.05 sec. This 
smoothness was evident during all maneuvers. However, tii 
simulation ran noticeably slower with the 0.025 increment, with no 
significant differences in response. As a result. 0.05 sec was chosen 
as the optimum time increment to balance stability of the numerical 


Scheme with program run-time. 
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Figure VII.9 
Archytas Simulation: Comparison of Time Increments 


For Identical Maneuvers 


2. Takeoff, transition to forward flight 
a. Thrust-to-weight envelope 

In order to characterize the effects of increasing weight 
on the vehicle, simulations were conducted at a range of thrust-to- 
weight ratios (T/W). In theory, any thrust-to-weight ratio greater 
than 1.0 is acceptable, given an extremely robust control system and 
the time to wait for a given flight condition to be reached. However. 
the reality pilots understand is that more thrust is better; the 


question is, “how much is enough?”. Several assumptions were made 





when considering a minimum T/W. 

The thrust was assumed to be a constant, as delineated 
in Section IV, and the weight was varied. Equivalent takeoff 
maneuvers were simulated at full thrust with weights of 115, 110, 
105 and 100 lb, corresponding to T/W values of 1.04, 1.09, 1.14, and 
1.2 respectively. The maneuver consisted of Starting in a vertical 
climb at 1 ft/sec, accelerating straight up to 5 ft/sec and then 
beginning a very gentle "“push-over", using 3 deg of vane deflection. 
The maneuver was considered complete when the vehicle had 
achieved stable forward flight with less than 10 deg angle of attack. 

The figures of merit for this T/W comparison were 
flight path and angle of attack. Flight paths were compared based on 
the altitude change from the point the "pushover" was begun. Angle 
of attack profiles for each T/W level were also compared. It was 
assumed that high-angle-of-attack (high-a) regions were to be 
avoided due to the difficult-to-control and unpredictable nature of 
flying qualities most airplanes exhibit at high a. Angles of attack 
above 20 deg were considered unacceptable for a vehicle controlled 
remotely by a person without the proprioceptive and visual cues that 
the pilot of a manned airplane has. 

Of the two figures of merit for this analysis, the angle of 
attack became the primary parameter for comparison. The flight 
paths at different T/W were all generally acceptable. In other 
words, at all T/W, owing to the gentle nature of the pushover, the 


vehicle did not lose altitude during the transition. The vehicle was 


only about 75 feet higher at the highest T/W than it was for a 
comparable maneuver at the lowest T/W. A comparison of the angle 
of attack profiles for the various T/W is shown in Figure VII.10. At 
T/W of less than 1.09 the vehicle exceeded the established angle-of- 
attack limit. 

Because the vehicle weight predicted by the analysis of 
Section VI was approximately 100 lb, a T/W value of 1.2 was utilized 


for all subsequent simulations. 
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Archytas Simulation: Thrust-To-Weight Comparison 


b. Takeoff transition tradeoffs 
Over 500 runs of program AWAI1 were performed for 


the transition from a hover to horizontal flight and several basic 


profiles emerged. These profiles and references for the Tesuiiiie 
plots are given in Table VII. For each maneuver the corresponding 
plots of downrange distance versus height ("X-H PLOT") and time 


history ("PARAMETER PLOT") are referenced. 


Tablewvllt 
ARCHYTAS SIMULATION: TAKEOFF TRANSITION PROFILES, 
DESCRIPTION AND LOCATION OF RESULTS 


MANEUVER DESCRIPTION() X-H | PARAMETER 
PROT PLOT 


l Low vane From hover, 3° Delv until | Figure] Appendix 
deflection 10° pitch attitude, then VII.11] H, Figure 
pushover pull up. ig 




































































from hover 
2 High vane From hover, 10° Delv Figure! Appendix 
ny deflection until O° pitch attitude, H, Figure 
pushover then H.2 
from hover pull-up to level flight. 
3 | Low speed, From hover accelerate Figure| Appendix 
low altitude Straight up until V=5 VII.13)| H,, Figitge 
climbing ft/sec, then 3° H.3 
gentle Delv until 45° pitch 
pushover attitude then pull-up to 
level fire: 
4 Vertical From hover, accelerate Appendix 
acceleration Straight up until V=20 H, Figuae 
to gentle ft/sec then” 3 Dely finial H4 
pushover 45° pitch gatuiude ye then 
accelerate and climb. 





Note: (1) "Delv" is vane deflection 


(2) 1S aspeed 


Figures VII.11 through VII.14 depict the results of a 
particular profile. The pull-up maneuver used consistently was a 
medium-effort pull-up defined by use of 10 deg of canard deflection 
and O deg of vane deflection. Appendix H, Figures H.1 through H.4, 
give detailed time histories of several parameters as well as control 
inputs. 

The first profile, in which the vehicle was pushed over 
using a small vane deflection, was characterized by over 100 ft of 
altitude loss and an angle of attack that peaked at 15 deg. The 
second profile, similar to the first except that the vane was deflected 
about three times as much, was marked by more altitude loss (more 
than 200 ft) and about the same angle-of-attack response. The third 
profile involved initiating a gentle pushover while climbing but at a 
relatively slow speed. The third profile's angle of attack increased 
both during the initial part of the transition and also during the pull- 
up. There was a slight loss in altitude from the apex of the 
transition. The fourth profile, similar to the third except the vehicle 
was allowed to accelerate to a higher speed prior to pushing over, 
was characterized by very low angle of attack and a continuous climb 
requiring less than 75 ft of vertical distance to complete. 

The fourth profile was the most repeatable and had the 
following advantages: 

* The angle of attack was low during the entire profile 
and was easily controlled with slight changes in vane deflection 


during the initial part of the transition and, as speed built, by canard 


deflection during the latter portion. Controlling angle of attack 
during the other profiles was not as carefree, requiring constant 
attention to ensure that a high angle of attack region was not 
enbene ae 

¢ Altitude required to complete the last profile was well 
less than 100 ft and it was all altitude gain . By contrast, the other 
profiles involved altitude /Josses up to 200 ft. In particular, 
performing high-vane-deflection pushovers similar to the first 
profile generally involved having to climb the vehicle in a hover to 
an extremely high altitude (up to 500 ft) to ensure enough distance 
below to complete the maneuver without striking the ground. The 
altitude required for these high-vane-deflection pushovers was very 
unpredictable. | 

Profiles 1, 2, and 3 had no apparent advantages and 


therefore the fourth profile was recommended. 
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Figure VII.11: Archytas Simulation: Takeoff Transition Profile #1 
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Figure VII.12:Archytas Simulation: Takeoff Transition Profile #2 
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Figure WVII.13: Archytas Simulation: Takeoff Transition Profile #3 
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Figure VII.14: Archytas Simulation: Takeoff Transition Profile #4 


130 


3. Horizontal flight 
a. Trim conditions 
The horizontal flight characteristics of the Archytas 
were simulated at two separate trim conditions listed in Table VII.II. 
These were considered the standard conditions for all horizontal 
flight analysis. Condition 1 was considered an average cruise 
condition for employment of the vehicle and Condition 2 was tested 


to compare with Aquila data at similar flight conditions. 


Table VIII 
ARCHYTAS SIMULATION HORIZONTAL FLIGHT TRIM CONDITIONS 


Pitch attitude Bae Nes 
| vane deftection Ss] So | 
Thrust SS) aL | 











a! 


b. Cruise flight short period characteristics 
The short period characteristics were explored in 
horizontal flight in Condition 1, as specified by Table VII.II, with a 
sinusoidal pitch doublet input. Figure VII.15 depicts the input, the 
response and an expanded view of the response along with the 
analysis of the dynamics. 
c. Predicted short period characteristics 
In order to compare the simulation’s behavior to that of 
the Aquila, the short period characteristics of the Archytas were 
explored in horizontal flight, at a slower speed: Condition 2 as listed 
in Table VII.II. The simulation was excited by a sinusoidal doublet 
near the natural frequency. Figure VII.16 depicts the input, the 
response and an expanded scale view of the response along with 
analysis of the dynamics. These are compared in Table VII.III to the 
Aquila’s short period characteristics, as derived from wind tunnel 


data and analytical methods [Ref 38}. 
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Figure VII.15: Predicted Short Period Response in Cruise Condition 1 
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Table VIL.III 
SHORT PERIOD CHARACTERISTICS: 
COMPARISON BETWEEN AQUILA AND ARCHYTAS 
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Natural frequency 
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Notes: (1) Wind tunnel and analytical analysis [Ref 28] 
(2) Simulation program AWAI1 
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Figure VII.16: Archytas Short Period Response at Condition 
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In view of the use of the Aquila's pitch damping 
coefficient in the Archytas simulation, the similarity of damping 
ratios between the two is a favorable indication of simulation 
fidelity. The difference between the two aircrafts period/frequency, 
a strong function of the derivative M,, can be attributed to the 15 
percent static margin of the Archytas compared to that of the Aquila, 
which has a static margin on the order of five percent [Ref.28]. 

d. Long _ period 

The long period or phugoid characteristics of the 
Archytas were simulated at both trim conditions listed in Table 
VIT.WI. The phugoid was excited by slowing the Archytas from a 
level flight trim condition by approximately 15 ft/sec and then 
returning the controls to the trim point setting [Ref. 39]. Figure 
VII.17 depicts the control input, the airspeed response and an 
expanded plot of the airspeed response with phugoid characteristics. 
Table VII.IVY shows a comparison of phugoid characteristics between 
Archytas simulation and Aquila. 

The simulation's phugoid period was within nine 
percent of the following approximation for the period [Ref. 36]: 


Ie 
reais = Von 





The simulation's long period characteristics also 


compared well to the information available for the similarly shaped 


Aquila. The only significant disparity between predictions for the 
two vehicles involved damping ratio. The difference for the two 
aircraft's damping ratios, a strong inverse function of the lift-to-drag 
ratio (L/D), may be explained by either a real or modelled difference 
in lift or drag. In as much as the lift model was based on a vortex 
lattice analysis that agreed well with data available for the Aquila, a 
likely reason for the difference in damping ratio was the drag model. 
However, both sets of long period characteristics are representative 
of aircraft flying at the simulation speed and will pose, at worst, only 
a nuisance to a control system or pilot attempting to counter the 
phugoid. 
Table VIIL.IV 
LONG PERIOD: COMPARISON BETWEEN AQUILA AND ARCHYTAS 


Gross weight (Ib) 100 


CG (percent MAC) 
Period (sec 


0.044 0.16 
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Natural frequency as] 0.40 
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Notes: (1) Wind tunnel and analytical analysis [Ref 28] 


(2) Simulation program AWAI1 


ve 


NPP YONI P GERMAN One 8 8 PPMP IPO Ata Me tate rate Mata Nata tata tatatanatat ease eattat that aha cPaaNataMantntn faaataa tate MON PPA ans MMen Naa HMMatahD adver Pe Dover Mehatate ham Mebane ata Mate PatenteAeMetaNna Aa atataeaeae 
*. 


Canard displaced until airspeed dropped 
by 15 ft/sec, then returned to neutral 


a for duration of test. 


x 
ro = 
a % 

£ 


BD 
2) 
UB 
an 
iP) 
a) 
= 
& 
5 
a 
8) 
a 
so) 
= 0 10 
VU 
2) 
0) 
SP 
P= 
=I 
> 
= 
iB) 
0) 
aw 
Si 
= 
<= 


< 





20 30 40 50 


Time,t,sec 


3 


Sotate a0 atete® 





oc 
Sc 








40 50 


>) 
ps 
oS 


20 
Time,t,sec 





Model: Archytas 
Method: Simulation 
Trim speed: 103 ft/s 






105 


See Figure VII.15 for explanation 
of terms and method of analysis 


Airspeed, V, ft/sec 





¢=0.16 4 
W 9= 0.398 rad/s | 


W ,, = 0.404 rad/s | 








Time,t,sec 


Figure VII.17: Simulation long period response at Condition 1 


138 


4. Transition from _ horizontal flight to vertical landing 
a. General 

The transition from horizontal flight to a hover, for the 
purpose of a vertical landing, was evaluated during approximately 
150 runs of the simulation program AWAI1. Two types of transitions 
were attempted, described as follows: 

¢ A level deceleration to hover parameters. From 
stabilized horizontal flight, the profile involved attempting to slow 
the vehicle without gaining altitude while progressively increasing 
the pitch attitude to 90 degrees. This profile necessitated the ability 
to generate high angles of attack during the deceleration to avoid 
gaining altitude as the pitch attitude was increased. 

¢ A “zooming transition to hover parameters. From 
stabilized horizontal flight the vehicle's pitch attitude was increased 
while maintaining unstalled angles of attack. This profile involved a 
climbing or "zooming" flight path. As the vehicle approached high 
pitch attitudes, the power was reduced from the level flight setting 
to allow the vehicle to decelerate. Pitch inputs were then made to 
Capture approximately 90 degrees for both pitch attitude and flight 
path angle. 

For standardization, the vehicle was considered to have 
entered a hover when the airspeed was less than 0.5 ft/sec, the pitch 
angle and flight path angle were within 5 degrees of 90 degrees and 


the pitch rate was less than 2 deg/sec. 
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b. Level deceleration profile 

As is common with canard designs, wherein the canard's 
aerodynamic loading is higher than the wing's, the Archytas 
simulation could not be slowed below approximately 65 ft/sec (38 
KTAS) in level flight. Attempts to further slow the airspeed resulted 
in a pitch attitude decrease with an accompanying increase in 
airspeed and loss of altitude--a classic canard stall. This occurred 
despite full-effort control inputs for both the canard and the vanes. 
Figure VII.18 shows a situational plot of an attempt to slow the 
vehicle with full effective canard and vane deflections. Appendix H, 
Figure H.5 displays the aircraft response and control inputs for this 
attempt. The canard was deflected to produce the maximum 
coefficient of lift for the canard, which occurred at approximately 15 
degrees relative angle of attack. As shown on Figure VII.5, the local 
angle of attack at the canard was a combination of the aircraft's angle 
of attack and the canard deflection. The vane was deflected to 30 
degrees, the maximum effective deflection. As shown in both Figures 
VII.18 and Appendix H, Figure H.5, the vehicle's pitch attitude 
simply could not be slowed below the 65 ft/sec range because the 
canard would stall prior to stalling of the wing. Figure VII.19 shows 
the balance of moments that prevented the pitch attitude from being 
increased. This figure shows that the canard and vane, both at 
maximum effort to pitch the nose up, were countered by the wing's 
nose down pitch moment. The high nose down pitching moment was 


due to the aircraft's high static margin coupled with the canard 
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design. Another characteristic of the level deceleration profile was 


the large downrange distance used to slow the vehicle. 
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Figure VII.18:Level Deceleration Attempt: Situational Plot 
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Figure VII.19 


Level Deceleration Attempt: Balance of Moments 


c. "Zooming" profile 
Attempts made using the second type of profile met 
with success. Figure VII.20 shows a transition first to horizontal 
flight and then back to a hover. The transition back to a hover began 
with the vehicle at approximately 70 ft/sec, at which time power 
was added and the nose raised with canard and vane as the vehicle 


climbed at angles of attack less than stall. As the pitch attitude 


neared 90 degrees power was reduced and an opposite vane/canard 
input was made to stop the pitch rate. Appendix H, Figure H.6 
depicts the aircraft response and control inputs made during this 
particular run. 

The altitude required to perform this type of maneuver 
varied from 350 to 900 feet, depending on the speed at which the 


zoom was initiated and the throttle scheduling used. 
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Figure VII.20 
Transition to Horizontal Flight 


Followed by "Zooming" Transition Back to Hover 
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VIII. CONCLUSIONS 


A. GENERAL 

The design and initial fabrication of a Vertical Attitude Takeoff 
and Landing (VATOL) Unmanned Air Vehicle (UAV) airframe was 
completed at the Naval Postgraduate School (NPS). The vehicle, 
called Archytas, was a combination of two existing UAV's--the AROD 
and Aquila--as well as locally manufactured components, including 
wing spars and support structure for a canard. The objective of 
creating Archytas was to provide a proof-of-concept platform for 
research to explore performance trade-offs and Stabiiiiy 
augmentation. A_ three-degree-of-freedom simulation was used as 
the focus of the design efforts, to validate design decisions made in 
the fields of propulsion, aerodynamics, structures and flight 
mechanics. Engine tests were conducted to determine thrust and 
control power. Structural components were designed, fabricated and 
then tested, making modifications where necessary to ensure 
sufficient airframe strength. A longitudinal control system was 
designed, validated by simulation and tested structurally. 


Specific conclusions, by section, follow. 
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B. SPECIFIC 
1. Engine tests 
An existing engine test stand was upgraded with improved 
force and moment measuring devices. The existing engine from the 
AROD was tested and shown to provide sufficient thrust and 
longitudinal control power to power the Archytas’ design. 
Specifically, thrust and control power were quantified as a function 
of rpm and throttle position. Additionally, several miscellaneous 
characteristics were explored including duct blockage, vane 
movement in the presence of propwash, the speed of the propwash, 
noise levels, fuel usage and engine dynamics. 
2. Wing spar design, fabrication and _ test 
The structural elements to connect the Aquila wings to the 
AROD duct, including wing spars and associated hardware, were 
designed, fabricated and tested on an individual basis. These 
individual tests revealed design deficiencies which were corrected 
and the components then retested with satisfactory results. A full 
configuration test was then conducted to simulate the flight loads on 
an actual vehicle. This full test demonstrated thatthe basic spar 
design was sound and identified deficiencies with the Aquila wing 
root structure. 
3. Canard design; support. structure fabrication and 
test 
A static longitudinal stability analysis was conducted to 


Choose the size and location of a canard. Detailed weight and balance 
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analysis as well as aerodynamic modeling dictated a shift in the 
location of an electronics pod from an aft center of gravity location to 
a forward one. Canard location and sizing analysis considered 
trimmability and controllability criteria over the Archytas’ expected 
airspeed range. 

Following canard design, a _ Structural analysis of a 
supporting boom structure was detailed. Empirical analysis of 
graphite booms was conducted and eventually abandoned in favor of 
an aluminum tube boom structure with graphite composite support. 
Aluminum tubing was chosen based on a combined strength and 
welght comparison. Boom support structural members, boom 
support hole in forward spar and boom base, were manufactured and 
tested with favorable results. 

4. Simulation 

Beginning with conservation of linear and angular 
momentum equations, a non-linear’ three-degree-of-freedom 
simulation program was created. The series of programs that 
constituted the simulation allowed for analysis of the Archytas’ 
longitudinal flight path and pitch axis. Simulation was used for 
validation and analysis of a wide variety of issues including: wing 
location and sizing; canard location and sizing; weight and balance 
and static margin optimization; vehicle’ thrust-to-weight ratio 
analysis; and transition profile prediction and analysis. Additionally, 
dynamic characteristics for horizontal flight were predicted and 


compared with theory and Aquila characteristics. 
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C. ARCHYTAS CONFIGURATION 

The airframe that results from the design work of this thesis is 
the Archytas. With the addition of an avionics suite for command 
and control, the Archytas will be a Vertical Attitude Takeoff and 
Landing (VATOL) Unmanned Air Vehicle (UAV) capable of hovering 


flight, horizontal flight and transition to and from horizontal flight. 
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IX. RECOMMENDATIONS 


A. GENERAL 

The Archytas program, comprised of students and _ professors 
from several departments at NPS, has so far been marked by a 
concerted effort to approach testing in a careful manner, often called 
“build-up” in the test community. This is evident in the testing of 
individual components and use of low risk platforms, including a 
remotely controlled car and conventional airplane, for validation of 
systems. The writer strongly encourages continuation of this 
laudable approach, particularly as the program enters flight test, 


with all its potential hazards. 


B. SPECIFIC 
1. Engine tests 

(a) The combination of engine tests and _ simulation 
indicated a satisfactory thrust-to-weight ratio for the Archytas. 
However, if the vehicle's weight rises beyond about 110 Ibf, further 
engine tests may be required to determine optimum thrust 
conditions. Specifically, a matrix of different blade angles, rpms, and 
carburetor mixtures could be performed to optimize thrust. 

(b) More accurate fuel flow measurements. perhaps by 
using a commercially available fuel flow sensor, could be taken 


during future engine tests. Such fuel flow data would be useful in 
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predicting vehicle range and endurance. 

(c) As discussed in thesis work by Lieutenant Commander 
Moran, several electronics failures were experienced during engine 
testing [Ref. 16]. Future electronics manufactured must be built in a 
hardy fashion. Redundancy in design is also encouraged where 
possible. 

2. Wing spar _ design 

(a) The spar design manufactured by the FRL appears to be 
Satisfactory. Modifications to the Aquila wing root will be required 
before a flight article can be assembled. These modifications might 
include: adding graphite and fiberglass material to the wing root spar 
to increase its crush strength, thereby allowing the mounting bolts to 
be tightened more than was feasible during testing; use of skin as 
structural members; and use of structural adhesive to supplement 
the present system of bolts and plates. 

(b) All design modifications should be tested to determine 
actual strength. Such tests should eventually include a full 
configuration test. 

(c) The Department of Aeronautics and Astronautics at NPS 
Should have a finite element analysis program capable of performing 
analysis on composite material intended for manufacture. The 
Department presently has a program called PAL2 but use of this- 
copy is for educational purposes and is specifically restricted from 
being used to assist in the design of any components to be 


manufactured [Ref. 40]. 
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3. Canard design 

A canard must be designed and manufactured to sizes 
specified in this thesis. In view of the fact that canards provide lift 
for the vehicle, a slightly cambered airfoil, such as a NACA 2412, is 
recommended. The recommended method of fabrication is_ hot- 
wiring of blue foam, then covering the foam with a layer of graphite 
weave and then a finish layer of fiberglass weave. Structural tests of 
a sample should be conducted. 

4. Simulation 

The fidelity of the simulation program may be improved as 
further information about the Archytas is determined. Appendix E 
discusses the assumptions made in formulation of the simulation 
program and is therefore a good source of potential improvements to 
the program. MHigh-angle-of-attack aerodynamics is perhaps the area 
that deserves the most attention. 

Most simulation runs were made using a _ thrust-to-weight 
ratio of 1.2. If the flight article's thrust-to-weight is less, the 
simulation should be repeated to determine acceptable profiles at a 
higher weight. 

Simulation results indicate several recommended Archytas 
flying techniques: 

¢ Takeoff transition should be performed by allowing the 
vehicle to accelerate nearly straight up, allowing speed to build to 
approximately 20 ft/sec and then gently pushing over with vane 


control. As the pitch attitude reaches 45 deg the vanes may be 
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relaxed to neutral, and the transition to level flight completed with 
canard control. 

¢ Short period characteristics in horizontal flight are 
predicted to have moderate damping and frequency. Care should be 
taken to avoid excitation of the short period, particularly with inputs 
near the natural frequency. 

¢ Phugoid characteristics in horizontal flight are predicted 
to be benign and should present no adverse flying qualities to the 
Archytas pilot. 

. Due to the vehicle's large static margin and its canard 
design, the transition from horizontal flight to hover will necessarily 
have to be made using a “zooming” profile. The simulation showed 
that a deceleration in level flight to a hover will not be possible. The 
preferred method will be to "zoom" the vehicle from horizontal flight, 
aS shown in section VII. 

5S. Weight control 

Thrust-to-weight ratio in any VTOL design is critical. The 
weight of candidate components, electronic or otherwise, must be 
high on the designer's list of priorities and continually stressed to 
team members. A minimum thrust-to-weight ratio of approximately 
1.09 was identified during simulation and, as a result, an absolute 
maximum weight of 110 Ib is recommended. Accurate records 
Should be kept of the weight of all components added to the 
Archytas to assist in a weight-reduction analysis, in the event the 


vehicle becomes overweight. 
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6. Lateral and directional flight mechanics 

This thesis essentially ignored the lateral and directional 
axes of the Archytas. Great care is recommended when considering 
the areas of horizontal flight directional stability and high-angle-of- 
attack lateral/directional stability and control. Wind tunnel tests 
may be useful in determining design requirements in these areas. 
Specifically, a study should be conducted to determine directional 
stability in the Archytas’ current configuration and then assess the 
need for addition of a lightweight vertical surface. Active control of 
sideslip through use of the AROD's vanes, essentially augmented 
directional stability, is possible but probably would not work well in 
the event of engine failure. 

7. Flight” test 

The author has some experience in flight test and so a few 
words on the subject are appropriate. Prior to approving any flight 
test, including tethered hover, a thorough test plan is recommended. 
This test plan should include a strong description of exactly what 
tests are to be performed and exactly what data are to be taken. 
The safety of test personnel and bystanders must be thoroughly 
considered and “what-if” questions asked and answered. A matrix of 
events, showing the order, description and purpose of each event as 
well as data required, will be a useful tool. Build-up of events should 
be considered and specifically addressed, in other words, low-risk 
events should precede high-risk ones. Checklists for operation of the 


vehicle and its systems are encouraged as a means of ensuring 





Standardization and continuity in the project. If possible, the test 
plan should be reviewed by someone with experience in UAV testing 
and/or flight test. The test plan should be approved by competent 
authority prior to the commencement of any test plan events and 
any deviations from the test plan must first receive approval from 
the same authority--the time for creative thinking is as you write 


the test plan, not as you execute it. 
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APPENDIX A: TABLES 
TABLE A.I 
AROD THRUST STAND TEST MATRIX 


EVENT | THRUST (1) DATA (2) REMARKS 
POSIT(DEG) 


STAND WITH KNOWN 
Test Team, Date, 
AROD number 
and 
Configuration 
(ignition 
SOUTCe, @cIc, 
eter. 










LOADS. 
Prestart. 






See note (A) 





0-100% by 10% | REMOVED TS, rpm, T, TSG, | Warm-up, basic cal. 
increments 


M, MSG, Pa, Ta, | First run up to full 

event number power to define 100%, 
note rpm. and TS, 
divide into increments. 


Posty-r lot. [Sand rpm 

vs. T and decide if any 
non-linearities exist. 
If linear, use similar 

TS’s for all events. 


VP Cals 










3 100% Vary all in same Same as event 2 | Attempt to find max 
direction and thrust vane position-- 
magnitude swirl effect. The 

position for max 
thrust will be called 
VPMT 

23,4 VPM. ole As above + To quantify Single vane 






















VPMT to VPMT+ max |distance from control power 
available throw in center of vane 
Positive direction, in |to x-axis of 

approximately 5 AROD. 
degree increments for 
each TS. Repeat in 

NEGATIVE direction. 











NOTE: discontinue 
when vane Stalls (i.e. 
thrust falls off) 
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5 100% ONLY 1/3: WPMT; 2/4: Double vane control 
VPMT to VPMT+ max power: vanes deflected 
available throw in so as to both cause 
Positive direction, in rolling moment in same 
approximately 5 directions 
degree increments for 
each TS. Repeat in NOTE: discontinue 
NEGATIVE direction. when vane Stalls (i.e. 

thrust falls of 


100% Repeat event 5 using Total vane control 
50% all four vanes, AT power. 
75%, in that THREE DIFFERENT 
order THRUST SETTINGS. NOTE: discontinue 
when vane stalls (i.e. 
— as —___}-thrust—falis—on | 
100%, chop to rpm _ initial, Attempt to characterize 
75% | rpm final, T thrust response during 
initial and engine deceleration. 
final, TS initial] Qualitatively 
and final, TIME | determine if the spool 
between rpm down sounds linear. 
initial and 
final. 
75% max rate 7 Attempt to characterize 
power thrust response during 
application to engine acceleration. 
100% | Qualitatively 
determine if the spool 
sounds linear. 
NOES: (A) Prior to start up check all AROD and thrust stand components for 
security. Ensure fire extinguisher is available. Brief test team on: Noise, FOD, 
and prop arc hazards. 
(1) Initially computer input value for throttle setting (TS) will be 
used to increment engine power. Following initial runs a decision will be 
made conceming linearity of TS vs. T (Thrust). If feasible, TS will be used for 
follow-on events. 
(2) Abbreviations: 
TS: Computer input Throttle Setting 
rpm: Revolutions per minute 
T: Thrust reading in LB read directly from scale. 
TSG: Thrust Scale geometry, i.e. Moment arm, slider location 
Sufficient to calculate actual thrust. 
M: Rolling moment reading in Ib read directly from Moment 
Scale: 
MSG: Moment Scale geometry, i.e. Moment am, slider location 
sufficient to calculate actual Moment. 
Pa: Ambient pressure 
Ta: Ambient temperature 
VS: Vane positions for each vane as input from computer. 
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TABLE A.2 
AROD ENGINE TESTS: 4/16/93 
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TABLE A.3 


AROD ENGINE TESTS 
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TABLE A.4 
COMMERCIALLY AVAILABLE ALUMINUM TUBE CONFIGURATIONS 


AND STRUCTURAL ANALYSIS SPREADSHEET 
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APPENDIX B: PHOTOGRAPHS OF FABRICATION AND TEST 
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Figure B.1: Archytas construction detail: through-strut bolts” 
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Figure B.2: Archytas construction detail: plate on inside of duct 
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Figure B.3: Archytas construction detail: spar grip plates 
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Figure B.4: Archytas structural test: individual spar 
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Figure B.5 


Archytas structrual test results: Aft spar failure 
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Figure B.6: Canard boom base structural test 
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Canard boom carry-through hole structural test 
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PROGRAM "CAMBER", MNODIFIED FROM PROGRAM OF SAME NANE AS USED IN AE3501 

C THIS VERSION INCLUDES THE DERIVATIVE OF TIHIE AQUILA UAV’S MEAN CAMBER 

Peta NIE NOTE. rtllIS IS@FOR THE UNITWESTED AIRFOIL. THE “ACTUAL. ATREFOIL 

™ WAS THREE DEGREES OF GEOMETRIC TWIST AS SPECIFIED IN AD-AO68 (REPORT 

C ON AQUILA). 

c 

ececvR R.B. STONEY NPS WRITTEN 2@DEC 1992 

- 
2 2 2 2 oe eo 
= ior hOGRAM ALLOWS YOU TO INPUT THE CONTROL POINTS OF THE VORTEX PANELS 
eer IND THE RESULTING AOA’S FOR A GIVEN CAMBER DISTRIBUTION. theo FLEE 
: fayeoe EASIT.Y MODIFIED FOR DIFFERENT DISTRIBUTION EQUATIONS. 


a A a a a 


FROGRAN CAMBER 

Peel 0) kel 20), 2C(120),02C(120) ,ALP(120) ,xQ0(120),XTQ(120) 
OPEN (UNIT=20,FILE='CAMBER.IN’ ,STATUS='’UNKNOWN’ ) 

OPEN (UNIT=21,FILE='’CAMBER.OUT’ ,STATUS=’ UNKNOWN’ ) 

PRitite- | INPUT THE NUMBER OF CONTROL POINTS ON PLANFORM:" 


PRINT *,’(CHORDWISE TIMES SPANWISE STATIONS)’ 
READ *,N 

PRINT *, INPUT THE Y-COORDINATE OF THE WING ROOT’ 
PRINT *,’0.0 IF PLANFORM IS WING ONLY:’ 

RERO *, YCORD 

DO si— tv. N 


Pepi oO, bee Ootl), XTO()),YCT) 
eee! ) abt s| YCORD) THEN 


NDATA = I 

END Le 

Britt <7 y (1 je, GORD NDATA,N 
END po 


Corn - 0.0 
NPFANEL - O 
NSTART - ] 
Dp 10 T-1,NDATA 
CORD = CORD + 2*ABS(XTQ(I)-XO(I)) 
Cert <) 9) scORD. y(14+1),Y¥(1) 
NPANEL = NPANEL + 1 
IF(Y(I+1).EQ. Y(1I)) GOTO 10 
XLE=XO(NSTART)-(XTO(NSTART)-XQO(NSTART))/2.0 
DO J=NSTART,NSTART+NPANEL-1 
PRINT*, J,NSTART,NPANEL, CORD 
XC(J) = ABS((XTQ(J)-XLE))/CORD 
(7 a ne eee ee ee eae e AOU T LA AIRFOIL*#*& &# ta ea tk kk ht 
Pet xec( J) bE. .«.3) THEN 
7 C( J eels 7 2-32. 154*XC( 3) 4+332.91*xC(J)**2 


+ = AdiemoexC i) ** 342190.35*xXC(J)**4 
ALE CI) =) “ATANUDZE()) ) 
ELSE 
MeGio i= 0eeo532-1.63076*x%C(J)+2.75505*xXC(J)**2 
+ -1.41616*xC(J)**3 
ALP(J)= ATAN(DZC(2)) 
END IF 
(: a a a a a ee ee ee ee a a oe 
END DO 
NSTART=,J 
NPANEL = O 
CORD = 0.0 
10 CONTINUE 


DO I = NDATA+1,N 


AGP) ) = 02:0 


END DO 
[REMEMBER THE SLOPE #& -AOA 
WRITE(Z),20) (-ALP(1),1=1,N) 
Z0 FORMAT(®@F1I0.4) 
END 
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Program "xplot W 


OQ 


Bob Stoney 3/24/93 Naval Postgraduate School 

A program to implement a survey of static longitudinal stability using 
formulas contained in Anderson's text, specifically eq'n 7.30 on pg 386 

This program analyzes a CANARD design, the Archytas UAV 

being developed at NPGS. The basic idea is to plot both the center 

of gravity and aerodynamic center versus canard CHORD (which is simply 

the canard AREA divided by it's span...in the archytas configuration 

the span is considered a constraint (to fit the canard between the "goalpost" 
canard boams). Where the lines cross is zero static longitudinal stability 
PIGerOuEDowlcre GL the Imtersection 1S varying degrees of static margin. 
(items marked with "****" are those that will be frequently changed by user) 
(items marked with @@@@@@ are those that are a vector) 

"FRL" is fuselage reference line, the leading edge of the mean aero chord 
(negative is forward), which is 9 inches aft of the top of the AROD duct. 

For the center of gravity portion, the following items' weights are accounted 


aw Pol aw 97° oP 


0 


ave 


NH oP gl? 


HO WW 


oie 


ae on al? 


% for: The AROD, the wings and hardware, the electronics, the canard boan, and 

% the canard itself. See the end of the program listing for a definition of the 
% weights. 

elear 

elle 

ola, Off 

——--—- Semele einpuesr rom user for configuration: —-=-—--— — 


GC ARAARARAANARRRARANRAARANRAAKRADRARAAAARKRARADRAAKRANRARAKRARAAARKARAARKRARRANRARARRKRERAKRKREDRRALRRAAAEAARARABRARAR 
[eo] 


disp (' Input config (l=aft pod; 2=forward pod) ') 
CONFIG=input (’ '); 


I DG 2) OO) GO DOG ALA AL AAD AACA L A ANN ARAAA NARANAARRAARARARARQNQSNAONRABRRAARNG 


— oo me veces Of Canaremererds to calculate Acg and xac:-—z---——---—- 
BARRARRARAAARAAARAAAARAAANRARAAKRARAAARAAARAARARARARAAARAAKAARRARAARAKAAAKRLKALALRARALRLABRARLKARARAKR KABA 
o 

n=100; % mumber of increments of canard chord. 

cc=linspace(.01,1.5,n); % Canard chord in feet @@@@ 


DAAARARAANRRAAAAARANRAARARARAAAANNRAARARANRAAAAKRARAAKRBRADRNRAARALRARAKRRARRARRAKRAKRKRALRL AL ALRALKAEALRAARREA AA 
be) 


#-—---—---—----—- AERODYNAMIC CENTER PORTION--—---—----—---—~---- 


ZBARARRARRARAKRARNRNAKRAAKRNRAARAARARARARARARAARAARAARAAAARRAARAARAAARRABRARARAKAKRRARRAARRALKALRABRAKRAA 
[o) 


4) 


% Uses equation 7.30 on pg 386 of Anderson's "Intro to Flight" 

disp('Input the "Canardlocation"--the distance in inches fram the duct to’) 
disp('canard c/4. REMEMBER: forward is negative....enter negative number’) 
ansl=input(' '); 
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Canardlocatiorn=ansl/12; % Distance fram duct lip to canard quarter chord 
Xacc=Canardlocation-9/12; % x position of canard's aero center from FRL ****** 
Xacwb=8/12; % x position of wing/body aero center (in feet) 
Des oy lZ: % Canard span in feet (limited by boams) a 
ac2d=5.3: % Canard's 2 dimensional lift curve slope (per radian) 
e-0.75; % efficiency factor 
ac3d=ac2d./(1.tac2d./(pi.*be./cc.*e)): sCanard's 3-D Ivft curve slope (7 7aa es 
deda=0.0; % d(epsilon-canard) /d(alpha)....the change in downwash of the 
% canard due to angle of attack change. Ref: Roskam part 6 fig 
% 8.67 page 274. 
awb=3.77075; % Wing/body lift curve slope (per radian) (fran Vortex Lattice) 
SCC" Ce, % Canard area in ft*2 @@@@@e 
pede % Main wing area in ft*2 (Used in vortex lattice program) 
Cr=3.324; % Root chord in ft (Used in vortex lattice program) 
Cen i. % Tip chord in ft (Used in vortex lattice program) 


Char=2/3* (Cr+Ct-Cr*Ct/ (Cr+ctj); % Wing mean aemedynamic Ghamd 


VH=Xacc.*Sc./ (cbar.*S): % @@RREE 
Xacabar=Xacwh/char + VH.*ac3d./awb; % eq'n 7.30 (here hn is same as Xacabar) 


DARA KAAAAKAAAA ANA AAAAAAAAAKRAARANKRANKRAAAAAAANRAAARAAAAAAARAAAAKRAAAAAARAAAAAAARAAAAAAAAAAAARDR 
a4 


$-—-— = —~—CENTER OF GRAVITY PORTION---—~--—----—----—---- 
g AAAAR AAA AAAAAAAAAAAAAAAAAAAAAAAAAAAAA AAA AA AAA RRR 
5 —==WEIGHTS: (in 1b) 
Wohallast=input ('Enter ballast weight in lb (located at boom tip)'); 
% This is a ballast to try and improve the static margin 
Warod=52; % AROD weight (updated: 2/7/93) See "equipment list" at prog. end 
Wfhuel=1.5*6.8; % Weight of fuel (eq'n shows conversion fram gals to lbf) 
% Electronics pod weight is target weight given to electronics design group 
% See equipment list at program end for desicription of actual AROD 
% electronics pod weights. 
if CONFIG=1, 
Wpod=9.0; % Aft mount configuration 
else 
Wpod=9.5; % forward mount configuration (slightly heavier due to 
% supports and aerodynamic fairings) 
end % if 
Wwings=23.2; % Weight of TWO WINGS see equipment list (updated:2/7/93 ) 
Xboammount=22/12; % Pemnament length of boam embedded in wing 
% rOOl, We tO lp or cic 
Boomlength=xXboanmountt+abs (Canardlocation); % Total canard boom in ft 
Nboam=2; % number of booms 
Boomden=0.312; % Weight of boom per linear foot from Aluminum tube design 
Wooams=Booml engt h* Boamden* Nboom; % Total weight of all boams 
Canardden=0.4; % Canard weight per square foot (Bluebird wing w/2 servos: .68) 
(Don's R/C plane w/no servos: 0.58) Took 75% for foam/glass construction. 
Woeanard=bc* cc*Canardden; % Total weight of the canard  @@@@@é@ 
Wservo=0.25; % Weight of servo (type:IC servo MS-74 7B) 
WIOTAL= (Warod+Wfuel +Wpod twWwings+WboamstWservotWbal last) *ones(cc)+Wcanard; %S@@@E 
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% === MOMENT ARMS: (ft). Reference is leading edge of MAC, fwd is negative== 
XOGAROD=2. 94/12; % For arod. 

Fuelconfig=l; a eCUcimeime laced between Stiuts **** 

% 2=AROD fuel storage location 
if Fuelconfig=1, 
XCQFUEL=-1/12; % for "external" fuel—no fuel in AROD body 
else 
XOQGFUEL=-6.75/12; % for use of AROD's fuel tanks 
end if 
ie eeneilesspodce.g. Gepends on coniiguration. For config 1 (aft pod) the c.g. 
is as presently configured on vehicle (April 1993). For config 2 it is 
assumed that the pod c.g. will be at the same location as that for the AROD. 
if CONFIG=1, 
XcgPOD=21/12; % Aft mount configuration in ft 
elseif CONFIG=2; 
XCgPOD=-26. 75/12; % forward mount configuration in ft 

end % if 
XCQWINGS=5/12; % for wings (determined by balancing wing w/camponents 2/7/93) 
XCQGBOOMS=Xboommount-Boamlength/2; % for booms 
XCQCANARD=Xacc*ones(cc)+cc/4; % for canard Gaga 
XCGSERVO=XcgCANARD; % for servo, assumes servo and canard cg's the same @@@é@ 
oa ace cc/4; % Assumes ballast is positioned on front of boom @@é@ 

MENTS (in a hH ===: 
Met bsterbal last .*xballast: % @GEeE 
Marod=Warod*XcgAROD; 
Mfuel=Wfuel*XcgFUEL; 
Mpod=wpod* XcaqPcD ; 
Mwi ngs=Wwings*XcQWINGS; 
Mbpoams=Wboams* XOgBOOMS ; 
Mcanard=-Weanard.*XcogCANARD; % Cocca cade 
Mservo=Wservo*XcgSERVO; % @GCeeee 
MPOTAL= (Marod+Mfue 1 +Mpodt}Mwings+HMiboams) *ones (cc) +Mcanardt+MservotMballast; % @@@@ 


ov? 


oP oi 


% =Total CENTER OF GRAVITY: 
XCQTOTAL=MIOTAL./WIOTAL; % @@@@GG 
XCQTOTALbar=XcgTOTAL. /cbar; % A/C's cg location normalized by wing MAC @@@@ 





ZAARKRARARARKRAKRARAADRARAKRAKRARAARKRAKRAAAALRADRAARARARARKRKRAKRARARARARARARKAAKRARAARKRAARAKRRARKRARARRARAL 


<a. .-> Be eee iene helo sO Ale) oe olAric MARGINS———-—————————— 
ZEARAR RA DRAADRKRA KARR ARARARAAKRARARAAADR AR ARAAARARBRARARADRARARARARAARRAARRARAKRARARRARRARAAKRARKRARALR 
plot flag5=0; plot flagl0=0;plotflagl5=0; % a plotting flag, see plotting section 
fem i=):n, % Note: second number must match linspace # points 
SM (i) =Xacabar (1) -XegTOTALbar (i); 2 Steele margin 
if abs (SM(i)-0.15)<=0.001, % to identify 15% static margin value 
SMLS=SM(i); % The 15% static margin values 


Xacabar]5=Xacabar (i); ae 
XcQTOTALbar15=XcgTOTALbar (i); % " 

eel 5=ce (i) 12; = 

pict flacis=1.; % to enable polyline plot later 
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acl 5=ac3d (i) ; aC) alpha at 15% sy 

Scl5=Se (i); % Canard area at 15% gy 

Cancd1 5=cc (i); % Canard Chord at 152 SM 

Weanardl 5=Weanard (j ); Sthe wr of the canard when SM=153 (put On plot) 

WIOTAL] S=WIOTAL (i) ; % the total weight when SMF15% (put on plot) 
elseif abs (SM(i)-0.1)<=0.001, id NEO) identify 10 Static margin value 

SMLO=SM(i); The 10% Static Margin values 

° sed 


SCgTOTALbar 1 0=XcgTOTALbar (4) - zo" 

CCLO=ce (i) *12- a 

Plotflagl0=1- % to enable Polyline plot later 

aclO=ac3d(i); Cl alpha at 103 SM 

Scl0=Se (i): *% Canard area at 10% gy 

Weanardl 0sWeanard (i) - sthe Wl of the canard when SM=10% (put on Plot) 

WIOTAL1O=WTOTAL (i) - % the total weight when SME10¢ (put on plot) 
elseif abs (SM (i) -0.05)<=0.001, 3 to ident ify 9% Static Margin 

SMO=SM(i); 3 The 5% Static Margin values 


*CQTOTALbar S=XcgTOTALbar (4) - Sve 
CC3=Ce (i) *12- a (Converted to inches, like Picw) 
Plot flag5=1; 3% tO enable Polyline Plot later 


aCo=ac3d (i): % Cl alpha at 53 SM 
SCcS=Se (i) % Canard area at 5% gy 
end %if 
end for 
ray cc 
Lo oe CG AND AERO CENTER ON A PLOT -— eee —---- 


Frid hae eee =" 02, Keceiaal 
Qrid; xlabel ('Canarg Chord, cc, Layee 


XS5=[ec5 CCS]; & for polyline function 
¥5=[Xacabar5 ACQTOTALbarS} ; en 
polyline (x5, v5, 1-1) % draw 5% Static Margin line 
rent (oe nae 5 XCOTOTALar5) /24 KotOR 
text (cc5, texty, "53 SM!) % put label On 5% SM line 
end %if 
098 Plotflagl0—=1, 
hold on 
X10=[cc10 CclO}; & for polyline function 
Y10=[Xacabar10 ACQTOTALbar 10} ; Se 


Polyline (X10,y10, 11 % plot 10% Static margin line 

rent (ool a O-NCOTOTAL bar 10} /2+xagtON a 

Cext (ccl0, texty, "192 obit) % put labe} on 10% SM line 
Clare 
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if plotflagl5—1, 
hold on 
wie -(Gelomeelo); % rOr polyline function 
¥15=[Xacabar15 XcgTOPALbar15]; oe ae 


poly lime (x1 Sav15,"——"} % plot 15% static margin line 

texty= (Xacabarl1 5-XcgTOTALbar15) /2+XcgTOTALbar15; 

nem (Colo, texty, 15> SM") % put label on 15% SM line 
end %if 


* Plotting "niceties": 
tex (cc (10) *12, Xacabar (10), 'A/C Neutral Point') 


text (cc (10) *12, XcgTOTALbar (10), 'A/C Center of Gravity') 


Ltop=sprintf (‘Ballast 
Ll=sprintf ('Wing area = %g',S);text(.45, .87,L1,'sc') 


L2=sprintf ('Wing/body a = %g', awb) ;text(.45, .84,L2, 'sc') 


ta Wwoallasta); tent (.45,.9, Ltop, *se*) 


meeooriner ('Can. ("fram duct) = %g’, Canardiocation*12Z) ;text(.45, .81,1L3, 'sc") 


L4=sprintf ('Xacwb = %g', Xacwb*12) - text (.45, .78,L4, 'sc') 
E>=sprintf ("Canard span = %gq',bc*l12) ;text(.45, .75,L5, 'sc") 
L6=sprintf ("XcgAROD = %q', XCgAROD*12) ; text (.45, .72,L6, 'sc') 
i-serimct(*XcgPOD = %q’, XogPOD* 12) ;text(.45, .69,L7, 'sc') 


if plotflaglS—1, 


L8=sprintf ('Xacabar (.15SM) =%g', Xacabar15) ;text (.45, .66,L8, 'sc') 


me-opranci ('CG (in) (.1SSM) = oq", XcqIOFALbar] 5*cbar*12) ; 


Uone(.45,-05,19,' Sc’) 
elseif plot flagl0O—1, 


L8=sprintf ('Xacabar(.10SM) = %q', Xacabarl10) ;text (.45, .66,1L8, 'sc') 


end 


Peee—sorintt('Config: %q',CONFIG) ;text(.7, .90,Rtop, 'sc') 


Rl=sprintf ("'Warod = %g',Warod);text(.7,.87,R1, 'sc') 
R2=sprintf (’Wfuel = %g',Wfuel) ;text(.7, .84,R2,'sc’) 
R3=sprintf ('Wpod = %g', Wood) ; text (.7,.81,R3,'sc') 


R4=sprintf ('Wwings = %g',Wwings) ; text (.7, .78,R4, 'sc') 
R5=sprintf ('Woooms = %g',Wboams) ; text (.7,.75,R5, 'sc') 


if plotflagl5—1, 


R6=sprintf ('Weanard(.15SM) = %g',Wcanardl15) ;text(.7,.72,R6,'sc') 
R7=sprintf (*Wtotal (.15SM) = %g',WIOTAL1S) ;text(.7,.69,R7,'sc') 
ee opmenirt("ce(.154M) = %qg',ccl5) > text (.7, .66,R8, *sc'} 


elseif plotflagl0=—=1, 


Ro=sprintf ('Wceanard(.1SM) = %g',Wcanard1l0) ; text (.7,.72,R6,'sc') 
R7=sprintf ('Wtotal (.1SM) = %g',WIOTAL10);text(.7,.69,R7,'sc') 


end 
if plotflagl0O—1, 


Semeebintl ( Co(.1054%) = %q",ccl0);text(.7, .63,R9, 'sc') 


end 


7s 
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SCAKRARAANRNNRANRAAAAAAAAAAAAAAAAAAAAAAABRAAAARANRAAAAAAARAARAAAAAAAAAAA 
~ 


$-—--- ——-PYrint OUL fom Prema. for 5,10,15% SM gives acta ecer se 
disp ('Configuration®) 
if CONFIG=1, 

disp('Pod Aft') 
elseif CONFIG=z2, 

disp ('Pod Fwd, gyros aft') 
elseif CONFIG=3, 

disp ('Pod Fwd,gyros fwd') 
end 
if plotflagl5—1, 
Gusp(’ ror ss ai) 


GErso«! ac3D Canard ("duct) c cherd (nie se (ae 2a) 
outl5=[acl5, Canardlocation* 12, eels, sella]; cise) 
end 


if plotflagl0O—=1, 
diset For 102-407) 


disp (' ac3D . “Canard ("aact) -e chord se fies ae 
out]0=[acl0, Canardlocation*12, cel0; Sell] ;disstout 16) 
end 


if plotflag5—1, 

Oise ( HOE saan} 

aise (* ac3sD Canard(“dict) Gc chord (in) se (sea 4a 
QUtS= (aco, Canardlocation*12, ces, seo | diseteurs) 

end 


BSAAAAAKRAAAAAARANRAAKRAAANRAAAAAAAAAARAAAABRAAAARAAAAANAAAAAAAARAAAAAANKAAAAAAAARAAAKRAAAAAAA 


= —-EQUIPMENT LIST (DEFINES WHAT IS INCLUDED IN WEIGHTS) 

AROD Weight includes (2/7/93): 

Airframe "Andy", in NPGS configuration (i.e. with no avionics between struts). 
*- Camplete engine and prop assembly (recently run) with intake 

filters installed on carbs and inlet guide vanes installed. 

* 4 vanes, 8 centerline supports (2nd set of 4 installed by LCDR Moran to 
steady aft mounted electronics pod) 

* One servo between struts (fuel management servo?). 7 1/2" aft of FRL 

* 4 vane servos. 

° 4 landing gear legs (in APOD configuration, not yet lenathened for Archytas) 
and camposite landing gear "hoop". 
, * NO: pod (or associated electronics), fuel, primer, roll cage, cr fOrmward 
OMe tr Sic Vice: 
® Weight: 50.5 1D 3) °630G (forsncon®) 
* CG position (fran FRL)> with hoco: Tie w/O hoop: Lisle 


dO _A\L 70 yO 


nN oP 


No WP 9 lO 


a?) >\0 
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POD WEIGHT includes (4/4/93): 

Pei thiatenolas Circuit Cards); skim (4" ring, 8” ring, 8" cap) 

¢ Blectronics: 3 ~7" cards (fiber obtic card removed); 1 MHz DC-DC converter; 
6 black power supplies in tray; camplete gyro/acclerameter package 

* Board that smooths generator output (the one with the big 

orange capacitor). 

° For forward mounting includes bullet supports (4), front cap (with fan and 
other electronics. 

WEIGHT: Aft mount: 10 lb Forward mount: 11.72 lb 

Woes! 6G: 1S 5 inches aft of forward lip of gyro section skin. 

See also: Figure entitled "Electronics pod dimensions--Forward pod config" 
(A MacDraw document) 


WINGS weight includes (2/7/93): 
° RH Aquila wing w/o hardware, RH shoulder; RH tip 


* Forward test spar, with 4 plates for holding end and 2 plates for 
attaching to shoulder's spar, plus 8 bolts to put it all together. 
* Aft test spar, with aft grips (grip weight estimated). 

TOTAL weight of two wings with associated camponents: 23.2 1b 


SS ———$——————————— a a aa 
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% Works only for 15% Static Margin..... 
% Must £1 ese min Sscotons 
% ##### indicates quantity to be changed by user 


DRA RARARAKRAKRARAAARARAANRARA ANA AKAN AARMAMAAKRA MANA MAAALN NAMA NAN MAA MANA NAMA MAN MANNRLRA KA AANA MANALI ALAA AALAND 


*% lengths of 0 and 13"; 


KRARARAARAKRAAARAAAAAAAAAKRAABRAANA AAA BRAARARMAK AMAA AARABRAAARAAKABRARA MABE AAA AMAA AMAA AAMAS 


% 
ela 
hela GFE 
ME ADLGer Ina lis=—17 
Gas % mumber of increments for incidence angle ####### 
eas % increment step size ######t 
initialincidence=-2; % in degrees HHHHH HH 


1c=(initialincidence+i-1)*pi/180; 
incidencedeg (1)=1c*180/pi; % The inciderce for this iteration, in deg 
% Find slope of Qn vs. alpha curve: 
1lc=Xacc- (XcgTOTALhbar15) *char; % Mament ann from total cg to canard 
% aero center (ft) 
Sells eel sy 1 ale a. % ccl5 is the canard chord length @ 15% SM (in inches) 
%$ bc is the canard span in ft (see “andersonlongxplot") 
% Sc 1S canard area in ft*2 
Vi=le*Sc157 iepar*S)- 
acsd=ac2d/ (i+ac2zd/ (p1*bc7 (ce15/12) *e)); 
hacvb=0.24;  % The aero center of wing/body non-dim by char 
Qmalpha (i) =awb* ( (XogTOTALbar15-hacwb) -VH*ac3d/awb) ; 
Cmacwh=0.00761; % From vortex lattice 
m0 (1) =<Qmacwb-VH*ac3d* ic; 
ir=50; % mumber of increments for next value: 
alphaABS=linspace (0, 15,m) *pi/180; 
Cmog (i, :)=Am0 (1) +Cmalpha (1) *alphaARs; 
(1) =alphaARS (20) *180/p1; % for labelling the incidence arngle 
V (1) =Cimeg (4,20) > % for labelling the incidence angle 
eni. & tor bleck 
disp (XcgTPOTALbar15-hacwb) 
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RAA AAAAKRAAAAAKRAAAKRAAARAARANSAAANANAALRAAAAAAAKAAAABRALAKAAANAAAAAALAIAA AA LAA AALS AIAN SAA MANA ANAK IRIN IN IS 


PLOT Oncg vs. AQA 
plot (alphaABS* 180/pi, Arcg) 
grid; xlabel (‘Absolute angle of attack, alpha a, deg'); 
ylabel('Pitching moment coefficient about c.g., Org’); 
title('AOA vs. Oreg for varying incidence angles') 
% Label the Oncg vs. AOA lines for each incidence: 
for ilies 
Label=sprintf('ic = %g', incidencedeg(1));text (x(1),y(1), Label) 
end 


a oe 


ow 


% Put important info om the plot: 
Mi=sorinitr( Contig = ty ,cUNF IG) ;text(.7,.9,U1, "sc”) 
L2=sprintf('Ballast = %7',Woallast) ;text(.7,.87,L2,‘sc') 
13=sprintf (‘h-hacwb = %', XcgTOTALbar15-hacwb) ;text(.7, .84,L3,'sc') 
LA=sprintf ('Weotal = %g',WIOTALL5) ; text (.7,.81,L4, 'sc') 


mS —colnb &( (ee (15S) = g'ce15) ; text (.7,..78,L5,'sc') 
L6=sprint£ ('Can. ("fran duct) = %',Canardlocation* 12) ;text(.7, .75,L6, 'sc') 
bv=sorintf£(*CG (an) (.15SM) = %', XcgTOTALbar15*char*12) ; 

EGxME C272) 2747, Sc’) 

% Put lines on the plot for Cruise (CR) amd Power Approach (PA) angle of attack: 
xlineCR=(2.28 2.28]; % 2.28 degrees is the Cruise angle of attack 
xlinePA=[(10.61 10.61); % 10.61 degrees is the PA angle of attack 
scales=axis; % "axis" will retum the axes presently in use 
yline=[{scales(3) scales(4)]; % Same for both CR amd FA 
polyline (xlinecR, yline) ;polyline(xlinePA,yline); % Draws vertical lines at 

% the desired CR and PA MA's 
text (xlineCR(1),yline (1), ‘Desired CR AQA') 
text (xlinePA(1),yline(1), ‘Desired PA AQA') 


else 
disp('Program incidence not run --15% S.M. point did not exist for this') 
disp ('config') 

end % 1f to see if 15% static margin point existed 


um nN RAIS RRAANAAAALN AAR AKRANRAANAAAAAAABALRAALRANALNKAAAAAAARANAARAAAAAAREAALNRALAALLAGA 


% FND PROGRAM “ INCIDENCE" 
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APPENDIXG= 
ASSUMPTIONS USED IN SIMULATION PROGRAM 


A. GENERAL 
This appendix contains many of the assumptions made during 

formulation of the AWA-series 3 DOF simulation programs. The 
format of this appendix is informal, stating the assumption, it's 
location in the program, justifying information and, in some cases, 
suggestions for further refining the assumption. The order of the 
assumptions roughly follows the order reached in the program 
AWAI. 


B. ASSUMPTIONS 


(1) ASSUMPTION: Wing area=29.2 ft42 

LOCATION: archytas3 function 

JUSTIFICATION: From vortex lattice, Aquila documentation and 
actual measurement of the wing. 

SUGGESTIONS FOR REFINEMENT: The duct will certainly add to 
the effective wing area--ducts are well known for their capability to 
provide lift. A study of duct aero might provide a better estimate for 
S: 
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(2) ASSUMPTION: Cdo=.06 

LOCATION: archytas function 

JUSTIFICATION: Prof Howard's estimate. 

SUGGESTIONS FOR REFINEMENT: Should account for fan, shroud, 


etc. Wind tunnel test 


(3) ASSUMPTION: Iyy=7.4 slugs-ft*2 

LOCATION: archytas function 

JUSTIFICATION: Aquila's was 10.4 according to ref 8. Reduced it 
25% for lighter archytas 

SUGGESTIONS FOR REFINEMENT: calculate using swing/time 


period method 


(4) ASSUMPTION: Cmgq=-1.47 

LOCATION:archytas function, used in control function for pitch 
damping 

JUSTIFICATION: This is Aquilas Cmq, as per Aquila documentation. 
SUGGESTIONS FOR REFINEMENT: Vortex lattice program, using 


canard and shroud. 


(S) ASSUMPTION: Canard and vane deflections are step inputs, no 
oscillations, no time constant. 

LOCATION: user input 

JUSTIFICATION: probably not true. Depends on servos. 


181 Appendix E 


SUGGESTIONS FOR REFINEMENT: once flight quality servos are 
obtained they could be tested under load to determine fidelity and 
time constant. May need to be modeled as a noise vector in the state 
equations. Note: LCDR Moran and CAPT Kuechenmeister have 


modeled the servo response for their thesis work. 


(6) ASSUMPTION: Sea level, standard day density 

LOCATION: main program 

JUSTIFICATION: simple, probably good for Monterey operations. 
SUGGESTIONS FOR REFINEMENT: could modify program for temp 
affects 


(7) ASSUMPTION: Aircraft coefficient of lift beyond stall as per 
profile discussed in Section VI. 

LOCATION: main program, after "Ydot vector initialization" block 
and again in iteration. 

JUSTIFICATION: 

SUGGESTIONS FOR REFINEMENT: I could find no high alpha 
information during my literature search for similar low sweep angle 
configurations. Wind tunnel tests would further refine this or a more 


thorough lit search. 


(8) ASSUMPTION: Drag model at high alpha is an assumed profile. 
LOCATION: main program 
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JUSTIFICATION: Drag varying as the square of alpha matches well 
with the pre-stall region and seemed reasonable for the post stall 
region (i.e. drag kept going up, quickly, past stall). 

SUGGESTIONS FOR REFINEMENT: Wind tunnel/lit search/CFD. 


Account for entire vehicle. 


(9) ASSUMPTION: Canard pitch authority is simply based on a Cl 
vs alpha (using overall vehicle alpha and an extrapolated NACA OOO9 
airfoil times dynamic pressure (same as vehicle's) times distance to 
c.g. 

LOCATION:control function 

JUSTIFICATION: Good first approximation. No downwash to wing 
from canard is probably a good approximation in that the duct may 
ingest the wake of the canard. 

SUGGESTIONS FOR REFINEMENT: Modify airfoil characteristics 
once configuration is known. Wind tunnel force measurement for Cl 
vs. alpha and flow visualization for canard--wing downwash (and 


flow energization) effects. Lit search? 


(10) ASSUMPTION: No contribution to pitching moment from alpha 
dot or q dot 

LOCATION: control function 

JUSTIFICATION: neglected. Is often neglected in conventionally 


configured aircraft where the effect would be more expected (i.e. 
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time lag of downwash transport) so I neglected it here, feeling that 
the canard's contribution to downwash over the wing would be small. 


SUGGESTIONS FOR REFINEMENT: 


(11) ASSUMPTION: Gyroscopic effects not considered 

LOCATION: 

JUSTIFICATION: 

SUGGESTIONS FOR REFINEMENT: Model them. CAPT 


Kuechenmeister has done this in his thesis work. 
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PROGRAM AWA 
LCDR R.B. Stoney 
Naval Postgraduate School 


Call function "Archytas3" 
to retrieve constants 












Solicit user for: 


# iteration 

DELTAT sO 
Weight Default 
Control Mode Vailable 





Suspend? 







SUSPEND? 
(for user to make 
control changes) 


Yes 


* 










Control Mode 
(options: 1=fixed q 
2=vane & canard 


=] 


Next page 


No pass 


S,K,Cdo,lyy,Lc,Lv,Lw, 
Sc,Sv,cbar,Cmq 






How many secs | 
between suspensions? 
Read: "UPDATE" 






Solicit: 
Canard deflection 
Vane deflection 
Throttle setting 


Next page 








APPEND EXE 





Read: desired q 





Call function Throttle”, which 
converts from user inputed 
throttle position to THRUST 





qdot=0 
Canard and vane 
deflection=0 





q=0 
qdot=0 
(initialize) 





g=32.174 

M=Weight/g 

PT=1 (counter initialized) 

UPDATECOUNT=1 (for program 
suspension to make inputs) 








Solicit initial position/attitude 
height 
speed 
pitch angle 
flight path angle 






x=0 (initial downrange posit) 
"Y" vector initialized: 
Y=[x,h,V,gamma ] 


YDOT VECTOR INITIALIZED: 

Ydot(1)=xdot,speed in x direction=eq'n VII.A 

Ydot(2)=rate of climb=eq'n VII.B 

Ydot(3 )=accel=0 (initialized) t 
Ydot(4)=time deriv of flight path angle=0O (initial) 











ALPHA 
row=density=0.00237688 (SSL) 
Q=dynamic pressure (CAUTION: "q" is pitch rate) :: A 
=0.5*row*Y(3)A2 look Van 
(if/then/else) 


ALPHA=angle of attack 
=theta-gamma=THETA-Y(4) 





Next page 
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"Look up” 
Drag coefficient 
(if/then/else) 


“Dp 


ALPHA 









WRME TO OUMRPUT FILES: 





OUTFILE(1,:)=[x,h, V, Vdot, THETA,ALPHA,counter #,T,xdot,hdot 


gamma,Energy height, energy height dot,q,qdot,canard defl, 
vane defl,time] 





EXTRA(1,:)=[counter#,gammadot,Q,row,Cl,Cd,x/h/V/gamma 
at half step,Mcgw,Mcgc,Mcgv,Mq,L, terma,termb,Q,Cl] 
**This is an output file used for troubleshooting 


STOPFLAG=1 
(this is a counter to end the program when 
the user-inputed # iterations is reached) 










While loop 


8 > 


While STOPFLAG~=999 


STOPFLAG=999 A> 


STOPFLAG~=999 









PT=PT+1 (counter) 


UPDATECOUNT=UPDATECOUNT=1 (suspend option) 


Next page 
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Pass in: 
Current Y(=[x,h,V,gamma] 
THETA,T,S,M,Q,CI,Cd, ALPHA 














CALL FUNCTION "YPRIMES" 





“yprimes" 
equations 
VIILA,B,C,G 









Passed back: 
Ydot(=[xdot,ydot, 
Vdot,gammadot ] 








Half step forward 
differencing (equation 
VII.XX) 






CALL FUNCTION "YPRIMES" 


Predictor, corrector 
and Richardson extrapolation 


(equation VII.XX) 





As before, write to output files: 
OUTFILE(PT,:) 
and 
EXTRA(PT,:) 















n<=0 ? 
(i.e. the vehicle 
hit the ground) 






We PRINT: "Ground impact" 


While loop 







STOPFLAG=999 
Next page (will kick us out 


of while loop) 
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Pass in: q,\DELTAT, THETA, T,ALPHA, {Y}, 
lyy,Q,S,Sc,Sv,cbar,Lc,Lv,Lw,row,Delc,De 


lv,Cma 
Call function "control3” Function "control3" 
Pass back: (see next page) 





qdot,q, THETA 


Yes 
PRINT: “iterations done” 


While loop 


No STOPFLAG=999 BD 


Write to “EXTRA” 





#Hiterations 
reached? 










Time to suspend? 
(to allow user to 
view data and/or 
make control inputs) 


" Plot for user to see 


vehicle flight path 
& 
- 
X 
Yes 


Want a plot? 
Offer a choice of plots 
















Want to change 
Delc, Delv or T? 





Input new Delc,Delv,T 






Whilé€ loop 
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Function ' [Function “controi3” } 
eee) 


THETA=THETA+q*DELTAT | \“ 





a 
Control mode 


2 (canard and vane 


control) 











Wing Moment 


Equation Vll.Ja 






Canard Moment 


Equation VIl.Jb 





Vane Moment 


Equation IV.E 


Pitch damping 


EQUATION VII. | 


and 
q=q+qdot*DELTAT 
THETA=THETA+q*DELTAT 


ENG ENOL Pass: THETA,q,qdot 


Equations VII.XX and 


Vil. Jd 





TO AWA1 


EAC APPENDIX F 











Offer choice of plots. 
Note: the user can ask for 
a plot of any parameter 

in the output file "OUTFILE” 
versus any other parameter 
in the same file 


After running AWA1, the user may then run: 


"Plotter"--gives 8 plots of the most important info 


“Movie"--replays the run 
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cr ce er cc cr cr ep cc rr re rm tr cr cs cr rs cr cr cr rr te ce cr ee ee es ce ee es ce ee ee ee ee ee ee ee ee ee ee 


SL AE SY SSS RR hm emmy Sm me cm mec me co cr ce ee ee ee ee ee ee es 


% LCDR Bob Stoney, Naval Postgraduate School, Monterey 

% Written for APPLE MACINTOSH ® 

+ Last update: 12 MAY 93 (Throttle model and vane model updated) 

% 29 May 93 (PROFILE question added at beginning of program 

% to optimize continuously updated plot for different 
2 flight profiles (hover versus horizontal flight) 


+ Three degree of freedom (pitch, x, z) simulation for ARCHYTAS Vertical 

s Attitude Takeoff and Landing (VATOL) UAV. Includes improvements over PMS3, 

in that the equations of motion (contained in "yprime4" and "control3" do not 
have any small angle assumptions. Also incorporates autoscaling plot feature 
for continously updated downrange distance vs. height plot. 

Thanks to Colin Cooper (aka "Mr MATLAB") of the EC department for help with 
making the graphics look nice. 

Note: "*******" means an area open to further programming (may be 
accompanied by amplifying information concerning work to be done) 


oP 3° lO 


ae 


HN oP 0 


fel 


clear % Clear all variables fram system--wipe the slate clean 
ere % Clear the carmand window 
clg % Clear the graph window 


hold off % To release any previous graph settings 


a RS Sm cc em ee te ee ee ee ee ee ee ee 


PROGRAM SCHEME 

The algorithm is performed as follows: (line 30) 
¢ Aircraft constants--gathers needed constants for particular aircraft 

fram a function. 

¢ Starting values—initializes model's position, orientation and velocity 
¢ Iteration—Utilizes predictor and corrector to increment desired 

% quantities. 

Predictor--first guess at the solution of the governing Dif FO's 

: Corrector--Corrected solutions to the governing Dif EO's 


HNO 0 lO 


yO 0 


ra) 


HO 


2 ¢ Output file management--sets up a matrix of the most important 
2 quantities, which are recorded at each iteration. 

% ¢ Aircraft control--Uses a function to control system inputs. 

z ¢ Final output--To let the user extract desired information. 
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* Functions—This program utilizes the following created functions: 
edrenyEass —wircratt data for the archytas UAV 
"yprimes4"--calculates y' solutions to coupled system equations 
"control 3"--provides aircraft control inputs 
"labeler"--provides labels for the plots 
RODE SS —Converes User imput Of throttle position {in %) to 

Thrust (in lbf) based on engine test data (see "TP vs. T") 
* Other programs to help plot inputs and response: 
"plotter"--plots control inputs and vehicle's response 
"Mplot"—used for troubleshooting. Shows moments into conservation 


AP g\? 


oP of 39 oP ne «0 


eles 


VO 


% of angular mamentum equation contained in "control4" 
% "Movie"--replays the vehicle's position and pitch 
% attitude during the most recent run or use "load workspace" 


mn a em rm me rm i a em i i SS Sra 


Ll=1;L2=1; L3=1;L4=1;L5=1; % Flags for updated plotting routine 

a a> oS PISO CNESAEIC, ASICS EGS i as cs 
{S,K, CDo, lyy, Lc, Lv, Lw, Sc, Sv, Cobar, nq] =archytas3; 

feo Prenvids function fOr Gescription of call variables. (line 60) 


aso ( ") 

disp('What type of maneuver? l=takeoff-transition or hover; 2=level flight") 
disp('or horizontal flight to landing transition') 

PROFILE=input (' '); , , 
a. SSS Nf) SUSE Sn 
disp(' i) 

disp('Input as a row vector: (#iterations,delta t, Weight (lbf),'); 
disp('Control mode (l=fixed pitch rate or 2=cont canard/vane),and'); 

disp ("suspend capability (l=on, O=off)]. Should be five inputs."); 

disp('OR TYPE 1 FOR DEFAULT VALUES=(500 0.05 100 2 1]') 

INITIAL=input (' '); 


if INITIAL—l, % DEFAULT VALUES: 


ITERS=500; SNo. of iterations 
DELTAT=0.05; Time increment between calculations 
Wt=100; tAircraft weight in lbf 
CONTMODE=2; Control Mode (l=const q; 2=user inputs Del cév) 
SUSPEND=1 ; sSuspend flag (to view/control during run) 
else, % (User specified values:) 
ITERS=INITIAL (1); %No. of iterations 
DELTAT=INITIAL (2) ; Time increment between calculations 
Wt=INITIAL (3) ; SAircraft weight in lbf 
CONIMODE=INITIAL (4) ; %Control Mode (l=constant. q; 2=user inputs Del c&v) 
SUSPEND=INITIAL (5) ; Suspend flag (to view/control during run) 


end 
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i Block to set flag ("UPDATE") to alert program to periodically 


Jeo a = halt program for view/control....or will allow uninterrmpted run 
if SUSPEND=1, % (line 90) 
disp('Input the #secs between suspensions'); 
UPDATE=input (' '); % Time period between view/control 
elseif SUSPEND—0, 
UPDATE=0; 
else 
disp('You did not enter a 0 or 1 for suspend option—default is 0') 
UPDATE=0; 
end 
a ae End UPDATE set. logp-===—-—-——-- =a a 
ona Control input DlOCK=s<<=<——====3——==—— > = ae (line sie 
if CONIMODE=1, * Constant pitch rate mode (mostly for troubleshooting) 
disp('(Input the desired pitch rate in deg/sec (down is neg) and TP]") 
qdea=input (' '); 
q=qdeg (1) *pi/180; % Pitch rate converted to rad/sec 
adot=0; % Because constant pitch rate 
Delc=0; % Canard angle set to O for output file 
Delv=0; % Vane angle set to 0 for output file 
TP=adeg (2); % Throttle position in percent 
else, 


if CONIMODE~=2, 
disp('Your input for control mode was not 1 or 2. =detaulU ioe se 
else 
disp('Input as row vector initial canard deflection (LE dcwn is +),') 
disp('Vane deflection (TE down is +) and Throttle position ') 
disp(' (percent—0O is idle 100 is full throttle): [Delc, Delve 


disp("OR TYPE 1 FOR DEFAULT: [-2Z- 0100) % see default loop below 
CONTINIT=input (' '); 
if CONTINIT=1, 
Delc=-2*pi/180; % Default canard deflection in rad(line 120) 
Delv=0; % Default vane deflection in rad 
TP=100; % Default throttle position in percent 
else 


Delc=CONTINIT(1)*pi/180; % Canard deflection converted to rad 
Delv=CONTINIT (2)*pi/180;  % Vane deflection converted to rad 
TP=CONTINIT (3) ; % Throttle position im perecen= 
end SCONTINIT=1? 
q=0; % Pitch rate set to 0 for first entry in output file 
qdot=0; % Pitch accel set to 0 for first entry in output file 
end SCONIMODE~=2? 
end %CONIMODE if,else 
[T}=THROTTLE (TP) ; % call function that converts Throttle position to thrust 


2 


-—--~—-—--———Find Control input block—----—--------—----—---—--------- 
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q=32.174; 

tEWt/g; 2Mass in slugs 

Jeli % Counter 

UPDATECOUNT=1 ; % Counter for program's option to suspend run for view/input 
ROW=0 . 00237688; % Density in slugs/ft*3 


*-——---The Y vector (camprised of 4 elements that constitute the dependent 
Svariables in the differential equations) is as follows: 
% Y(1)=x position in ft, Y(2)=height in ft, Y(3)=Velocity ft/sec, 
5% Y(4)=gamma (flight path angle in radians). Or: [Y]=[x,h,V, gamma] Transpose 
% This vector will be overwritten throughout the program, with required values 
% (line 145) 
% being retained in the OUTFILE matrix (to be introduced below). 
disp (' m) 
disp('Input as row vector (deg) (hi,Vi, gamma (i)--deg, and theta(i)--deg]}') 
disp('Note: Vi must not be zero or a DIVIDE BY ZERO error will occur.') 
eso (OR TYPE I FOR DEFAULT: [15 .5 89.5 90)]") 
Petes 2=inout (' '); 
)i=0: SiMe x postion iy ft 
if INITIAL2=1, % i.e. default selected by user. 
PNOGPTALZ (1)=15; Y(2)=INITIAL2Z (1); % These are all spelled out 
PNET IALZ (Z)=0.5; Y (3)=INITIALZ (2) ; % because they are used in a 
INITIAL2 (3)=89.5; Y (4)=INITIAL2 (3)*pi/180; % configuration definition later 
INITIAL2 (4)=90; | THETA=INITIALZ2 (4) *pi/180; 


else 
Y (2)=INITIALZ2 (1); % Initial height in ft 
i o)— INET LALZ (2) > Teimielel velocity im ft/sec 


cee EAE (3)>89.5 & INETIALZ (3)<90.5, 
INITIAL2 (3)=89.5; % to prevent instability in numerical scheme due to 
% Singularity in gamma dot equation contained in function "yprimes" 
end 
Y (4) =INITIAL2 (3) *pi/180; % Initial flight path angle in rad 
THETA=INITIALZ2 (4) *pi/180; % Initial Pitch attitude in rad 
end % Default if loop 
eS Siem COtOme nae cel Ola os 


tna ech On tial i zZation——_——--—_———=——_---—----—_--—- 

%——---The YDOT vector is as follows: 

% YDOT(1)=Groundspeed in ft/s, YDOT(2)=Rate of climb in ft/s, 

6 YDOT(3)=Acceleration along velocity vector in ft/s*2 

+ YDOT(4)=time rate of change of flight path angle in radians/sec 

% Or; [YDOT]=[Xdot, Hdot, Vdot, gamma dot] transpose 

= This vector will be overwritten throughout the program, with required values 
@)Seing retained in the CUTFILE matrix (to be introduced below). 


mor (1 )=Y (3) *cos(Y (Qlae Initial groundspeed in ft/s 
fer (2) =Y (3) *sin(y (4)); = Initial rate of climb in. ft/s 
YDOT (3) =0.0; % Initial acceleration along vel vector in ft/s*2 
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NN: Ole % Initial flt path angle rate in rad/sec 


% Record configuration: 
CONFIG=(Wt, DELTAT, INITIAL2 (1) , INITIAL2 (2) , INITIAL2 (3) , INITIAL2 (4) }; 


0-052 ROW ree % Dynamic pressure 
ALPHA=THETA-Y (4); %AQA in radians 
5 From "Final Archytas Cl vs. AOA" plot: 
if (ALPHA<=-1.57) | (ALPHA>=1.57), 
CL=0; 
disp('Warning: WING aoa greater than +/- 90 deg') 
elseif ALPHA>-1.57 & ALPHA<-0.72, 
CL=-1 .083-0.69*ALPHA; 
elseif ALPHA>=-0.72 & ALPHA<=-0.285, 
CL=2.81+33.21*ALPHA+105. 38*ALPHA*2+1 38 .1*ALPHA*3+64. 9*ALPHA*4; 
elseif ALPHA>=-0.285 & ALPHA<=0.29, 
CL=9.52e-2+3.38* ALPHA; 
elseif ALPHA>0O.29 & ALPHA<0.9, 
CL=-2.0374+25.54*ALPHA-72.15*ALPHA*2+82 .59*ALPHA*3-33.4*ALPHA%4; 
elseif ALPHA>=0.9 & ALPHA<1.57, 
CL=1 .875-1.194*ALPHA; 
end 


CD=0 .0608+0.0543*ALPHA+0.910*ALPHA%2; % Coefficient of drag as fcn of aoa 
* for post-stall drag increases which CD=f(CL*2) doesn't account for. 

ES=Y (3) 2/(2"g) 42); % Energy neioht 

ESCOt=Y (3) * Yor (3) 7 ot wmeer()- % Time rate of change of specific eneraqy—for 


Cae ie eS a —-BEGIN FILE MANAGEMENT --—-——————- nee 
The outfile consists of the following: 

OQUIF ILE (row#, :)=[{X,H,V, Vdot, Theta, Alpha, Counter#, Thrust, Groundspeed, Rate 

of climb, flight path angle, energy height, time rate of change 
3 of Eng Ht,pitch rate,pitch accel, canard defl,vane cefl, time, TP} 
OUTFILE (1, :)={¥ (1), Y (2), Y(3) , YOOT (3) , THETA*180/p1, ALPHA 130/ pie 
YDOT (1), YDOT (2) , ¥ (4) *180/pi, Es, Esdot, q*180/pi, q@aot*180/p1, Delc =e G7 aie 
Delv*180/pi, DELTAT*PT, TP) ; % Note continuation marks (...) 
EXTRA (PT, :)=(PT, YDOI (4),0, ROW, CL,CD,0,0,07070, 00 One 
% Troubleshooting info (0's for same items ‘cause undefined so far) 

XHPLOTCOUNTER=0; % Used in continuously update x-h plot 


WO oP oO 


> oi? 


LOG APPENDIX=S 


fo] 


STOPFLAG=1 ; 


while STOPFLAG~=999, % i.e. terminates loop when STOPFLAG=999 (assigned) 
Pu=Pitl % Increment counter 
UPDATECOUNT=UPDATECOUNT +1 ; % Increment "suspend" option counter 
XHPLOTCOUNTER=XHPLOTCOUNTER+ 1 ; % Used in continuously update »-h plot 
ALPHA=THETA-Y (4) ; SAQA in radians 


* Predictor and intermediate values 
Peet = Jorimesd(7, TRE IA, &,S,M,O,CL, CD, ALPHA); 3 Call to the D.Ecn Function 


nore): 4, 
i) = (DEL TAr~YDOT (1); SEuler first step forward predictor 
YHALF (1) =Y (1) +DELTAT/2* YDOT (i) ; % Euler half step forward predictor 
end send for 


OF 1=1:4, 
Za) HALE (1) +DELITAT/2*YDOT (1)-> Richardson extrapolation 
Y (1) =2* Y2 (1) -Y1 (i); ‘Correction made 
end send for 
Me (3) 27 (2°o) +Y (2); SEnergy height 
Esdot=Y (3) *YDOT (3) /g+YDOT (2); %New TRC of energy height 
O=0. D* ROW*Y (3) “2; % Dynamic pressure 
ALPHA=THETA~Y (4) ; SAQA in radians 


% From "Final Archytas Cl vs. AQA" plot: 
if (ALPHA<=-1.57) | (ALPHA>=1.57), 
CL=0; 
disp('Warning: WING aoa exceeds +/- 90 deg') 
eceir ALPHA>—1.57 & ALPHA<-—0.72, 
CL=-1.083-0.69*ALPHA; 
elseif ALPHA>=-0.72 & ALPHA<=-0.285, 
CL=2.81+33.21*ALPHA+105. 38*ALPHA*2+138.1*ALPHA*3+ 64. 9*ALPHA” 4; 
elseif ALPHA>=-0.285 & ALPHA<=0.29, 
CL=9.52e-2+3.38*ALPHA; 
elseif ALPHA>O.29 & ALPHA<0.9, 
CL=-2.037+25.54*ALPHA-72.15*ALPHA*2+82 . 59* ALPHA” 3-33. 4*ALPHA* 4; 
elseif ALPHA>=0.9 & ALPHA<1.57, 
Cl] 875-1 -194* ALPHA; 
end 


CD=0 .0608+0 .0543*ALPHA+0. 910*ALPHA’2; % Coefficient of drag as fen of aoa 
%* for post-stall drag increases which CD=f(CL*2) doesn't account for. 
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Cease = —Output file update (add latest iterat to”) ae 


OQUTFILE (PT, :)=[Y (1), ¥ (2), Ys) ae) , THETA*180/pi, ALPHA*180/p1, Pl 1) eee 
YDOT (1) , YDOT (2), ¥ (4) *180/pi, Es, Esdot, q-180/p1, adot* 1607p bas eee 
Delv*180/pi, DELTAT*PT, TP] ; 


32 ee Aircraft control—----~—~--—---—----—---—_-=== 
% ***Constant q or Delc/Delv control, depending on CONIMODE selection*** 
if Y(2)}<=0 %Stops @ ground impact 
disp('Ground Impact--and that is bad!') 
STOPFLAG=999; %Stops while loop 
end send if 


[THETA, g, qdot, Mcgw, Mcgc, Mcgv, Mq, L, Q, CL] =control4 (q, DELTAT, ... 
THETA, T, ALPHA, Y, Iyy,Q, S, Sc, Sv, Coar, Lc, Lv, Lw, ROW, Delc, Delv, Cng, CONT aaa 
=Ca I SGonerol Tuncelen 


1f PT=ITERS, % l.e. if # of iterations is through 
STOPFLAG=999; % Kicks out of while loop 
disp('All iterations campleted') 

end %end if 


FXTRA (PT, :)=[PT, YDOT (4) ,Q, ROW, CL, CD, YHALE (1), YHALF (2); Y BALE (3 eee 


YHALF (4) ,Meqw, Mege, Magu, Ma, L, 0, CL]; 
% Extra info for troubleshooting 





eee ee ea aa —Continuously updated plot for user referenC@———s2s557 ae een 


i929,0,0.0,0,0.0,0,0,0.0,0,0.0.0.0,0.0,0,0.0.0.0.0,0,0.0,0.0.0,0.0.0.0.0.0,0.0,0,0.0.0,0.0.0,0,0.0,0.0,0.0,0,0.0,0,0,0.0,0,0.0,0.0.0,0,0,0,0,0.0,0,0,0,0,0,0,0.4 
TIME=.05; % Put in here so that the plot wouldn't update every time step...slow 


if PROFILE=1, % i.e. takeoff transition 


This plot places a circle ("polymark") at the vehicle's x/hn position and 


& draws a line ("polyline") at an angle corresponding to a/c's pitch attituges 
if (XHPLOTCOUNTER*DELTAT)>=TIME, % 1.€. Causes this plot to update eve TIME F 
’ j ran) 
if ((Y(1)<100 & ¥(2)<100) & LI=1), % This if loop sets plot axes for update 
hold off | | ; 
V=(-25,100,0,100) ; axis (V); 2 Axis limits (xmin, xmax, ymin, ymax) 
plot (QUIFILE (2,1), GgTPiLE(:,2)7 oe | | 
“label ('Downrange distance, x, ft'):ylabel('Height, y, ft’) 


L1=0; % set the plot flag to zero to prevent redrawing of plot 
hold on | 7 
elseif (((Y(1)>=100 & Y(1)<200) | (Y¥(2)>=100 & ¥ (2)<200))) & L2—1), 
hold ort | | . 
V=(-25, 200, 0, 200), axis( 47 2 Axis limits (xmin, xmax, ymin, ymax) 


plot (OUTFILE(:,1),OUTFELE G2) ;.0" | 
vlabel ('Downrange distance, x, ft');ylabel ("Height, y, 
1.2=0 ; 

hold on 


ne) 
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elseif (((Y(1)>=200 & Y(1)<400) | (Y(2)>=200 & Y(2)<400)) & (L3=1))), 
hold ort 
V=(=25, 400, 0,400]; axis (V); % Axis limits (xmin, »max, ymin, ymax) 
poet eri nine, |), UIPILE(:,2),'o") 
xlabel ('Downrange distance, x, ft');ylabel('Height, y, ft') 


L3=0; 
hold on 

elseif (((Y¥(1)>=400) | (Y¥(2)>=400)) & (L4=1)), 
bola orf 


i275, 1000, 0, 1000); axis (Vv) ; % Axis limits (xmin, xmax, ymin, ymax) 
woe (OUlEIEE (=, 1) ,OUTFILE(:, 2), 0") 
xlabel ("Downrange distance, x, ft');ylabel('Height, y, ft') 


L4=0; 
hold on 
elseif (((Y¥(1)>=1000) | (¥(2)>=1000)) & (L5=1)), 
ola off 
(0 2000,.0,2000); axis(V); % Axis limits (xmin, xmax, ymin, ymax) 


lor (QuGe ihe (o, 1), QUITE ELE (:, Z)y to!) 

xlabel ("'Downrange distance, x, ft');ylabel('Height, y, ft') 

LS=0; 

hold on 
end % end if 

Selou(OUlE ILE (:,1),OUTFILE(:,2)) 

polymark (OUTFILE (PT, 1) ,OUTFILE(PT,2),'0') 

Gx=( (V(2)-V (1)) 750) *cos (THETA) ; dy=( (V (4) -V (3) ) /50) *sin (THETA) ; 

polyline ({OUTFILE (PT, 1)-dx OUTFILE (PT, 1)+dx], (OUTFILE(PT,2)-dy ... 

QUIET UE (Et 2) +dy],"=w ) % This shows the vehicle's pitch attitude 

XHPLOTCOUNTER=0; % To reset 


end % end if for "(XHPLOTCOUNTER*DELTAT) >=0.1" 
a PreellE Z25(em Ze OG Lancing transition) SECTION OF CONTINUOUS PLOT-—-----— 
elseif PROFILE=2, % i.e. horizontal flight optimized display 
if (XHPLOTCOUNTER*DELTAT)>=TIME, % i.e. causes this plot to update every TIME 
re (3)>=40, % Velocity > 40? If so, leaves "big picture plot", otherwise 
% focuses on vehicle for control to hover.... 

if ((Y¥(1)<1000 & Y(2)<1000) & Li=1) 

mola off 

| —25, |000,0, 1000] ;axis(V); % Axis limits (xmin, xmax, ymin, ymax) 

BloL (@WIEILE(:, 1), OUIFILE(;:,2),'o') 

xlabel ('Downrange distance, x, ft');ylabel('Height, y, ft') 

L1=0; %* set the plot flag to zero to prevent redrawing of plot 


told on 


ibe), SPEEND I 2G 


elseif (((Y¥(1)>=1000 & Y(1)<2000) | (Y(2)>=10000e ry (2)-2000 ee 
held oft 
V=([—25, 2000, 0, 2000] axis) % Axis limits (»min, xmax, ymin, ymax) 
plot (OUTFILE (: 7 1) Cure IEE (2) 6) 
xlabel ("Downrange distance, x, ft"): ylabel(feighn eee 
L2=0 ; 
DOlewen 
elseif (((Y(1)>=2000 & Y(1)<4000) | (Y¥(2)>=2000 & Y(2)<4000)) & (L3=1)), 
heleeors 
V=[-25, 4000, 0,4000]; axis(V); % Axis limits (xmin, xmax, ymin, ymax) 
plot (QUTPILE (71) CULE IEE eZ) el 
xlabel ("Downrange distance, x, £t');ylabel("Heighr, =e) 
L3=0; 
hold on 
for this if: "if ((Y(1)<1000 & Y(2)<1000) & Ll=1)" 
if (Y(3)<40) & (L4=1), 
Soe Ore 
V=[Y (1) -100, y (1) +100,0,1.5*Y(2)); axis(V); % x +/— 100 and 0 to Ii sheovame 
plot (QUIPILE(:, 1) Ure ih (+715 46") 
xlabel ("Downrange distance, x, ft');ylabel('Height, y, ft') 
L4=0; 
hold: on 
end % end if for Airspeed > 40? (has elseif too) 
polymark (OUTP ILE (PT, 1), OUTEILE (ET, 2)3 or) 
d= ( (V(2)-V (1) ) /50) *cos (THETA) ; dy=( (V (4) -V (3) ) /50) *sin (THETA); 
polyline ( [OUTFILE (PT, 1) —dx oes: 1) +c], (QUTPILE(PI 2 Ge 
OUTFILE (PT, 2)+dy}, '-w % This shows the vehicle's oien attitude 
XHPLOTCOUNTER=0 ; % To reset a allow plot to work in another .1 sec 
end %* end if for ™(XHPLOTCOUNTER*DELTAT) >=0.1" 
end % if for which PROFILE 


MD O 
pare 
‘OF 


5-9,0,0,0,0,0. 0.0.0. 0.0.0.0,0.0.0.0.0.0.0.0.0.0.0.0.0.6.0.0.0.0.0.0.0.0.0.0.0,00.0.0.0.0.0.0.0.0.06 0,0,0,0,0.0.0,0,0,0.0,0.0,0,0.0,0,0.0.0.0.0,0,0,0.0.4 
% —-----—-—--— END Continuously updated plot for user referencC--—-——— =a aaa 
9,0,0,0.0,0.0,0.0.0,0.0,0.0.0.0.0.0.0,0,0.0.0.0,0.0,0.0,0,0,0,0.0,0.0,0,0,0,0.0,0,0,0,0,0,0,0,0,0.0,0.0,0.0,0,0,0.0,0.0 9,0, 0.9.0.9,6.0,0.0,. 


aa ae Suspend loop (to view progress/update control)--—--- 
if UPDATECOUNT*DELTAT>=UPDATE, % i.e. "Time to suspend?” 
UPDATEQOUNT=0; reset it for next suspend 
ans=input ('Program suspended, want a plot? l=yes O=no'); 
while ans~=0, 
holdsers; 
axais( (t02- 354) aris: 
Lil=1; L2=1; L3=1; L4=1; 
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disp("Ihnese are the quantities stored:') 


aigp i=] x posit 6: Alpha 11: Flt Path Angle 16: Delc') 
disp('2: height 7: Point# 12: Energy Ht (Es) 17: Delv') 
(speed Se tnrust Poe ESdoOt 18: time’) 


arson 74: Accel 9: Groundspd 14: q Tier | Rosi. ) 
disp( 5: Theta HOV rave Or Cle iSscdot 7) 
% Ask user for desired x and y plot values: 
xanswer=input ('Type #for what you want on x-axis'); 
yanswer=input ('Type #for what you want on y-axis'); 
if xanswer<l | xanswer>19 | yanswer<1 | yanswer>19, % (i.e. misentry) 
xanswer=18; yanswer=18; % just plots time vs. time as error indication 
end 


plot (OUTFILE(: . xanswer) .OUTFILE (:.vanswer) ) 
[XLABEL, YLABEL]=labeler (xanswer, yansSwer) ; 


xlabel (XLABEL) ; ylabel (YLABEL) ; 
title({'Weight: ',num2str (Wt), 'lbf')) 
pause 
% Ask user if they want another plot: 
ans=input (‘Another plot? Q=NO'); 
end % end suspend while loop 
if CONIMODE=2, 
MAT2=([Delc*180/pi, Delv*180/pi, TP]; 
disp('Present values for Delc,Delv (in deg) & TP:'); 
disp (MATZ) 
ans2=input ('Want to change these settings? O=No'); 
if ans2~-0, 
NEWCONT=input ("Enter as row matrix, new [(Delc(deg),Delv (deg) ,TP(%)]}'); 
NEWCONTSize=sSize (NEWCONT) ; 
Check=NEWCONTSi ze (2); 
a£ check~—3, 
disp ('You misentered the new Control--default is previous input' ) 
else, 
Del c=NEWCONT (1) *pi/180; 
Del v=NEWCONT (2) *pi/180; 
TP=NEWCONT (3) ; 


end % this if block prevents program termination by user misentry. 
end tans 2 if 
elseif CONIMODE=1, 
disp('Present value for pitch rate (deg/sec) & TP:') 
disp (q*180/pi) ;disp(T); 
ans3=input (‘Would you like to change? O=no'); 
if ans3~=0, 
NEWOT=input (‘Input new q, in deg/sec (down is neg) & TP as { }") 
G=NEWOT (1) *pi/180; 
TP=NEWOT (2) ; 
end % ans3 if 
end % CONIMODE=2 or 1? 
(T] =THROTTLE (TP) ; % call to function "THROTTLE" 
end % end new control input if loop 
end 3% end suspend if loop 
— Peet ecco oe — 
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end $end the iteration while loop 
% Display configuration to user: 

disp ('Configuration: ') 

fommat camnpact 











disp(' Weight Delta T hi Vi Gammai ‘«‘Thetai~ = ='') 
disp (CONFIG) 
format short 
$-—--__—----—---——-Final Output —---—-------- -—---——-- 
ansl=input (‘Want a plot? l=yes O=no'); Kick out flag 
while ansl~=0, 

hola off; 


axis({1 2 3 4]);axis; 
L1=1; L2=1; L3=1;L4=1; 
disp('These are the quantities stored:') 


disp(’l: X posit” 6: Alpna 11: Flt Path Angle 16: Delc') 
disp(’2:; heigmte 7. Beimey 12: Energy Ht (Es) Il?) Delve 
disp('3: speed 8: Thrust 13%) Esa@oe 18: time') 
disp('4: Accel 9; Groundspa 14. 3caas 19: Thrtl Posua 


disp('5: Theta 10:Rate of Cl 15: qdot *) 

% Ask user for desired x and y plot values 
xanswer=input ("Type #for what you want on x-axis'); 
yanswer=input ('Type #for what you want on y-axis’); 
disp('Note: will have to add own axis labels, etc’) %******at this point 
plot (OUTFILE (:, xanswer) , OUTFILE (:, yanswer) ) 
[XLABEL, YLABEL]=labeler (xanswer, yanswer) ; 
x label (XLABEL) ; ylabel (YLABEL) ; 
title (['Weight: ',num2str (Wt), 'lbf')) 
pause 

% Ask user if they want another plot: 
ansl=input (‘Another plot? O=NO'); 

end % this while loop 


$-—__——__—----Final instructions (help for storing and retrieving Gata) > ae 
disp ('Recammended carmands to save data:') 

disp ('>>filename-OUTFILE; ') 

disp ('>>filenameC=CONMFIG; (to keep a record of the configuration) ') 
disp('>>save filename filename;') 

disp('>>save filenameC filenamec;') 

disp('Then, when you want to use it again:') 

disp('>>load filename or filenamecC (filename is a 18 by #iters matrix for a 
particular set of conditions—specified on filenameC') 

% NOTE: also useful are programs "filesaver” and “resurrect” 


%-—--NOTE: See also "AWA assumptions" (Microsoft word document) 
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x0 ye jo 


met On| 5,4, DO, Lyy, LC, Lv, Lw, Sc, Sv, Coar, OQnq, xcg]=archytas3 (); 

*Aircraft data function *******Note these are all WAGS at this point*** 
% Used by program PMS3 

% LCDR Bob Stoney 

% Last update: 8 May 93 


S29. 2; % Wing Area in ft*2 

e=0.8; AE Geiciency factor 

AR=4.25; % Aspect ratio (vortex lattice output) 

K=1/ (pi*AR*e) ; % Constant for CD calculation (CD=CDo+K (CL“2) ) 

CDo=0. 06; % Rick's guess | . 

Iyy=7.4; % Moment of inertia (slugs-ft*2) Aquila=10.4 rept AD-A068 345 

Lc=4. 66; % Absolute Distance in ft fram canard's aero center to a/c c.g. 

ie) 5; * Distance in ft fram Vane's aero center to a/c c.g. **x*Refine 

Lw=0.758; % Distance in ft from wing's aero ctr to a/c c.g. (Fram program 
Pe olot!) 

Bk wR kk kt kk kK KK CHECK Salen «Ae a ox 29 6 FD De KEK AKARE A AERAA RARE KR ERE S 

Se—2.6; % Area of canard in ft*2 (From program "Xplot") 

Sv=1.111; mepEea ee two (2! !') vanes mm ft°2 (Each is 8’ by 10") 

er—3,33; e Rhoce Cherd in £t 

Gc]. 92; is (hale) felptepge| (28 pM ee 

lambda=Ct/CR; % Taper ratio 

Char=2/3*CR* (1 + lambda*2/ (1+lambda) ); % MAC 

eiye——-1.47; % Pitch damping (1/rad) Ref: Aquila rept AD-A068 345***** 

% END eyes eam 


yO 
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% 
———————— re al rent artny ee EEE eee 
or 


A campanion program to "AWAL", the 3 DOF simulator. 

This function receives the user input throttle position (TP) (in %—Q=idle, 
100=full throttle) and converts it into thrust (T) in lbf, based on data fran 
engine tests on 4/16/93 and shown curve fit on the cricket graph file 

TOEUSe vse le 


0 of of oP of 


% Bob Stoney, Naval Postgraduate School, created: 5/12/93 
% Last update: 5/12/93 


function [T]=THROTTLE (TP) ; 


flag; 
while flag=1, 
if Te>100 | Te<0, 
disp ('THE THROTTLE POSITION YOU INPUT WAS <0 OR >100, RE-INPUT 0-100") 
TP=rEpue( ") 5 
flag=1; 
else 
T=21.215-.40269*TP+.10932*TP*2-2.1907e—-3*TP*3+1 .7515e—-5*TP*4-5.1479e—-8*TP* 5; 
flag=0; % i.e. will kick us out of while loop and back to main prog. 
end % TP out of bounds if loop 
end 





————— ee eG program "THROTTLE" 
————— 
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nw 


function "yprimes4" 


aw 


0 


LCDR Bob Stoney, NPGS. 

Contains the equations defining conservation of linear manentum, 
used in program AWA-series 

Written: Aug 92 

Last update: 4 May 93 


RB) \a 


nO {Sy 


+y,0 


5 Umprovement over previous editions: 

gamma dot equation (YDOT(4)) doesn't have the small angle assumption of 
previous programs. 

; See AWAL for description of this function 


ae 0 


70 


function [YDOT])=yprimes4 (Y, THETA,T,S,M,Q,CL,CD,ALPHA) {passed out) (passed in) 


% The ydot vector is as follows: [xdot, hdot, vdot, gammadot] 


gre 


This function calculates the y solutions to the coupled system equations: 
YDOT (1) =Y (3) *cos (Y (4) ); 
YDOT (2) =Y (3) *sin(Y (4) ); 
uj © nisl Gi ys % Drag 
32.174; % gravitational constant 
YDOT (3) =T* cos (THETA-Y (4) ) /M-D/M-g* sin (Y (4) ); 
L=0* S*CL; JEG ee. 
Wee — (eam ian(y (4) ))/ (M*Y (3) ): 
TERM2=T* (tan (Y (4) ) *cos (ALPHA) +Sin (ALPHA) ) / (M*Y (3) )? 
ePMs——/ (¥ (3) *costy (4) )); 
TERM4=-YDOT (3) *tan (Y (4) ) /Y (3); 
YDOT (4) =TERM1+TERM2+TERM3+TERMG4 ; 


a0 


>. 


Fat 
= 
ay 
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fo} 


eee 


% LCDR Bob Stoney. Last update: 12 May 93. (vane control power re-modeled based 
% on engine test data). 

% This function is for aircraft control. It is essentially a model of Newton's 
% second law applied to conservation of angular momentum around the pitch axis. 
% This function is a campanion to AWAL, a point mass simulator for the Archytas 


% Change over control3: Improved Vane control power mocel 
% and Cl vs. AOA model. 


% CAUTION: "Q" is dynamic pressure. "gq" is pitch rate. 


function[THETAOUT, gout, qdotout, Mcgw, Mcgc, Mcgv, Mg, L,Q,CLJ=control4... 
(q, DELTAT, THETA, T, ALPHA, Y, LYY, 0; 5, 5G, ov, 2 ee 
Char, Lc, Lv, Lw, ROW, Delc, Delv, Qnq, CONIMODE) ; 


% ----Constant pitch rate section: 
if CONIMODE=1, 
THETAOUT=THETA+q*DELTAT; % Note: for CONIMODE=1 q,qdot have already been 


% specified 
qout=q; 
qdotout=0; % Constant pitch rate. 
Mcgw=0; % These effects not included. Passed for troubleshooting. 
Mcgc=0; % These effects not included. Passed for troubleshooting. 
Mcgv=0; % These effects not included. Passed for troubleshooting. 
MoO? % These effects not included. Passed for troubleshooting. 

$-—--— sing Canard and vane contro! . Saas 


else, i.e. if Delc and Delv specified 


bowing Contribucion, 
SGLEL: 
% Fram “Final Archytas Cl vs. AOA" plot: 
if (ALPHA<=>-1.57) | (ALPHA>]=1.57), 
CL=0; 
disp('Warning: WING aoa greater than +/- 90 deg’) 
elseif ALPHA>-1.57 & ALPHA<-0.72, 
CL=-1.083-0.69*ALPHA; 
elseif ALPHA>=-0.72 & ALPHA<=-0.285, 
CL=2.81+33.21*ALPHA+105. 38*ALPHA*2+138.1*ALPHA*3+64.9*ALPHA” 4; 
elseif ALPHA>=-0.285 & ALPHA<=0.29, 
CL=9.52e-2+3.38* ALPHA; 
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elseif ALPHA>0.29 & ALPHA<0.9, 

CL=—2 .037+25. 54*ALPHA-72 .15*ALPHA*2+82. 59*ALPHA*3-33. 4*ALPHA* 4: 
elseif ALPHA>=0.9 & ALPHA<1.57, 

Si-l ble-l. 1 94*ALPHA, 
end 


% Qnac contribution: (Ref: NACA2412 in Anderson.) 
Cmac=-0.05; 


o\e 


Add Lift and Cmac contributions: 


t=O S*CL; 
Mogw= -L*cos (ALPHA)*Lw + OQmac*Q0*S*Char; 


%-——--Canard contribution: %Presently using fit fram main wing to show high aoa. 
% From "Final Archytas Cl vs. AQA" plot: 
if ((ALPHA-Delc)<=-1.57) | ((ALPHA-Delc)>=1.57), 
Clc=0; 
disp('Warning: Canard aoa greater than +/- 90 deg') 
elseif (ALPHA-Delc)>-1.57 & (ALPHA-Delc)<-0.72, 
Clc=-1 .083-0.69* (ALPHA-Delc); 
elseif (ALPHA-Delc)>=-0.72 & (ALPHA-Delc)<=-0.285, 
Clc=2.81+33.21* (ALPHA-Delc) +105. 38* (ALPHA-Delc)*2... 
+138.1* (ALPHA-Delc) *3+64.9* (ALPHA-Delc) *4; 
elseif (ALPHA-Delc)>=-0.285 & (ALPHA-Delc)<=0.29, 
Clco=9.52e-2+3.38* (ALPHA-Delc) ; 
elseif (ALPHA-Delc)>0.29 & (ALPHA-Delc)<0.9, 
Clc=-2 .037+25.54* (ALPHA-~Delc) -72.15* (ALPHA-Delc)*2... 
+82 .59* (ALPHA-Delc) “3-33. 4* (ALPHA-Delc) *4; 
elseif (ALPHA-Delc)>=0.9 & (ALPHA~-Delc)<1.57, 
Clc=1.875-1.194* (ALPHA-Delc) ; 
end senator Clic loop 
Mcegc=0* Sc*Clc*Lce*cos (ALPHA) ; 


$—-—-Vane contribution: 
% Note: "Delv", the vane deflection, is assumed to be 0 at O control 
% power. Any deviation, plus or minus, produces a pitching mament. 
Vanemode1=2; % Used to campare the two different types of 
% vane control power models. 1="old"—-NACAQ0009 data with 
% velocity measurements fram engine tests. 2="new"—moament 
% data from engine tests (only applies for 90-120 lb but left 
% as is for now) converted to pitch direction. 
% NOTE: Vanemodel=2 was used for all runs used for analysis in the thesis 
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if Vanemodel=1, 
if (Delv)<=.56, 
Clv=0; 
disp('Warning: Vane local aoa less than -.56 (low end)') 
elseif ((Delv)>-.56) & ((Delv)<=-.23), 
Clv=1 .0419+22.5* (Delv) +67.754* (Delv) *2+56. 417* (Delv) *3; 
elseif ((Delv) )>—.23 Saaieel) 2) 
Clyv=2.2772e-20+5 .8809* (Delv) ; 
elseif (Delv)>=.2 & (Delv)<.65, 
Clv=-3.2096e-2+12. 433* (Delv) = 
~37.257* (Delv) *2+27.979* (Delv)“3; 
elseif (Delv)>.65, 
Clyv=0; %***Beyond specified alpha there's no lift*** 
disp('Warning:Vane local aoa gtr .65 (high end)') 
end tend if for Gly doe 
* Fran Thrust tests, the following relationships apply to convert 
% fram thrust to rpm and then to velocity at the vane (Vv) (ft/sec): 
RPME) 7/1. 2t76.65/*T-0. 262055 
Vv=36. 47-2 .016e-3*RPM+1. 867e-6*REM*2; % £ (AOA, V2) ***** 
Mogv——Lv*0. S*ROW* Vv 2*Sv* Clas *% Pitching Moment due to vane 
elseif Vanemodel=2, 
Newvanelift=2* (0.267*Delv*180/pi); % must convert to degrees 
% ‘cause that's how curve was fit. 
Mcgv=-Newvanelift*Lv; 
end % if for vanemodel type 
> ae eC Gann: 
Ma-OngO*S* Gar 2/2" i eee 


%-——Total: Pitch acceleration (by conservation of angular momentum): 


qdotout= (Megw+Mege+Megv) /LyytMqrgq; %q dot to send to main program 
qout=q+qdotout*DELTAT; %q to send to main program 


A ee we = 
et a NS eee — a 


THETAOUT=THETA+qout *DELTAT; *pitch angle to send to main program 
end tend if 


% See also: AWA assumptions (Microsoft word document) 
De ee 
eee eee ed 
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Function "labeler" 


WP ae 40 


This function provides axis labelling for AWA series of programs. 
[ep Reboo Stoney,’ Naval Postgraduate schoo! 
- Last update: 12 MAY 93 (addition of throttle position) 


no 5P iP 


function (LABEL, YLABEL}]=labeler (xanswer, yanswer) 

* The pass in list contains xanswer and yanswer, which is what 
% the user asked to be plotted on an x-y plot. The pass out 
variables XLABEL and YLABEL are text strings to be placed on 
the desired plot. 


B) oa) BS 


if xanswer=1, 

XLABEL='Downrange distance, x, ft'; 
elseif xanswer= 2, 

XLABEL='Height, y, ft': 
elseif “answer=3, 

Po be l= Speed, V, ft/sec’; 
elseif xanswer—4, 

XLAREL='Acceleration, Vdot, ft/sec*2'; 
elseif xanswer==5, 

XLABEL="'Pitch attitude, theta, deg’; 
elseif xanswer—=6, 

XLABEL="'Angle of attack, alpha, deg’; 
elseif xanswer=’7, 

MLABEL='Point number of iteration’; 
elseif xanswer=8, 

Peo erie Limurst, 1, Lbf'; 
elseif “answer=39, 

XLABEL='Groundspeed, ft/sec'; 
elseif manswer—10, 

APE Rate of Glimb, ft/sec’: 
elseif xanswer=—l1l, 

XLABEL='"Flt path angle, gamma, deg'; 
elseif manswer=12, 

POFBEL='Energy height, Es, ft’; 
elseif xanswer=13, 

Pa 7PRi="Bs dot, ft/sec’; 
elseif mhanswer==14, 

XLABEL="Pitch rate, q, deg/sec'; 
elseif xanswer—15, 

XLABFRL="Pitch acceleration, qdot, deg/sec*2'; 
elseif xanswer=16, 

XLAREL='Canard deflection (TE down +),Delc, deq'; 
elseif sanswer-=17, 

XIABEl="Vane deflection (TE up +),Delv, deg’; 
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elselt sans.wer—1S5, 
MUABREL=" vine te see | 
elseif xanswer—le, 


XLABEL='"Throttle pesitilOne le, eereeii=a 


end % end if loop for “LABEL 


if yanswer—1, 

YLAREL='Downrange distance, x) fe — 
elseif yanswer==2, 

YIABEL= "He loin eee 
elseif yanswer=3, 

YLABEI='Speed, V, ft/sec’; 


elseif yanswer—4, 

YLABEL='Acceleration, Vdot, ft/sec*2'; 
elseif yanswer=5, 

YLABEL="Pitch attitude, ihetawa@eq = 
elseif yanswer=6, 

YLABEL='Angle of attack, alpha, deg': 
;elseif yanswer=7, 

YLABEL='Point number of iteration'; 
elseif yanswer=8, 

WIABEL— Inurseeal.. ler: 
elseif yanswer==9, 

YLABEL='"Groundspeed, ft/sec'; 
elseif yanswer=10, 

YLABEL='Rate of Climb, ft/sec'; 
elseif yanswer=11, 

YLABEL='Flt path angle, gamma, deg'; 
elseif yanswer—=12, 

YLABEL='Eneragy height, Es, ft'; 
elseif yanswer=13, 

YIABEL="Es dor, ft/sec’: 
elseif yanswer—14, 

YLABEL='Pitch rate, q, deg/sec'; 
elseif yanswer—=15, 


YIABEL='"Pitch acceleration, qdot, deg/sec*2'; 


elseif yanswer=16, 


YLABEL='Canard deflection (TE dewn +),Delc, deg'; 


elseif yanswer=17, 


YLABEL="Vane deflection (TE up +),Delv, deg’; 


elseif yanswer=18, 
YLABEL="Time, €, sec’; 
elseif yanswer—=19, 


YLABEL='"Throttle position, Te aerpercenr = 


end % end if loop for XLABEL 
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end function "labeler" 
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= Program used in conjunction with the AWA series programs. 

% Uses subplot to plot most important parameters on a two 

% sheets of paper. 

% To use, simply type plotter at the MATLAB cammand carrot >> 
PeLciR bes Stoney, NPGS 

= Created 12 SEP 92 


cla 
hold off 
— Aircraft parameters: Alpha, speed, theta, q-—--——-- 


eMeelou(7 1), plot (OQUIPILE(:,18) ,OUTFELE(:, 6) ) 

“label ('Time, t, sec');ylabel('Angle of attack, alpha, deg');grid 
mre et |’ Weiaqht—’, mimzstr (wt), ' lbft *)) 

py OL(277), Plot (OUTT ILE(:,18),OUIFILE(:, 3) ) 

xlabel('Time, t, sec');ylabel ('Speed, V, ft/sec');grid 

Se elee( 273) polot (OUIFILE(:, 18) ,OUIFILE(:,5)) 

xlabel ('Time, t, sec');ylabel ('Pitch attitude, theta, deg');grid 
Sriemiiet (774), plot (OUIFILE(:,18) ,OUITILE(:,14)) 

xlabel ('Time, t, sec');ylabel('Pitch rate, q, deg/sec');grid 


pause 
clg 

hold off 

— Control parameters: Canard, vane, thrust, and x vs. h----—---- 


pee Ou 221), plot (OUIFILE (:,18) ,QUIFILE(:,16) } 

pabel("Yime, t, sec’');ylabel (‘Canard defl., Delc, deg')-grid 
title({'Weight=', num2str (Wt),' lbf')) 

eemolot (222), plot (OUIFILE(;, 18) ,OUTFILE(:,17)) 

xlabel('Time, t, sec');ylabel ('Vane deflection, Delv, deg'):grid 
subplot (223), plot (OUTFILE (:,18), OUTFILE(:,19)) 

mlebel (Time, t, sec’) ;vylabel('Throttle posit, TP, %');grid 
subplot (224), plot (OUTFILE(:,1),OUTFILE(:,2)) 

xlabel ('Downrange distance, x, ft');ylabel('Height, h, ft');grid 
pause 

Isola off 


io 


27) program “plotters 


0 


0 
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——— == SS MISC, program Tove" 


S 
ray 








« 'Movie' LCDR Bob Stoney, NPS. 1 March 93 
* Last update: 30 May 93 

*« Companion program to AWA1 that shows the vehicle's position and pitch 
attitude during the most recent run. 


alg 
hold off 


ansl=size (KEEPMARK) ; 
ansla=ansl (1); 


disp('What axes would you like on the movie? [xmin, xmax, ymin, ymax]') 
disp('Enter 1 for default to last axis used on AWAl run.') 
movieaxis=input(' '); 

if movieaxis=—l1, 


axis (V); % 1.e. DEFAULT 
else, 
axis ([movieaxis (1),movieaxis (2) ,movieaxis (3),movieaxis(4)]):; 
% Axis limits ([xmin, xmax, ymin, ymax)) 
end 


PLOC(OUTE TEL), OUP ELE (i. 2) so.) 

“label ("Downrange distance, x, ft');ylabel('Altitude, h, ft') 
Li=sprint£ (‘Delta T = %q', DELIAl) ; Gexeie2,.. 32,1), sew 
L2=sprintf ("Weight = %o°,Wt)+text (127.8 7, se) 


STEPSIZE=5; % Interval of plotting the vehicle's position and pitch attitude 


for PT=1:STEPSIZE: ansla; 
= Scale the pitch attitude to the last scale used by main program 
xScaler= (Vkeep (ansla, 2) -Vkeep (ansla, 1) ) / (Vkeep (PT, 2) -Vkeep (PT, 1) )? 
%to scale pitch attitude line in x direction to last plot coordinates 
c=KEEPdxdy (PT, 1) *XScaler; 
YScaler= (Vkeep (ansla, 4) -Vkeep (ans1la, 3) ) / (Vkeep (PT, 4) -Vkeep (PT, 3)); 
Sto scale pitch attitude line in y direction to last plot coordinates 
Gy=KEEPdxdy (PT, 2) *YScaler; 
* end scaling 
polymark (OUTFILE (PT,1),OUTRILE (ET 2) es) 
polyline ((QUIFILE (PT, 1)-dk OUTER ILE (eae 
, (OUTFILE (PT,2)=dy CUIPILE (PT, 2)7 a) eae 
% This shows the vehicle's pitch attitude 


———— rig program Wile! 
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a0 


\O 


ee SO, program 4c OGG) ad 


0 


xe 


ae een MOlee— = = 

Program to aid in troubleshooting the control function. 

Simply generates a plot of the Moments about the cg 

due to Wing, canard, vane and pitch damping vs. time. 

Also shows (on another plot) the theta, q and qdot results of 

the control function. 

To use: type "Mplot” at the command carrot. 

LCDR Bob Stoney 

Created 13 SEP 92 

clg 

hold off 

subplot (211) 

eee io (ey oper x RA? 04), '—" ,OUTPILE (:,18),EXTRA(:, 12), *--‘,.... 
Oe wee rota 5) '*' OUIFILE (:,18), EXTRA(:,14),°0O") 

“label ('Time,t,sec') ;ylabel ('MWing- C-—- V* q oo');grid 

eeere (('Weight=",num2str(Wt),' lbf ')) 

subplot (212) 

fe oor tre tS) OULF me (2,5), '=",OUTFILE(:,18) ,OQUTFILE(:,14),'——',.... 
COUPE TE, bS) OUTEILE (+,15),°:") 

Plabel ('Time,t,sec') ;ylabel (‘theta- q-- qdot..');grid 

eicle((‘'Weight=',numZstr (Wt),' lbf')) 

use 


NO NW pO VIO 


yo 


FW go 


\O 


ov? tig 


yo 


Cl program "Mp1 OF Lg 


oe 


Dies APPENDIA G 


oO x0 


OD program “Ti besaver’ 


% "filesaver"...this program is used to save information™fram 
« program AWAL. 

= LCDR R.B. Stoney, Naval Postgraduate School 

z May 1993 


RO V=ONTFILE; 

R2 7E=EXTRA;: 

R2 1C=COF IG; 
save R27 R27; 
save R27E R27E; 
save R27C R27C; 
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He 


Bhd 


et SG program “FeESUrrecc 


ye 


ve 


Beesurrect’....a program to load the results of a run of program 
AWA] that were saved. See also "filesaver" 

ICDR R.B. Stoney, Naval Postgraduate School 

May 1993 

user: change the number on each line for the appropriate file 

+ Note: can use "save workspace" if desired but that takes up a LOT 
of memory 


vO yO 


NO aC 


no NO 


load R24 > “OUTFILE” matrix 

load R24C > "CONFIG" matrix 

load R24E % the "EXTRA" matrix 
OUTFILE=R24 ; 

WL=R24C (1); 

disp('File number R24 has been recovered') 


HO 9 


——————————e rn ond program "TeESUITrecL” 
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Figure H.2 
Archytas simulation: "High vane deflection pushover from hover’ 
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Figure H.4 
Archytas simulation: “Vertical acceleration to gentle pushover’ 


Appendix H 


Pitch attitude, theta, degAngle of attack, alpha, deg 


Canard defl., Delc, deg 


Throttle posit, TP, % 





Speed, V, ft/sec 


6 : 
0 2 10 Js 





Pitch rate, gq, deg/sec 


0 2 10 IES 


Time, t, sec Time, t, sec 


Weight=100 lbf 








0 5 10 IS 


Time, t, sec Time, t, sec 


Vane deflection, Delv, deg 


Height, h, ft 





0 5 10 IS 0 500 1000 
Timertesec Downrange distance, x, ft 
Bigunewiae 
Archytas simulation: maximum effort slow flight, unsuccessful 


attempt 
at level transition from horizontal flight to hover 


Appendix H 


tO 
tO 
= 


Pitch attitude, theta, deg Angle of attack, alpha, deg 


Canard defl., Delc, deg 


Throttle posit, TP, % 





Speed, V, ft/sec 


i@) 

-D) 

7 9) 

“eD 

oD] 

aa) 

=> 

3) 

cs 

= 

<= 

0 5 10 15 & 
Time, t, sec 














oD 
10 Weight=100 Ibf 4 
| 3 
ON eo as hy 2 Me ee je a 
2: aes 
_|@ ee: rere Sale oo = 
| . 2 
YO ———— as o 
0 5 10 IS 3 
hime stusec > 
100 
s * : 
Sy | = 
| oD 
cae 
Q ee a 
0 5 10 15 0 200 400 600 
ime st, SCC Downrange distance, x, ft 
Figure H.6 


Archytas simulation: Transition to horizontal flight 
followed by "zooming" transition back to hover 
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